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1.  Objectives 


There  is  an  increasing  demand  in  today's  electronics  systems  for  greater  complexity 
and  higher  performance.  With  data  rates  increasing  into  the  lOO’s  of  MHz  to  multi-GHz 
range  (X-Band),  issues  arise  such  as  thermal  management,  reliability,  miniaturization,  and 
manufacturability  which  must  be  resolved  through  the  development  of  new  electronic 
packaging  technologies.  A  key  problem  in  such  systems  is  increasing  the  density  of  the 
electronic  components,  not  only  to  reduce  size  and  weight,  but  in  many  cases  to  increase 
performance  as  well.  One  attractive  method  of  accomplishing  this  in  complex  systems  is 
multichip  module  (MCM)  packaging  technology.  MCM  technology  permits  high- 
density  packaging  of  portions  of  high  speed  digital  electronic  systems  and  combined 
analo^digital  mixed-mode  functions  into  a  Mixed  Mode  Module  (MMM).  What  is 
required  is  a  mixed-mode  module  (MMM)  technology  designed  to  provide  a  medium  for 
efficiently  integrating  a  number  of  technologies  in  high  density  modules  with  the  goal  of 
significantly  impacting  identified  high  leverage  products  and  markets,  both  commercial 
and  defense,  by  reducing  size,  weight  and  cost,  and  increasing  performance.  Strategic 
Analysis,  Inc.  (SA)  has  endeavored  to  achieve  this  goal  through  its  Air  Force  program 
“Research  On  Mixed  Mode  Modules”. 

After  careful  evaluation  and  based  on  results  achieved  in  this  program,  SA 
determined  that  additional  work  would  be  required  in  evaluating  3D  interconnects  for 
MMMs.  SA  performed  analysis  on  interconnects  for  various  3D  MMM  architectures, 
and  summarized  the  results  in  Progress  Report  #1:  Research  on  Mixed  Mode  Modules, 
dated  September  1,  1996.  As  a  result  of  this  effort,  SA  determined  that  research  in  this 
topic  would  prove  of  great  value  to  the  Air  Force  in  analyzing  the  various  architectures 
that  are  needed  for  3D  MMM  structures.  As  part  of  Task  9  of  the  contract,  SA 
determined  that  3D  mixed-mode  systems  will  greatly  reduce  size  and  weight  (with 
potential  reductions  of  cost  and  significant  increases  in  performance)  for  a  wide  range  of 
military  and  commercial  electronic  applications. 

SA  also  determined  that  further  reductions  can  be  realized  in  the  size  and  weight 
of  3D  and  2D  MMMs  by  integrating  passive  components  into  the  substrates  of  these 
modules.  In  some  modules  discrete  passive  devices  (resistors,  capacitors,  inductors)  can 
take  up  over  50%  of  the  substrate  area  of  a  module.  Through  elimination  of  solder  joints 
and  reduction  in  the  number  of  assembled  components,  this  will  also  enable  substantial 
improvements  in  the  reliability  of  these  modules  for  the  military  markets.  SA  proposed 
an  additional  task  #10  for  determining  the  viability  of  developing  technologies  and 
materials  that  can  be  used  for  processing  passive  devices  into  the  substrates  of  3D 
MMMs.  This  study  also  evaluated  the  feasibility  of  using  these  processes  in  current  2D 
MMMs.  The  objective  of  the  research  efforts  undertaken  by  SA  in  this  contract  were  to 
aid  in  the  advancement  of  basic  research  in  key  areas  of  electronics  miniaturization  that 
may  prove  crucial  in  meeting  the  challenges  of  future  Air  Force  warfighting. 
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2. 


STATUS  OF  EFFORT 


SA  performed  research  on  Mixed  Mode  Module  (MMM)  technology  spanning  the 
range  from  applications  and  mixed-mode  systems  requirements  to  the  baseline 
technologies  which  must  be  put  into  place  to  realize  MMMs.  SA  examined  promising 
military  applications,  and  analyzed  the  critical  issues  relating  to  the  design,  fabrication, 
and  manufacturing  of  MMMs  for  advanced  electronic  systems.  SA  has  defined  one  such 
example  of  a  system  that  incorporates  both  analog  and  digital  components,  the  Global 
Positioning  Satellite  (GPS)  Receiver.  Strategic  Analysis,  Inc.  also  investigated  the 
requirements  for  design  tools  for  simulating  an  MMM  with  special  attention  paid  to 
insuring  analog  and  digital  signal  integrity  and  optimizing  the  amount  of  silicon  or  other 
device  area  on  a  high  density  MMM  host  substrate.  SA  has  determined  that  significant 
effort  still  remains  to  be  done  in  the  area  of  CAD  development  for  mixed  signal 
technology  so  that  analog  and  digital  circuits  can  be  simulated  at  the  same  time. 

The  ultimate  packaging  density  of  the  MMM  and  its  maximum  size  depend  on  the 
patterning  technology.  SA  has  evaluated  some  of  the  technological  advances  in  cost- 
effective  patterning  and  via  drilling  for  MMM  substrates.  SA  has  also  performed  an 
analysis  of  various  interconnect  strategies  for  determining  the  architecture  that  would  be 
most  beneficial  for  MMM  systems.  One  solution  to  thermal  management  problems 
resulting  from  high  density  MMMs,  which  SA  has  investigated,  is  to  use  diamond 
substrates  which  have  high  thermal  conductivity  capabilities.  While  the  technology  for 
MMMs  is  being  developed,  one  of  the  things  that  cannot  be  overlooked  is  the  integration 
of  MMMs  into  systems.  SA  evaluated  the  requirements  for  an  infrastructure  to  produce 
cost-effective  MMMs.  The  system  house  responsibilities  for  driving  the  MMM 
technology  in  advanced  electronic  systems  was  analyzed.  Mixed  mode  system  packages 
are  being  produced  presently,  but  they  are  not  the  product  of  well-developed  design 
methodology.  Based  on  careful  evaluation  and  the  results  achieved  in  this  program,  SA 
determined  that  research  in  evaluating  3D  interconnects  for  MMMs  would  prove  of  great 
value  to  the  Air  Force  The  principal  problem  that  low-cost  3D  electronics  packaging  will 
address  is  the  critical  need  in  a  wide  range  of  military  systems  to  greatly  reduce  the  size, 
weight,  and  volume  that  various  vehicles  and  platforms  must  carry. 

In  current  practice,  implementing  complete  mixed-mode  systems  or  sub-systems 
requires  mounting  individual  packaged  devices,  MCMs,  shielded  or  unshielded  analog 
modules,  and  numerous  passive  components  for  filters  and  RF  processing  on  a  variety  of 
substrate  types.  This  approach  constrains  the  size,  weight  and  performance  capabilities 
of  defense  applications.  In  order  to  meet  military  requirements,  the  development  of 
materials  and  processes  for  integrating  discrete  passives  (capacitors,  resistors,  and 
inductors)  directly  into  the  MMM  substrate,  is  important  since  embedded  passives 
technology  can  reduce  substrate  area  by  as  much  as  ten  times  while  providing  increased 
performance  and  reliability  at  1/4  of  the  cost.  In  addition,  many  secondary  issues 
associated  with  MMM  such  as  including  cost  effectiveness,  yield,  testability,  and 
reliability  receive  little  attention  during  the  development  stage  of  a  technology,  though 
their  impact  is  as  important  as  first  level  research  questions.  SA  directed  its  research 
effort  at  the  understanding  the  system  impact  of  high  frequency,  embedded  passives  and 
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secondary  issues  associated  with  MMM  technology  earlier  in  the  development  cycle  to 
more  effectively  direct  primary  research  and  understand  applications  limitations  prior  to 
instead  of  after  project  completion.  Efforts  focused  on  experimentation  and 
microstructure  analysis  of  early  material  candidates  for  embedded  passives  for  high 
frequency  applications.  The  results  of  this  research  can  provide  direct  input  into  the 
design,  planning  and  development  of  next  generation  MMM  systems  for  defense 
applications. 
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3. 


ACCOMPLISHMENTS 


3.1.  Task  One:  Identify  and  define  electronic  systems  that  will  benefit  from  Mixed 
Mode  Module  technology. 


The  portable  communications  and  wireless  market  is  taking  advantage  of  Mixed 
Mode  technology  for  increased  functionality  and  performance  while  reducing  size  and 
costs.  The  military  will  take  advantage  of  the  new  technologies  for  such  applications  as 
smart  weapons  and  systems  that  bring  real-time  intelligence  to  the  battlefield.  An 
example  of  an  important  military  application  of  Mixed  Mode  technology  is  the  Global 
Positioning  Systems  (GPS),  a  worldwide  satellite  navigation  system  developed  by  the 
DoD.  GPS  receivers  have  application  to  both  commercial  and  military  systems  and  are 
also  expected  to  take  over  the  civil  aviation  market  in  the  next  ten  years.  A  typical  high 
performance  GPS  system  using  advanced  MMM  technology  measures  only  about  6"  x  3" 
X  9",  consurnes  several  watts  of  power,  and  weighs  several  pounds.  The  GPS  architecture 
allows  for  partitioning  into  multi-use  modules  for  a  variety  of  requirements.  A  modular 
approach  to  MMMs  provides  savings  both  in  NRE  design  costs  and  unit  costs  because 
of  the  large  combined  applications  of  the  modules.  MMMs  also  reduce  the  size  and 
power  of  the  GPS  receiver  to  significantly  increase  the  total  GPS  market  potential.  Using 
a  mixed  mode  analog/digital  module  design,  it  is  possible  to  reduce  the  overall  size  of  the 
current  GPS  until  it  fits  in  a  package  the  size  of  a  \sTistwatch.  In  achieving  this  goal  it  is 
necessary  to  also  incorporate  the  integral  components  such  as  discrete  resistors, 
capacitors  and  inductors  into  the  substrate.  This  technology  is  discussed  further  in  Task 
Ten. 

3.2.  Task  Two:  Characterize  the  types  of  electronic  modules  suitable  for  integration 
on  a  Mixed  Mode  Module.  Evaluate  the  suitability  and  level  of  increased  performance 
with  incorporation  of  optical  interconnects  for  realizing  external  I/O  functions. 


Data  transmission  is  increasingly  the  bottleneck  of  modem  communications.  In 
future  systems,  optical  interconnects  between  and  within  electronic  cabinets  are  expected 
to  be  cost  effective  for  both  commercial  and  military  processors.  Furthermore,  some 
critical  I/O  bottlenecks  in  current  embedded  processors  might  be  eliminated  only  by  using 
optical  interconnects.  Optical  I/O  applications  of  module-level  packaging  include  parallel 
and  serial  cabinet-to-cabinet  intercoimects  for  computer  netw'orking,  workstation-to- 
peripheral  interconnects,  and  commercial  supercomputers.  Optical  intercoimects  using 
ffee-air  or  single-mode  fiber  optical  waveguides  offer  extremely  low  dispersions. 
Interconnect  lengths  are  on  the  order  of  tens  of  kilometers  for  wide  signal  bandwidths. 
With  low  dispersion  graded-index  multimode  fibers,  interconnect  lengths  are  on  the  order 
of  several  hundreds  of  meters,  much  greater  than  with  normal  metal  transmission  lines. 
The  disadvantage  with  this  technology  is  the  requirement  for  transducers  at  both  ends. 
Optical  transmitters  and  receivers  are  generally  expensive  and  power  consuming  and  add 
substantial  signal  propagation  delay  compared  to  their  normal  metal  coimterparts.  The 


4 


relatively  low  signal  amplitudes  from  the  photodetector  in  the  receiver  must  be  amplified 
to  digital  signal  levels  with  transistors  of  limited  gain-bandwidth  product  which  leads  to 
significant  optical  receiver  propagation  delays.  For  intra-module  interconnects  this  is  a 
serious  technical  concern. 

One  factor  that  could  have  a  profound  effect  on  the  implementation  of  board-level 
optical  interconnects  is  the  vertical-cavity  surface-emitting  laser  (VCSEL)  which  has 
demonstrated  significant  performance  gains  and  size  reductions.  Several  VCSEL-based 
transceiver  modules  are  now  being  offered  commercially.  Hewlett-Packard  Co.  has  a 
106.2-Mbit/second  transceiver,  packaged  in  a  1  x  9  configuration,  as  well  as  a  Gigabit 
Ethernet  transceiver.  Vixel  Corp.  is  bringing  out  a  1 -Gbit/s  VCSEL-based  component  for 
fiber-optic  applications;  and  Sun  Microsystems,  using  VCSEL  technology,  has  put 
together  what  it  calls  the  fastest  ATM  networking  board  available  today. 

Early  lateral-emitting  diode  lasers  had  an  awkward,  5-micron  x  1,000-micron 
geometry.  The  latest  VCSELs  have  a  round  aperture.  In  addition,  the  diodes  radiate  light 
vertically,  making  them  convenient  for  packaging  in  array  configurations.  New  devices  are 
being  manufactured  with  greater  than  99  percent  yield,  and  the  newer  oxide-contained 
diodes  have  pushed  the  threshold  current  from  around  1,000  micro  A  to  below  10  micro 
A.  The  current  research  goal  is  to  push  the  diodes  to  4  square  microns  and  develop 
packaging  and  mounting  techniques  to  integrate  them  with  CMOS  VLSI  chips. 

Another  significant  breakthrough  is  the  advent  of  high-performance  plastic  fiber. 
Graded-index  fibers  made  from  perflourinated  polymers  have  realized  a  200-times 
improvement  in  bandwidth  loss  over  previous  plastic-fiber  technology  and  could  cut 
interconnect  costs  by  a  factor  of  10.  The  high-efficiency  fiber,  coupled  with  efficient 
VCSELs,  could  reduce  the  power  cost  of  interconnect  at  the  backplane  from  100  mW  to 
around  1  mW. 

3.3.  Task  Three:  Investigate  the  requirements  for  design  tools  simulating  a  Mixed 
Mode  Module  with  special  attention  paid  to  insuring  analog  and  digital  signal 
integrity. 


Noise  problems  in  MMMs  result  from  the  large  digital  signals  affecting  the  analog 
circuits  due  to  direct  inductive  pickup  from  poorly  shielded  signals,  to  the  high  inductance 
exhibited  by  power  and  ground  feeds,  and  to  the  inductance  in  the  power  and  ground 
planes  from  through  vias,  wire  bonds,  package  leads  and  lead  frames,  and  the  balls  in  ball 
grid  array  packages.  Metal  layers  in  MMMs  are  too  thin  compared  with  their  skin 
depths  to  shield  ac  magnetic  fields  at  the  near  audio  frequencies  found  in  slower  speed 
applications.  Design  tools  that  account  for  these  effects  are  not  available  and  are  essential 
in  the  long  run.  There  has  been  limited  success  in  theoretically  modeling  the 
electromagnetic  behavior  of  power  and  ground  planes.  Diagnostic  test  approaches  to  the 
characterization  of  power  and  ground  plane  noise  are  primitive.  Such  test  capabiliU’  will 
be  necessary  as  an  aid  in  the  current  "build  and  fix"  approaches,  then  to  test  the  emerging 
physical  models,  and  in  the  long  run  to  provide  inputs  to  CAD  development. 
Commercial  CAD  vendors  have  not  yet  focused  on  the  mixed  signal  problem,  due  to  a 
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lack  of  resources  and  the  pressure  of  more  immediate  problems.  These  common  CAD 
and  measurement  needs  may  require  a  cooperative,  coordinated  collaborative  approach. 

3.4.  Task  Four:  Evaluate  alternative  substrate  materials  for  Mixed  Mode  Module 
host  substrate  use,  including  extensions  of  PWB/MCM-L  type  technology  to  the  region 
of  higher  interconnect  and  via  densities  required for  Mixed  Mode  Modules. 


As  MMM  circuit  complexities  increase,  the  need  for  heat  dissipation  becomes 
greater.  MMMs  with  ceramic  substrates  offer  a  high  thermal  performance  but  at  a  high 
cost.  This  technology  is  used  exclusively  for  high  end  niche  markets.  MMMs  with 
flexible  laminate  substrates  are  a  lower  cost  alternative,  but  are  limited  to  low  power 
commercial  applications  due  to  a  high  thermal  resistance.  SA  has  evaluated  several 
technologies  that  offer  a  lower  cost  alternative  to  the  ceramic  substrate  MMM. 

David  Samoff  has  developed  a  double-sided,  low  temperature  co-fired  ceramic  on 
metal  core  technology  (LTCC-M).  The  metal  core  is  copper/moly/copper.  The 
advantages  of  this  technology  are  that  it  has  more  thermal  conductivity  than  a  laminate 
board,  the  metal  core  prevents  xy  shrinkage  of  the  ceramic,  and  fewer  processing  steps  are 
required.  LTCC-M  involves  hours  rather  than  weeks  of  processing  time.  A  custom  glass 
bonding  layer  is  used  to  attach  the  ceramic  to  the  metal  core  during  the  cerarmc  firing 
process  which  prevents  the  ceramic  from  shrinking  in  the  lateral  plane  (<0.5  ±  0.1% 
versus  10-15%).  The  glass-ceramic  has  the  lowest  reported  dielectric  loss  to  40  GHz.  A 
low  cost  method  for  laser  drilling  seven  mil  feedthroughs  in  the  metal  core  is  used.  Silver 
thick  film  conductors  for  vias  and  traces  compatible  with  the  dielectric  system  formulated 
have  been  developed  for  this  technology.  An  ability  to  fabricate  thin  film  wiring 
structures  on  LTCC-M  boards  has  been  demonstrated.  A  cost  analysis  carried  out  for  the 
LTCC-M  technology  shows  it  to  be  cost  competitive  with  MCM-L. 

AlliedSignal  Electronics  Technology  Center  has  developed  a  chip-on-board  MMM 
technology  that  enables  both  high  circuit  density  and  high  power  density  at  a  significantly 
lower  cost  than  ceramic.  The  substrate  is  Viagrad,  a  flexible  polyimide  laminate  2  mils 
thick  which  is  a  product  of  Sheldahl.  Viagrad  contains  laser  ablated,  copper  plated  vias  in 
a  regular  grid  pattern  spaced  every  0.5  mm.  The  vias  are  used  for  both  trace-routing  and 
thermal  conduits  for  heat  to  flow  from  the  circuitry'  on  top  through  the  polyimide  to  a 
heat  sink.  The  via  structure  lowers  the  thermal  resistance  of  the  polyimide  to  as  low  as 
0.05  C°AV/in"  as  compared  to  0.4  C°AV/in^  for  standard  polyimide  materials.  Viagrad  is 
also  capable  of  operating  at  ambient  temperatures  in  excess  of  200°C  without  degrading. 
The  substrate  can  be  attached  directly  to  a  heat  sink  using  a  3  mil  thermoplastic  adhesive. 
The  heat  sink  not  only  serves  as  a  cooling  path  for  the  electronics,  but  also  acts  as  the 
mounting  surface  for  the  flex  laminate  substrate  during  assembly.  A  low  stress  epoxy  is 
used  to  encapuslate  the  electrical  components  and  provides  a  good  moisture  and 
environmental  barrier. 
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3.5.  Task  Five:  Investigate  large  area  patterning  and  via  drilling  approaches  for 
Mixed  Mode  Modules. 

For  production  quantities  of  Multi  chip  Module  (MCM)  substrates,  mechanical 
drilling  of  vias  represents  neither  a  practical  nor  an  economic  option.  Prices  for  drill  bits 
smaller  than  250  pm  (10  mil)  in  diameter  range  up  to  $15  each  -  with  a  realistic  maximum 
number  of  hits  between  500  and  1000.  Hit  rates  of  less  than  one  hole/sec  and  drill 
breakages  translate  into  very  slow  processing  and  low  yields.  Reactive  ion  etching  (RIE) 
presents  a  possible  alternative.  RIE  generates  heat,  however,  which  can  cause  dimensional 
stability  problems,  and  has  limited  spatial  resolution.  Wet  photochemistry  methods  have 
also  been  investigated,  but  the  approach  has  not  proven  economically  viable,  it  is  a  low- 
yield  (high-cost)  process.  IC  stepper  technology  could  be  considered  for  via-drilling,  but 
is  not  viable  for  low-cost  MCM  substrates. 

A  laser  beam,  on  the  other  hand,  is  a  noncontacting  (i.e.,  zero-wear)  tool  that  can 
directly  drill  up  to  hundreds  of  vias  each  second.  The  unique  properties  of  lasers  make 
them  ideal  tools  for  micromachining  a  wide  range  of  materials.  Lasers  have  been  used  in 
the  semiconductor  industry  for  a  variety  of  tasks  including  photolithography  and  package 
welding.  Intense  laser  light  can  be  focused  to  a  small  spot  (a  few  microns  or  less),  so  that 
large  amounts  of  energy  can  be  applied  to  precise  locations,  allowing  materials  to  be 
drilled,  cut,  scribed,  annealed,  or  welded.  The  specific  interaction  between  a  laser  beam 
and  a  processed  material  depends  on  the  characteristics  of  the  material  and  the  wavelength 
of  the  laser.  Visible  and  infrared  laser  beams,  such  as  those  from  Nd:  YAG  (1.06  pm)  and 
CO2  lasers  (wavelength  of  10.6  pm),  can  provide  intense  local  heating.  Targeted  material 
is  essentially  boiled  off  or  vaporized.  Unfortunately,  the  intense  heating  often  causes 
thermal  damage  to  surrounding  material.  The  extent  of  the  damage  is  defined  by  the  "heat 
affected  zone"  (HAZ).  The  output  of  a  Nd:  YAG  laser  can  be  frequency  quadrupled  to  the 
deep  UV  (266  nm);  but  that  process  is  not  efficient,  and  the  lower-power  beam  cannot 
support  practical  throughput  rates. 

With  ultraviolet  lasers,  the  process  is  quite  different.  When  a  UV  photon  is 
absorbed  by  a  material  such  as  polyimide,  its  energy  is  not  converted  into  vibrations  of 
the  polyimide  molecules  (i.e.,  heat).  Instead,  it  directly  excites  the  electrons  forming  the 
molecular  bonds  that  hold  the  material  together.  Material  is  thus  atomized  in  a  process 
called  ablation.  This  instant  (cold)  process  produces  clean  edges  and  little  or  no  HAZ 
even  in  very  delicate  materials.  The  most  powerful  commercially  available  UV  lasers  are 
called  excimer  lasers.  An  excimer  laser  is  well-suited  to  precision  tasks,  such  as  drilling  via 
sites.  A  high-voltage  discharge  generates  pulsed  laser  light  from  a  mixture  of  inert  gas 
(helium  or  neon  mixed  with  argon,  krypton,  or  xenon)  and  a  halogen  gas  (fluorine  or 
hydrogen  chloride).  A  typical  industrial  excimer  laser  produces  pulse  energies  of  500  mJ, 
with  pulse  durations  of  only  10-40  nsec.  A  rapid  stream  of  such  high  energy  pulses  is 
ideal  for  percussive  drilling  applications.  The  laser  "pecks"  through  the  material  at  a 
constant  rate,  and  hole  depth  can  be  precisely  controlled  by  the  number  of  pulses 
delivered.  An  excimer  laser  allows  fast  drilling  because  it  delivers  repetition  rates  up  to 
500  Hz  with  overall  (time  averaged)  output  powers  as  high  as  200  W.  The  most  widely 
used  wavelengths  are  193,  248,  308,  and  351  nm.  The  energy  density  required  to  ablate 
polyimide  and  other  polymers  with  308-nm  light  is  300  mJ/cm^,  while  gold  and  copper 
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require  nearly  5  J/cm^  for  ablation.  Thus,  with  the  power  level  at  the  work  surface  set  to 
around  500  mJ/cm^,  holes  can  be  drilled  completely  through  the  polyimide  (down  to  the 
pads  and  ground  planes)  without  any  risk  of  damage  to  the  conductors. 

Masks  for  polymer  substrate  drilling  are  either  dielectric  or  chrome  patterned  on 
quartz.  They  absorb  or  reflect  the  laser  beam,  except  at  the  hole  locations.  Excimer  lasers 
have  a  finite  operating  cost  per  photon;  so  the  wasted  laser  light  translates  directly  into 
increased  processing  costs.  Also,  while  a  given  laser  has  a  maximum  output  power,  the 
light  fluence  reaching  the  surface  of  the  parts  largely  determines  the  maximum  throughput. 
Wasted  laser  light  raises  both  the  fixed  and  variable  production  costs.  Companies  are 
beginning  to  use  holographic  masks  and  micromachined  Fresnel-zone  plate  lenses  as  a 
viable  alternative  to  conventional  masks  for  substrate  machining.  Holographic  masks  are 
already  well-proven  for  micromachining  on  a  much  smaller  scale  in  microlithography  for 
integrated  circuits.  They  have  an  optical  efficiency  greater  than  80%,  since  they  redirect 
the  unwanted  light  to  the  desired  locations  instead  of  blocking  it.  Use  of  holographic 
masks  will  produce  much  higher  throughputs.  Processing  speed  depends  on  the  substrate 
thickness  and  the  laser  power.  With  50  pm  polyimide,  users  report  that  current  masking 
technology  allows  a  25  x  25-mm  area  to  be  processed  in  about  10  sec.  Cost  depends  on 
substrate  complexity  and  thickness,  but,  for  a  25  x  25  x  50  pm  polyimide,  the  cost  to  drill 
the  vias  ranges  from  less  than  a  cent  to  3-4  cents/substrate.  Holographic  masks  will 
decrease  those  costs  by  a  factor  of  about  five. 


3.6.  Task  Six:  Define  and  compare  alternative  interconnect  approaches  for  Mixed 
Mode  Modules. 


SA  has  examined  the  capabilities  of  normal  metal  electronic  interconnects  in  meeting 
requirements  for  2-D  and  3-D  MMM  applications.  Electrical  interconnects  greatly 
reduce  requirements  for  transducers  or  signal  conditioning  electronics  at  either  the  drive  or 
receiver  end  and  therefore  offer  minimal  excess  signal  propagation  delays  when  compared 
to  the  current  state  of  other  approaches.  Normal  (non-superconducting)  metal 
interconnects  are  easily  fabricated  in  high  densities  at  low  cost  using  chemical  processes  in 
MCM,  IC  or  PWB  processes.  The  principal  disadvantage  of  normal  metal  interconnects 
is  that  at  high  frequencies  the  signal  losses  increase.  Reducing  the  conductor  size  to 
achieve  the  high  interconnect  densities  required  for  MMMs  results  in  a  sharp  increase  in 
the  signal  losses  per  unit  length.  Severe  signal  loss  and  distortion  result  when  ohmic  line 
resistances  are  not  small  in  comparison  to  the  characteristic  impedance.  Copper 
conductor  sizes  for  MMM  type  line  lengths  (1"  to  1  foot)  are  in  the  10  pm  to  30  pm 
range. 

HTSC  materials  avoid  the  ohmic  resistance  problem  since  superconductors  have  no 
measurable  dc  resistance  and  stirface  resistance  values  that  are  negligibly  small  at  most 
frequencies  for  digital  systems.  Conductor  sizes  of  1.5  pm  to  3  pm  suffice  for  HTSC 
interconnects  which  makes  it  possible  to  achieve  extremely  high  interconnect  densities  in 
2-D  MMMs.  HTSC  have  frequency-independent  penetration  depth  and  do  not  show  the 
frequency  dependent  impedance  and  dispersion  effects  seen  in  normal  metal  lines.  HTSC 
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materials  also  offer  good  metal  shielding  and  ground  path  conductivity  which  results  in 
very  low  crosstalk  with  proper  conductor  geometries.  HTSC  interconnects  can  be  used 
for  long,  extremely  wide  bandwidth  interconnects  as  well  as  short  intra-MMM 
interconnects.  The  principal  disadvantage  is  the  need  for  cryogenic  cooling  to 
approximately  80°K  and  the  lack  of  mature  low-cost  fabrication  process  for 
manufacturing  high-density  multi-layer  HTSC  interconnects. 

Optical  interconnects  using  free-air  or  single-mode  fiber  optical  waveguides  offer 
extremely  low  dispersions.  Interconnect  lengths  are  on  the  order  of  tens  of  kilometers  for 
wide  signal  bandwidths.  With  low  dispersion  graded-index  multimode  fibers,  interconnect 
lengths  are  on  the  order  of  several  hundreds  of  meters,  much  greater  than  with  normal 
metal  transmission  lines.  The  disadvantage  with  this  technology  is  the  requirement  for 
transducers  at  both  ends.  Optical  transmitters  and  receivers  are  generally  expensive  and 
power  consuming  and  add  substantial  signal  propagation  delay  compared  to  their  normal 
metal  counterparts.  The  relatively  low  signal  amplitudes  from  the  photodetector  in  the 
receiver  must  be  amplified  to  digital  signal  levels  with  transistors  of  limited  gain- 
bandwidth  product  wliich  leads  to  significant  optical  receiver  propagation  delays.  For 
intra-module  interconnects  this  is  a  serious  technical  concern. 


3. 7.  Task  Seven:  Research  and  evaluate  thermal  management  concepts  that  achieve 
practical  heat  extraction  from  high  density  Mixed  Mode  Modules. 


A  key  problem  in  implementing  a  high  density  MMM  structure  is  developing 
methods  for  removing  the  heat  from  the  cube.  One  way  to  extract  heat  convectively  from 
an  intercoimected  MMM  package  is  to  use  an  electrically  insulating  super  high  thermal 
conductivity  substrate  material  (like  diamond)  to  conduct  the  heat  out  laterally  through 
the  substrate  material  to  heat  sinks  at  opposing  faces.  Artificial  diamond  substrates 
produced  by  state  of  the  art  chemical  vapor  deposition  (CVD)  techniques,  have  an 
extremely  high  thermal  conductivity  (700-1700  W/m°C)  and  exhibit  electrical  insulating 
characteristics.  Diamond  also  possesses  a  low  coefficient  of  thermal  expansion  (CTE) 
(1.5  ppm/°C)  which  is  close  to  Ae  CTE  value  of  Silicon  (2.6  ppm/°C).  Using  diamond  as 
substrates  to  conduct  the  heat  out  of  Ae  MMM  structure  frees  the  designer  to  fully 
utilize  all  of  Ae  space  between  modules  to  implement  a  high  area  density  array  of  vertical 
interconnects  between  all  adjacent  modules  in  the  stack  to  mmimize  interconnect  delays. 

One  approach  Aat  SA  has  evaluated  uses  1mm  thick  synthetic  diamond  substrates 
with  VLSI  chips  bonded  to  Aeir  upper  surface  and  a  high-density  multi-layer  Cu/BCB  or 
Cu/Polyimide  X-Y  interconnect  mat  on  the  bottom  side.  A  large  number  (on  Ae  order  of 
10,000)  of  4  mil  diameter  vias  are  laser-drilled  through  Ae  substrate  and  metallized  for 
passing  signals  from  Ae  chips  to  Ae  controlled  impedance  X-Y  interconnect  in  the 
bottom.  20  mil  diameter  fuzz  buttons  are  retained  m  Ae  spacer  boards  between  the 
diamond  substrates.  SA  has  determined  that  as  long  as  Ae  vertical,  Z  interconnects  and 
the  horizontal,  X-Y  MMM  intercoimects  have  Ae  same  characteristic  impedance,  mixed 
vertical-horizontal  runs  should  cause  no  signal  integrity  problem  for  implementing  a  3-D 
diamond  cube  MMM. 
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Thermal  conduction  models  for  both  (MMM-L)  and  additive  MMMs  (MMM-D), 
using  diamond  as  the  thermal  management  substrate,  have  been  completed.  These  models 
assess  the  feasibility  of  conducting  heat  laterally  through  the  diamond  substrates  to  heat 
sinks  located  at  the  board  edges.  This  frees  the  systems  designer  to  fully  use  all  of  the 
active  area  of  the  populated  substrate  to  implement  a  high  density  array  of  vertical 
interconnects.  These  models  have  illustrated  that  up  to  500  watts  can  be  dissipated  from 
a  typical  10  cm  x  10  cm  diamond  MMM  (while  maintaining  a  viable  chip  junction 
temperature)  by  using  forced  liquid  convection  technology  or  e\en  an  advanced  phase- 
change  heat  exchange  approach  such  as  spray  cooling.  Demonstration  MMM 
brassboards  interfaced  with  these  selected  board-edge  coolers  are  being  developed  to 
experimentally  verify  these  thermal  modeling  results. 

Coupled  with  their  silicon-matched  bonding  properties  and  electrical  insulator 
characteristics,  artificial  diamond  substrates  are  an  enabling  material  technology  for 
implementing  extremely  high  density  microelectronics  packaging.  Signal  latency  can  be 
reduced  by  a  factor  of  3  over  conventional  2-D  MMMs  through  the  use  of  an  area  array 
interconnected  3-D  configuration.  Applications  vary  from  supercomputers  through 
mainframes  and  high  performance  workstations.  The  technical  benefits  of  this  technology 
include:  (1)  minimization  of  interconnect  delays  in  supercomputers  operating  above  500 
MHz,  (2)  enhancing  inter-processor  I/O  bandwidths  in  MPPs  and  scaleable  processors, 
and  (3)  providing  for  resource  additions  to  workstations  such  as  vector  and  graphics 
processors.  From  a  military  standpoint,  this  technology  can  lead  to  new  operational 
capabilities  in  signal  and  image  reconnaissance  systems.  From  a  commercial  viewpoint, 
this  technology  may  accelerate  supercomputer  and  various  other  computer  product 
developments  approaching  TERAOPS  performance  capability. 


3.8.  Task  Eight:  Analyze  the  system  house  responsibility  for  driving  MMM 
technology  in  advanced  electronic  systems. 


There  is  currently  ongoing  research  by  the  Government,  Universities  and  industry  in 
the  area  of  Mixed  Mode  Modules.  The  significant  advantages  that  MMM  technology 
offers  in  size,  cost  and  performance  lend  it  to  being  integrated  into  high  performance 
systems.  Mixed  Mode  systems  are  becoming  increasingly  important  in  the  electronics 
industry  with  increased  reliance  on  analog  sensors  combined  with  ADCs  and  digital  signal 
processing  or  control.  Electronic  packaging  will  have  to  be  available  to  support  these 
systems.  Lowest  assembly  cost,  lowest  size/weight  and  highest  performance  all  dictate 
use  of  mixed  mode  modules.  Examples  of  military  and  commercial  systems  include 
wireless  communication,  global  positioning  systems,  tactical  information  assistants, 
electronic  warfare  systems,  automotive  collision  avoidance  systems  and  smart  sensors  of 
various  kinds. 

Mixed  mode  system  packages  are  being  produced  presently,  but  it  is  not  the 
product  of  well  developed  design  methodology.  Unfortunately  the  lack  of  a  cohesive 
infrastructure  is  hindering  the  integration  of  MMMs  into  high  volume,  viable  systems. 
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There  is  a  need  for  coordination  among  the  developers  and  the  users  of  this  technology  to 
ensure  that  the  passive  components,  optimized  materials,  low  cost  manufacturing 
techniques  and  more  general  CAD  tools  being  developed  will  in  fact  meet  the  particular 
needs  of  mixed  mode  systems.  One  of  the  major  issues  for  developing  MMMs  will  be  a 
determination  of  who  is  responsible  for  integrating  all  of  the  parts  onto  one  substrate. 
The  device  designers  are  only  familiar  with  their  particular  technology,  and  the  substrate 
manufacturers  are  concerned  with  producing  only  that  part  of  the  MMM.  The  system 
houses  will  drive  the  rest  of  the  MMM  infrastructure.  System  designers  who  are 
responsible  for  integrating  all  of  the  devices  on  an  MMM  substrate  will  need  to  become 
familiar  with  this  technology  before  the  pieces  are  in  place.  This  technology  is  in  the 
early  stages  of  development,  but  one  of  the  key  considerations  in  its  progress  is  to  avoid 
the  problem  of  integrating  too  many  devices  together  on  one  MMM  host  substrate.  The 
more  integration,  the  narrower  the  market  for  the  MMM  module  due  to  higher  NRE  costs 
and  a  longer  learning  curve. 


3. 9.  Task  Nine:  Evaluate  various  architectures  for  3D  MMM  structures. 

The  following  is  a  summary  of  the  benefits  that  3D  packaging  would  bring  to 
military  systems  and  outlines  how  additional  effort  in  this  area  could  make  low-cost  3D 
electronics  a  practical  reality  for  both  military  and  commercial  applications.  The 
principal  problem  that  low-cost  3-D  electronics  packaging  will  address  is  the  critical  need 
in  a  wide  range  of  military  systems  to  greatly  reduce  the  size,  weight,  and  volume  that 
various  vehicles  and  platforms  must  carry.  This  need  is  particularly  obvious  in  the  case 
of  aircraft,  manned  and  unmanned,  but  applies  equally  well  to  platforms  and  vehicles  all 
across  the  military.  A  reasonable  goal  for  3D  electronics  by  the  year  1999  would  be  the 
capability  to  shrink  electronics  for  various  military  applications  by  at  least  20x  in  volume 
and  20x  in  weight  in  comparison  to  conventional  packaging  such  as  VME  cards.  The  size 
reduction  will  also  substantially  reduce  the  average  signal  wire  length  and  therefore  cause  a 
reduction  in  the  digital  output  drive  power  requirements  fi-om  the  IC  chips. 

IC  devices  are  virtually  two-dimensional,  and  when  spread  in  a  single  plane  begin  to 
cover  substantial  areas  which  leads  to  longer  wire  length,  increased  latency  and  increased 
output  power  drive.  In  3-D  packaging  a  very  small  (~2mm)  vertical  stacking  pitch 
between  boards  is  achieved  by  mounting  the  bare  die  on  thin  MCM-type  substrates  and 
mounting  these  substrates  as  close  together  as  the  height  of  the  die  will  allow.  The  result 
is  a  “cube”  type  of  form  factor,  which  has  a  small  size  and  a  convenient  shape,  and 
because  of  the  low  overall  dimensions,  tends  to  offer  a  very  low  weight  overhead. 
Interconnections  between  the  boards  are  implemented  by  means  of  area  arrays  of  high 
density  vertical  interconnects.  Effective  use  of  these  vertical  interconnects  between 
boards  is  made  possible  by  high  densities  of  through  vias  from  the  top  to  the  bottom 
sides  of  the  various  MMM  boards  making  up  the  3D  stack.  The  results  of  this  3D 
electronics  packaging  approach  is  improvement  in  size,  weight,  volume,  power,  signal 
delay  and  bandwidths  which  translates  into  better  system  performance.  One  inherent 
problem  which  results  from  the  high  volume  density  of  chips  is  that  3-D  electronics 
packaging  represents  a  serious  challenge  in  thermal  management.  Though  some  promising 
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technologies  are  available  (e.g.,  diamond  substrates  as  a  means  of  distributing  the  heat), 
this  is  an  area  that  vvdll  require  further  investigation. 

In  the  3-D  MMM  packaging  concept  that  SA  evaluated,  double  rows  of  IC  chips 
are  flip-chip  mounted  to  the  MMM  substrate.  In  between  these  double  rows  are  patches 
of  high  density  area  interconnect  pads  which  contact  the  inter-board  Z-interconnect 
media.  While  one  approach  to  thermal  management  is  to  use  diamond  substrate  MMMs, 
the  arrangement  of  the  IC  die,  and  vertical  interconnects  into  rows  allows  for  alternative 
methods  of  thermal  management.  Through-the-bulk  spray  cooling  down  the  chip  rows 
(between  Z-interconnect  areas)  could  be  used  as  could  thin  cold  plates  sandwiched 
between  the  backs  of  chips  and  the  adjacent  boards  with  cutouts  for  the  Z-interconnects. 
It  is  important  to  note  that  in  general,  the  spray  cooling  requirements  for  the  3-D  MMMs 
will  be  substantially  lower  than  for  conventional  packaging  cooling.  Because  of  the 
shorter  interconnect  lengths  in  3-D  MMM  packaging,  the  output  drive  power 
requirements  of  the  IC  chips  is  reduced  by  >2x  over  conventional  2-D  packaging.  The 
very  small  fluid  volumes  in  spray  cooling  lead  to  very  low  pump  power  requirements. 
Liquid  spray  cooling  can  reduce  the  cooling  weight  and  size  for  a  3-D  MMM  by  16:1. 

Another  critical  area,  besides  thermal  management  is  the  development  of  high 
density,  low  cost  Z-interconnects.  What  is  required  is  a  combination  of  high  compliance, 
with  uniform,  reliable  low  electrical  contact  resistances.  The  compliance  can  be  achieved 
by  making  the  conductors  out  of  polymer  fibers,  with  only  a  thin  metal  coating  around 
the  body  of  the  fibers.  The  fibers  would  be  suspended  in  a  sort  of  random  pattern  in  an 
elastomer  matrix.  The  fiber  matrix  could  be  sliced  (perpendicular  to  the  fiber  direction)  in 
~1  mm  thick  sheets.  The  protruding  ends  of  the  polymer  fibers  would  be  metallized  and 
nickel  plated.  In  order  to  achieve  reliable  low  metal-to-metal  contact  resistance,  even 
when  the  mating  force  between  the  fiber  ends  and  the  mating  metal  pads  on  the  boards  is 
low,  particle  interconnects  are  used  on  the  fiber  ends.  Very  tiny  diamond  particles  are 
incorporated  in  the  plated  metal  on  the  fiber  ends  to  achieve  the  penetration  of  the  surface 
oxide  layers  on  the  mating  metal  pads.  Because  of  the  large  density  of  the  Z  interconnect 
fibers  within  the  matrix,  the  medium  would  not  have  to  be  aligned  to  the  mating  pads  on 
the  two  boards  being  mated,  and  therefore  there  is  no  need  to  customize  the  inter-layer 
interconnect  media  to  the  particular  pad  pattern  on  the  boards  being  mated.  This 
intercoimection  approach  is  an  inexpensive,  high  volume  commodity. 


3.10.  Task  Ten:  Evaluate  various  technologies  and  processes  that  can  be  used  to 
integrate  passive  devices  into  Mixed  Mode  Modules. 


In  current  practice,  implementing  complete  mixed-mode  systems  or  sub-systems 
requires  mounting  individual  packaged  devices,  MCMs,  shielded  or  unshielded  analog 
modules,  and  numerous  passive  components  for  filters  and  RF  processing  on  a  variety  of 
substrate  types.  This  approach  constrains  the  size,  weight  and  performance  capabilities 
of  defense  applications.  In  order  to  meet  military  requirements,  the  development  of 
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materials  and  processes  for  integrating  discrete  passives  (capacitors,  resistors,  and 
inductors)  directly  into  the  MMM  substrate,  is  important  since  embedded  passives 
technology  can  reduce  substrate  area  by  as  much  as  ten  times  while  providing  increased 
performance  and  reliability  at  1/4  of  the  cost.  In  addition,  many  secondary  issues 
associated  with  MMMs  such  as  cost  effectiveness,  yield,  testability,  and  reliability 
receive  little  attention  during  the  development  stage  of  a  technology,  though  their  impact 
is  as  important  as  first  level  research  questions. 

A  number  of  DARPA  programs  have  focused  on  technologies  for  the 
implementation  of  mixed-signal  MCMs.  MLxed-signal  MCM  technology  is  understood 
to  be  one  in  which  most  of  the  passive  components  can  be  fabricated  on  the  MCM 
substrate  itself.  The  advantages  include  reduced  cost  and  much  higher  densities  than 
conventional  surface  mount  technologies  can  provide.  The  current  focus  in  mixed-signal 
MCMs  is  for  mixed  RF/digital  circuitry.  The  first-order  focus  on  MCM  embedded 
resistor  and  capacitor  test  structures  is  on  evaluating  the  layer  properties  (sheet  resistance 
or  C/A),  since  it  is  easy  to  calculate  the  resistance  of  embedded  resistors  from  the  sheet 
resistance  and  the  capacitiance  of  embedded  capacitors  from  the  C/A  of  the  dielectric 
layer.  The  characteristics  of  inductors  is  deducible  from  known  materials  characteristics 
and  known  physics  (Maxwell’s  equations),  but  figuring  out  the  inductance,  let  alone  Q(f) 
for  an  embedded  inductor  is  not  at  all  simple.  The  quantitative  details  of  how  a  bent  piece 
of  wire  turns  into  an  inductor  is  still  pretty  much  a  mystery.  Appendix  A  provides  more 
detail  on  how  to  design/analyze  MCM  embedded  inductors,  including  meeting  MCM 
inductor  design  goals  efficiently  and  quickly. 
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4.  PERSONNEL  SUPPORTED 

Associate  -  4955.5  hours 
Scientist/Engineer  -  3487.5  hours 
Sr.  Scient/Engin  II  -  627  hours 

5.  PUBLICATIONS:  None 

6.  Interactions/Transitions 

6.1.  Participation  at  Meetings,  Conferences,  Seminars 

DARPA  Electronic  Packaging  and  Interconnect  Review  Meeting,  February  28- 
March  3,  1995,  at  the  Crystal  Gateway  Marriot  in  Crystal  City 

DARPA  Electronic  Packaging  and  Interconnect  Review  Meeting,  March  13-14, 

1 996,  at  the  Lansdowne  Convention  Center. 

6.2.  Consultative  and  advisory  functions: 

SA  has  provided  input  to  DARPA  Program  Managers  on  the  3D  Electronics 
information  that  is  provided  as  part  of  this  report.  This  information  was  used  by 
Dr.  James  Murphy,  Program  Manger,  DARPA/ETO,  in  a  briefing  to  the  Director 
of  DARPA,  Larry  Lynn  on  July  7,  1995. 

SA  provided  input  to  DARPA  Program  Managers  on  Embedded  Passives  for  a 
workshop  presentation  on  Embedded  Passives,  March  31,  1996 

6.3.  Transitions: 

As  a  result  of  research  performed  imder  this  contract  on  3D  interconnects,  SA  was 
able  to  provide  input  to  DARPA  that  will  lead  to  future  technology  applications 
using  3D  electronics. 

SA  has  also  provided  input  to  DARPA  in  a  new  program  using  MMMs  for  all- 
digital  receivers  and  developing  embedded  passive  technologies  for  the  MMMs. 

7.  Discoveries,  inventions,  Patent  disclosures:  none 

8.  HONORS/AWARDS:  None 


14 


APPENDIX  A 


Importance  of  Inductors  in  Mixed  Signal  MCMs 


15 


Importance  of  Inductors  in  Mixed-Signal  MCMs 

A  number  of  DARPA  programs  are  focused  on  technologies  for  the 
implementation  of  mixed-signal  MCMs.  The  MCM-L  consortium  (Dr.  Frank  Patten 
and  Dr.  Robert  Parker)  has  as  its  initial  target  product  an  rf  pocket  pager  product 
which  mixes  rf  circuitry  with  digital  signal  processing  chips.  Dr.  James  Murphy  has 
a  new  mixed-signal  foundry  program  focused  on  having  mixed-signal  MCM 
capabilities  commercially  available  in  the  high-density  MCM-D  technology  for  both 
military  and  commercial  applications.  Additionally,  Frank  Patten's  high-temperature 
superconductor  technologies  are  aimed  principally  at  high-density  rf/ microwave  or 
mixed-signal  applications. 

Mixed-signal  MCM  technology  is  imderstood  to  be  one  in  which  most  of  the 
myriad  of  passive  components  (resistors,  capacitors  and  inductors)  required  in 
typical  analog/rf  designs  can  be  fabricated  on  the  MCM  substrate  itself.  The 
advantages  of  doing  so  (as  pointed  out  by  Dr.  Nick  Neclario  in  several  recent  talks) 
are  reduced  cost  and  much  higher  densities  than  conventional  "chip  shooter"  surface 
moimt  technology  (SMT)  can  provide. 

In  addition  to  the  capability  for  the  fabrication  of  the  passive  elements 
themselves,  the  most  serious  problems  with  mixed-signal  MCMs  is  crosstalk 
between  the  high-level  digital  circuits  and  sensitive  analog/rf  circuitry.  The  most 
serious  problem  involves  coupling  through  power  supply  planes,  which,  for  best 
margin  and  highest  density,  necessitates  having  bypass  capacitor  planes  to  decouple 
the  power  planes  (and  correspondingly  reduce  the  ac  impedance  of  the  power 
planes).  This  is  already  available  in  MCM-D  technology;  the  standard  nCHIP  MCM 
processes  have  available  50  nF/cm^  bypass  capacitor  planes,  usable  both  for  power 
supply  bypass  functions  and  for  implementing  the  embedded  capacitors.  (Since 
nCHIP  is  the  prime  contractor  on  Jim  Murphy's  DARPA  mixed-signal  MCM 
foundry  program,  this  technical  issue  seems  imder  control  there.)  There  are  a 
number  of  efforts  under  way  to  bring  bypass  capacitor  planes  and  effective  (high 
capacitance  per  unit  area  [C/A])  embedded  capacitor  technologies  to  MCM-L. 
Within  Frank  Patten  and  Bob  Parker's  MCM-L  consortium,  there  is  an  embedded 
passives  effort  with  the  University  of  Arkansas  HiDEC  facility,  with  RPI  (Gene 
Rymaszewski)  developing  a  high  C/A  reactively-sputtered  capacitor  deposition 
technology  suitable  for  MCM-L.  Separately,  Jim  Murphy  has  a  3-M  program  for 
embedded  capacitors  in  MCM-Ls.  Hence,  while  a  high  capacitance  (high  C/ A) 
embedded  capacitor  technology  is  very  important  for  mixed-signal  MCMs,  this  issue 
is  being  treated  elsewhere,  an  will  not  be  discussed  further  here.  Similarly,  suitable 
embedded  resistor  technologies,  both  for  MCM-D  and  for  MCM-L,  are  being 
pursued  under  the  programs  mentioned  above,  and  elsewhere,  and  will  not  be 
discussed  further  here. 

As  noted  above,  the  focus  of  most  of  the  applications  interest,  both  military 
and  commercial,  in  mixed-signal  MCMs  is  for  mixed  rf/ digital  circuitry,  with  rf 
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frsquencies  up  to  about  1.9  GHz  of  principal  interest.  A  typical  characteristic  of  rf 
circuit  design  (as  compared  to  baseband  analog  designs)  is  that  the  operating  rf  or  if 
frequencies,  fs,  are  typically  much  larger  than  the  signal  bandwidths,  Afs*  Whereas 
in  baseband  analog  designs  the  signal  frequencies  are  generally  limited  by  the  R-C 
time  constant  at  the  signal  nodes,  this  is  not  the  case  in  rf  designs.  For  example,  if  the 
capacitance  at  a  circuit  node  is  Cjv  and  the  real  part  of  the  circmt  impedance  at  the 
node  is  Rn,  then  the  3dB  signal  frequency  will  be  of  the  order  of  (fs)-3dB=l/2JtRn  Cn- 
With  typical  device  and  interconnect  node  capacitances,  baseband  signal  frequencies 
are  generally  limited  to  some  tens  of  MHz  with  typical  device  and  interconnect 
technologies.  In  rf  designs,  the  operating  frequencies,  fs,  are  typically  far  above  this 
range,  and  it  is  the  signal  bandwidth.  Afs,  not  fs  itself,  which  is  subject  to  the  R-C 
time  constant  limit  cited  above.  This  rf  design  "miracle"  of  having  only  the  signal 
bandwidth,  not  the  signal  frequency  itself,  is  made  possible  by  the  use  of  inductors 
as  energy  storage  devices  in  the  circuits. 

The  key  to  this  rf  ’bandpass"  circuit  design  is  the  fact  that  while  capacitance  is 
ubiquitous,  both  in  the  active  devices  themselves  and  in  the  circuit  interconnects, 
and  the  capacitive  susceptance,  as  frequency  is  increased,  rapidly  dominates  the 
impedance  of  signal  nodes,  in  fact  the  capacitance  dissipates  no  energy  in  the  circuit 
and  generates  no  noise.  With  the  proper  use  of  inductors,  for  any  given  signal 
frequency,  fs,  it  is  possible  to  set  up  an  energy  exchange  between  lossless  inductors 
and  the  circuit  capacitance  such  that  the  capacitive  susceptance  effectively 
disappears,  with  respect  to  both  signal  amplitude  and  noise  considerations.  The 

reason  for  this  is  that  the  capacitive  susceptance  is  +j(]oC,  as  compared  to  -j/ coL  for 
the  inductive  susceptance  (or  correspondingly,  the  impedances  are  -j/ (oC  and  +jcoL 

respectively,  where  as  usual,  co=27if).  Because  these  have  opposite  signs  because  of 
the  different  way  in  which  capacitors  and  inductors  store  energy,  they  can  "cancel 
each  other  out".  Because  of  their  different  frequency  dependencies,  however,  they 
can  cancel  exactly  at  only  a  single  frequency,  fs,  or  "adequately"  within  a  range  of 
frequencies  within  the  bandwidth,  Afs,  around  this  center  frequency. 

Note  that  underlying  the  effectiveness  of  this  rf  bandpass  circuit  design 
concept  is  that  the  inductors  and  capacitors  are  energy  storage  devices,  as  opposed 
to  energy  dissipation  devices  (like  resistors).  If  both  inductors  and  capacitors  are 
lossless,  they  can  function  perfectly  in  these  bandpass  circuit  functions  (from  both 
the  frequency  response  and  circuit  noise  standpoints).  To  the  extent  that  they 
dissipate  a  significant  fraction  of  the  energy  they  are  supposed  to  be  just  storing,  the 
frequency  response  and  efficiency  of  the  rf  circuits  will  be  compromised  and  thermal 
noise  will  be  added  to  the  signals.  As  discussed  later  in  the  "Dependence  of 
Bandwidth,  Q  and  Noise  on  Rac  of  Inductor"  section  (Eqs.  33  -  46),  the  figure  of 
merit  for  inductors  and  capacitors  is  their  "quality  factor"  or  Q,  defined  (Eq.  33)  at  a 
frequency  f  as  Q(f)=7t  (maximmn  energy  stored)/ (energy  dissipated  per  cycle).  For 
nearly-ideal  rf  performance,  designers  want  Q  values  as  high  as  possible. 
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The  losses  in  both  the  energy  storage  materials  (e.g.,  dielectric  losses  in 
capacitors  or  hysteresis  losses  in  high  permeability  magnetic  core  materials)  and  in 
the  metal  conductors  associated  with  them  (IR  losses)  serve  to  limit  the  values  of  Q 
obtainable.  In  practice,  because  the  conductor  geometries  associated  with  capacitor 
structures  tend  to  be  highly  parallel  (or  very  low  resistance),  and  low  loss  tangent 
dielectric  materials  are  available  capable  of  g;iving  high  Q  values,  particularly  at 
modest  rf  frequencies.  The  situation  is  much  more  difficult  for  inductors.  Most  of  the 
inductors  used  in  rf  applications  are  "air  core",  meaning  that  the  active  magnetic 
"core"  region  in  which  most  of  the  magnetic  fields  reside  are  either  air  or  some  other 
dielectric  material  with  unity  relative  permeability  (|Xr=l)-  ^Vhile  these  materials 
have  no  hysteresis  losses,  rather  long  metal  thin  conductor  structures  are  required  to 
achieve  the  desired  inductance,  L,  values.  As  a  consequence,  the  ac  skin-effect 
resistance,  Rac/  values  for  the  inductors  tend  to  be  substantial,  and  hence  the  Q 

values  (Q=0)L/Rac  for  inductors)  tend  to  be  mediocre,  particularly  at  lower 
frequencies.  As  will  be  discussed  later,  this  problem  of  poor  inductor  Q  tends  to  be 
aggravated  as  the  size  of  the  inductor  structure  is  reduced.  Since  part  of  the  intent  of 
mixed-signal  MCM  technologies  is  to  be  able  to  be  able  to  increase  circuit 
densities/reduce  size  by  making  the  embedded  components  small  (keeping  the 
MCMs  small  is  also  a  part  of  the  cost  reduction  strategy),  making  MCM  embedded 
inductors  with  satisfactory  values  of  Q  (particularly  at  lower  rf  and  if  frequencies) 
within  size  constraints  is  a  challenging  undertaking.  However,  the  essential  role  that 
(reasonably)  high-Q  inductors  play  in  rf  circuit  design  makes  tackling  this  challenge 
very  important. 


Status  of  Embedded  Inductor  Design/Analysis  Capabilities 

I  have  attended  reviews  of  a  number  of  the  DARPA  mixed-signal  MCM 
programs  (e.g.,  Frank  Patten  and  Bob  Parker's  MCM-L  embedded  passives  reviews, 
and  others)  and  have  noticed  a  curious  characteristic  which  is  striking  similar 
between  them.  That  characteristic  is  the  highly  empirical  approach  to  the  whole 
issue  of  MCM  embedded  inductor  design  and  analysis,  in  sharp  contrast  to  that  used 
for  embedded  resistors  or  capacitors.  In  embedded  resistors,  everybody  knows  that 
the  terminal  resistance  will  be  given  as  the  product  of  the  number  of  squares  in  the 
resistor  layout  tomes  the  sheet  resistance  of  the  resistor  layer.  If  different  sized  or 
shaped  structures  are  fabricated,  it  is  understood  that  it  is  for  the  purpose  of 
evaluating  second-order  fabrication  factors  like  contact  resistance,  edge  effects  in 
resistor  pattern  etching,  etc.,  not  to  imderstand  how  a  resistor  works.  The  same  is 

true  for  capacitors.  Once  the  capacitance  per  unit  area,  C/ A=eo£r/ td/  bas  been 
evaluated,  everyone  knows  what  (apart  from  edge  fringing  effects)  what  the 
capacitance  of  a  given  structure  will  be.  When  different  sizes  or  shapes  are 
fabricated,  it  is  for  the  purpose  of  evaluating  second-order  fabrication  or  yield 
characteristics,  or  for  purposes  of  creating  structures  for  evaluation  of  operation  in 
different  regimes  of  frequency,  not  to  figure  out  what  the  capacitance  will  be. 
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While  the  first-order  focus  on  MCM  embedded  resistor  and  capacitor  test 
structures  is  on  evaluating  the  basic  layer  properties  (sheet  resistance  or  C/A),  that 
has  not  been  at  all  the  case  for  the  embedded  inductor  test  structures.  In  fact,  since 
the  basic  inductor  critical  materials  characteristics  are  well  known  (lossless  r=l  for 
the  core  and  reasonably  well  known  metal  resistivity  and  thickness  for  conductors), 
there  would  really  appear  to  be  no  need  for  such  first-order  inductor  test  structures 
at  all.  That  is,  the  characteristics  of  the  inductors  should  be  deducible  from  the 
known  materials  characteristics  and  the  known  physics  (Maxwell's  Equations).  In 
fact,  the  initial  inductor  test  structures  in  these  programs  were  done  with  mimmal 
a  priori  knowledge  of  what  to  expect  for  even  the  inductance,  let  alone  Q,  of  the 
fabricated  embedded  inductors. 

The  reason  for  this  is  quite  simple.  Whereas  it  is  extremely  simple  to  calculate 
the  resistance  of  embedded  resistors  from  the  sheet  resistance  of  a  resistor  layer,  and 
capacitance  of  capacitors  from  the  G/A  of  the  dielectric  layer,  figuring  out  the 
inductance,  let  alone  Q(f)  for  an  embedded  inductor  is  not  at  ^  simple.  While  most 
electrical  engineers  have  a  pretty  clear  quantitative  grasp  of  how  capacitors  and 
resistors  work,  the  quantitative  details  just  how  a  bent  piece  of  wire  turns  into 
inductance  remains  pretty  much  a  mystery  to  most  people.  The  analysis  work 
described  in  this  report  was  done,  and  the  report  written,  to  help  the  people  who 
need  to  know  just  how  to  design/analyze  MCM  embedded  inductors  to  avoid 
wasting  a  lot  of  time  building  test  patterns  (which,  in  the  end,  only  serve  to  verify 
Maxwell's  Equations)  zero  in  on  MCM  inductor  design  goals  efficiently  and  quickly. 
It  should  be  pointed  out  that,  because  of  the  substantial  complexity  of  the 
design/analysis  process  for  some  types  of  MCM  embedded  inductor  structures 
(particularly  planar  structures,  in  whidi  the  conductor  orientation  is  unfavorable,  or 
cases  when  other  conductors  [e.g.,  ground  planes  or  other  circuitry]  lie  in  proximity 
to  the  inductor  structure),  it  is  good  practice  to  verify  the  inductance  and  Q(f)  (plus 
relevant  crosstalk,  etc.)  performance  of  the  final  designs  with  realistic  test  patterns 
which  match  the  intended  application  as  close  as  possible. 


Implementation  of  Inductors  in  Mixed-Signal  MCMs 

The  approach  adopted  for  this  work,  in  order  to  focus  this  effort  on  the  real 
problems  of  mixed-signal  MCM  design,  was  to  conduct  frequent  meetings  with  the 
engineers  involved  in  actual  mixed-signal  MCM  work.  Because  the  DARPA  mixed- 
signal  MCM  foundry  effort  is  just  starting,  I  focused  most  of  my  effort  on  meeting 
with  engineers  at  nCHIP  to  evaluate  both  the  customer  requirements  (passed 
through  from  inputs  from  their  system  customers)  and  the  economic /process 
constraints  on  mixed-signal  MCM-Ds.  While  some  of  the  detailed  focus  is  on 
MCM-D  technology  issues,  in  fact,  the  analyses  are  equally  applicable  to  MCM-L  or 
other  technologies. 

One  of  the  key  questions  that  came  up  in  early  discussions  the  choice  between 
of  the  implementation  of  embedded  inductors  in  planar  spiral  form  versus  nCHIP’s 
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wirebond  implementation  of  solenoidal  inductors.  Some  aspects  of  this  question, 
such  as  the  much  greater  susceptibility  of  the  planar  spiral  inductors  to 
loss/ degradation  of  Q  due  to  substrate  conductance  if  non-highly  insulating  were 
considered  (and  won’t  be  discussed  here),  but  assuming  highly  insulating  substrates, 
issues  of  relative  area  requirements,  performance  capabilities  [e.g.,  Q  and  self¬ 
resonant  frequency],  crosstalk,  etc.  between  spirals  and  solenoids  are  of  importance. 
In  the  meetings,  it  was  pointed  out  that  since  all  planar  spirals  must  have  their  axes 
in  the  same  direction  (the  Z-direction,  assuming  the  MCM  substrate  is  the  X-Y  plane) 
the  crosstalk  between  them  will  be  maximum.  On  the  other  hand,  with  a  wirebond 
solenoid  embedded  inductor,  the  magnetic  dipole  axis  can  be  oriented  anywhere  in 
the  X-Y  plane,  so  it  is  possible  to  properly  orient  two,  even  fairly  closely  spaced, 
inductors  to  have  zero  crosstalk  (and  probably,  with  some  spacing  constraints,  to 
find  zero  crosstalk  orientations  for  three  inductors).  When  this  was  discussed  in  the 
meetings,  it  was  noted  that  for  this  to  be  useful,  we  have  to  have  the  means  to 
calculate  the  orientations  that  will  give  zero  crosstalk  values.  I  agreed  and  accepted 
the  task  of  developing  the  capability  to  calculate  the  magnetic  fields,  etc.  for  these 
structures  that  allow  such  questions  to  be  answered. 

Another  very  good  question  that  came  up  in  the  meetings  involves  the 
thickness  required  for  the  substrate  metallization  portion  of  the  loops  in  the 
wirebond  solenoids,  or  the  size  of  the  wirebond  wire  itself,  or  the  related  issue  of 
how  to  calculate  the  ac  conductor  (skin  effect)  losses  in  these  structures.  The  issue  is 
that,  in  the  idealized  view  of  a  solenoid  inductor,  the  magnetic  field  is  strong  inside 
the  solenoid,  and  essentially  zero  outside.  If  that  is  the  case,  then  the  high-frequency 
surface  current  on  the  conductors  (since  the  surface  current  density  must  equal  the 
tangential  H-field  at  the  surface)  would  be  near  zero  on  the  outside  of  the  conductor, 
and  hence  the  conductor  loss  would  be  essentially  twice  as  large  as  you  would 
normally  think,  since  only  half  of  the  periphery  contributes  to  conduction.  Also,  this 
would  mean  that  the  substrate  metal  might  as  well  be  made  only  half  as  many  skin 
depths  thick  as  would  ordinarily  be  specified,  since  the  bottom  side  doesn't 
contribute  to  high-frequency  conduction  anyway.  (As  will  be  shown  later,  this  turns 
out  not  to  be  the  case.) 

Another  design  issue  relating  to  the  wirebond  solenoid  MCM  embedded 
inductors  is  what  tum-to-tum  spacing  (or  pitch,  S,  of  the  loops  making  up  the 
solenoid)  is  optimum.  If  the  pitch  is  very  tight  (S  small),  the  magnetic  coupling 
between  loops  will  be  tighter,  improving  the  inductance,  L,  to  number  of  turns,  N, 
ratio,  and  hence  the  Q  should  be  improved.  However,  small  S  also  implies  greater 
turn-to-turn  capacitance,  and  hence  lower  self-resonant  frequency.  Also,  as  I 
suggested  in  the  kick-off  meeting  for  the  program,  the  use  of  ribbon  bonds  (instead 
of  normal  roimd  wire  wirebonds)  offers  an  ideal  way  to  sharply  reduce  resistance  of 
the  inductor  loops  without  increasing  eddy  current  losses  in  the  conductors  (because 
the  wide  direction  of  the  conductor  is  parallel  to  the  H-field,  instead  of 
perpendicular  to  it  as  in  a  planar  spiral).  However,  there  wull  be  some  sacrifice  in 
how  closely  the  turns  can  be  spaced  using  ribbon  bonds,  as  opposed  to  wirebonds. 
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3.5.  Task  Five:  Investigate  large  area  patterning  and  via  drilling  approaches  for 
Mixed  Mode  Modules. 

For  production  quantities  of  Multichip  Module  (MCM)  substrates,  mechanical 
drilling  of  vias  represents  neither  a  practical  nor  an  economic  option.  Prices  for  drill  bits 
smaller  than  250  pm  (1 0  mil)  in  diameter  range  up  to  $  1 5  each  -  with  a  realistic  maximum 
number  of  hits  between  500  and  1000.  Hit  rates  of  less  than  one  hole/sec  and  drill 
breakages  translate  into  very  slow  processing  and  low  yields.  Reactive  ion  etching  (RIE) 
presents  a  possible  alternative.  RIE  generates  heat,  however,  which  can  cause  dimensional 
stability  problems,  and  has  limited  spatial  resolution.  Wet  photochemistry  methods  have 
also  been  investigated,  but  the  approach  has  not  proven  economically  viable,  it  is  a  low- 
yield  (high-cost)  process.  IC  stepper  technology  could  be  considered  for  via-drilling,  but 
is  not  viable  for  low-cost  MCM  substrates. 

A  laser  beam,  on  the  other  hand,  is  a  noncontacting  (i.e.,  zero-wear)  tool  that  can 
directly  drill  up  to  hundreds  of  vias  each  second.  The  unique  properties  of  lasers  make 
them  ideal  tools  for  micromachining  a  wide  range  of  materials.  Lasers  have  been  used  in 
the  semiconductor  industry  for  a  variety  of  tasks  including  photolithography  and  package 
welding.  Intense  laser  light  can  be  focused  to  a  small  spot  (a  few  microns  or  less),  so  that 
large  amoxmts  of  energy  can  be  applied  to  precise  locations,  allowing  materials  to  be 
drilled,  cut,  scribed,  annealed,  or  welded.  The  specific  interaction  between  a  laser  beam 
and  a  processed  material  depends  on  the  characteristics  of  the  material  and  the  wavelength 
of  the  laser.  Visible  and  infrared  laser  beams,  such  as  those  from  Nd:  YAG  (1 .06  pm)  and 
CO2  lasers  (wavelength  of  10.6  pm),  can  provide  intense  local  heating.  Targeted  material 
is  essentially  boiled  off  or  vaporized.  Unfortunately,  the  intense  heating  often  causes 
thermal  damage  to  surrounding  material.  The  extent  of  the  damage  is  defined  by  the  "heat 
affected  zone"  (HAZ).  The  output  of  a  Nd:YAG  laser  can  be  frequency  quadrupled  to  the 
deep  UV  (266  nm);  but  that  process  is  not  efficient,  and  the  lower-power  beam  cannot 
support  practical  throughput  rates. 

With  ultraviolet  lasers,  the  process  is  quite  different.  When  a  UV  photon  is 
absorbed  by  a  material  such  as  polyimide,  its  energy  is  not  converted  into  vibrations  of 
the  polyimide  molecules  (i.e.,  heat).  Instead,  it  directly  excites  the  electrons  forming  the 
molecular  bonds  that  hold  the  material  together.  Material  is  thus  atomized  in  a  process 
called  ablation.  This  instant  (cold)  process  produces  clean  edges  and  little  or  no  HAZ 
even  in  very  delicate  materials.  The  most  powerful  commercially  available  UV  lasers  are 
called  excimer  lasers.  An  excimer  laser  is  well-suited  to  precision  tasks,  such  as  drilling  via 
sites.  A  high-voltage  discharge  generates  pulsed  laser  light  from  a  mixture  of  inert  gas 
(helium  or  neon  mbced  with  argon,  krypton,  or  xenon)  and  a  halogen  gas  (fluorine  or 
hydrogen  chloride).  A  typical  industrial  excimer  laser  produces  pulse  energies  of  500  mJ, 
with  pulse  durations  of  only  10-40  nsec.  A  rapid  stream  of  such  high  energy  pulses  is 
ideal  for  percussive  drilling  applications.  The  laser  "pecks"  through  the  material  at  a 
constant  rate,  and  hole  depth  can  be  precisely  controlled  by  the  number  of  pulses 
delivered.  An  excimer  laser  allows  fast  drilling  because  it  delivers  repetition  rates  up  to 
500  Hz  with  overall  (time  averaged)  output  powers  as  high  as  200  W.  The  most  widely 
used  wavelengths  are  193,  248,  308,  and  351  nm.  The  energy  density  required  to  ablate 
polyimide  and  other  polymers  with  308-nm  light  is  300  mJ/cm",  while  gold  and  copper 
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require  nearly  5  J/cm^  for  ablation.  Thus,  with  the  power  level  at  the  work  surface  set  to 
around  500  mJ/cm^,  holes  can  be  drilled  completely  through  the  polyimide  (down  to  the 
pads  and  ground  planes)  without  any  risk  of  damage  to  the  conductors. 

Masks  for  polymer  substrate  drilling  are  either  dielectric  or  chrome  patterned  on 
quartz.  They  absorb  or  reflect  the  laser  beam,  except  at  the  hole  locations.  Excimer  lasers 
have  a  finite  operating  cost  per  photon;  so  the  wasted  laser  light  translates  directly  into 
increased  processing  costs.  Also,  while  a  given  laser  has  a  maximum  output  power,  the 
light  fluence  reaching  the  surface  of  the  parts  largely  determines  the  maximum  throughput. 
Wasted  laser  light  raises  both  the  fixed  and  variable  production  costs.  Companies  are 
beginning  to  use  holographic  masks  and  micromachined  Fresnel-zone  plate  lenses  as  a 
viable  alternative  to  conventional  masks  for  substrate  machining.  Holographic  masks  are 
already  well-proven  for  micromachining  on  a  much  smaller  scale  in  microlithography  for 
integrated  circuits.  They  have  an  optical  efficiency  greater  than  80%,  since  they  redirect 
the  unwanted  light  to  the  desired  locations  instead  of  blocking  it.  Use  of  holographic 
masks  will  produce  much  higher  throughputs.  Processing  speed  depends  on  the  substrate 
thickness  and  the  laser  power.  With  50  pm  polyimide,  users  report  that  current  masking 
technology  allows  a  25  x  25-mm  area  to  be  processed  in  about  10  sec.  Cost  depends  on 
substrate  complexity  and  thickness,  but,  for  a  25  x  25  x  50  pm  polyimide,  the  cost  to  drill 
the  vias  ranges  from  less  than  a  cent  to  3-4  cents/substrate.  Holographic  masks  will 
decrease  those  costs  by  a  factor  of  about  five. 


3.6.  Task  Six:  Define  and  compare  alternative  interconnect  approaches  for  Mixed 
Mode  Modules. 


SA  has  examined  the  capabilities  of  normal  metal  electronic  interconnects  in  meeting 
requirements  for  2-D  and  3-D  MMM  applications.  Electrical  interconnects  greatly 
reduce  requirements  for  transducers  or  signal  conditioning  electronics  at  either  the  drive  or 
receiver  end  and  therefore  offer  minimal  excess  signal  propagation  delays  when  compared 
to  the  current  state  of  other  approaches.  Normal  (non-superconducting)  metal 
interconnects  are  easily  fabricated  in  high  densities  at  low  cost  using  chemical  processes  in 
MCM,  IC  or  PWB  processes.  The  principal  disadvantage  of  normal  metal  interconnects 
is  that  at  high  frequencies  the  signal  losses  increase.  Reducing  the  conductor  size  to 
achieve  the  high  interconnect  densities  required  for  MMMs  results  in  a  sharp  increase  in 
the  signal  losses  per  unit  length.  Severe  signal  loss  and  distortion  result  when  ohmic  line 
resistances  are  not  small  in  comparison  to  the  characteristic  impedance.  Copper 
conductor  sizes  for  MMM  type  line  lengths  (1"  to  1  foot)  are  in  the  10  pm  to  30  pm 
range. 

HTSC  materials  avoid  the  ohmic  resistance  problem  since  superconductors  have  no 
measurable  dc  resistance  and  surface  resistance  values  that  are  negligibly  small  at  most 
frequencies  for  digital  systems.  Conductor  sizes  of  1 .5  pm  to  3  pm  suffice  for  HTSC 
interconnects  which  makes  it  possible  to  achieve  extremely  high  interconnect  densities  in 
2-D  MMMs.  HTSC  have  frequency-independent  penetration  depth  and  do  not  show  the 
frequency  dependent  impedance  and  dispersion  effects  seen  in  normal  metal  lines.  HTSC 
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materials  also  offer  good  metal  shielding  and  ground  path  conductivity  which  results  in 
very  low  crosstalk  with  proper  conductor  geometries.  HTSC  interconnects  can  be  used 
for  long,  extremely  wide  bandwidth  interconnects  as  well  as  short  intra-MMM 
interconnects.  The  principal  disadvantage  is  the  need  for  cryogenic  cooling  to 
approximately  80°K  and  the  lack  of  mature  low-cost  fabrication  process  for 
manufacturing  high-density  multi-layer  HTSC  interconnects. 

Optical  interconnects  using  free-air  or  single-mode  fiber  optical  waveguides  offer 
extremely  low  dispersions.  Interconnect  lengths  are  on  the  order  of  tens  of  kilometers  for 
wide  signal  bandwidths.  With  low  dispersion  graded-index  multimode  fibers,  interconnect 
lengths  are  on  the  order  of  several  hundreds  of  meters,  much  greater  than  with  normal 
metal  transmission  lines.  The  disadvantage  with  this  technology  is  the  requirement  for 
transducers  at  both  ends.  Optical  transmitters  and  receivers  are  generally  expensive  and 
power  consuming  and  add  substantial  signal  propagation  delay  compared  to  their  normal 
metal  counterparts.  The  relatively  low  signal  amplitudes  from  the  photodetector  in  the 
receiver  must  be  amplified  to  digital  signal  levels  with  transistors  of  limited  gain- 
bandwidth  product  which  leads  to  significant  optical  receiver  propagation  delays.  For 
intra-module  interconnects  this  is  a  serious  technical  concern. 


S.  7.  Task  Seven:  Research  and  evaluate  thermal  management  concepts  that  achieve 
practical  heat  extraction  from  high  density  Mixed  Mode  Modules. 


A  key  problem  in  implementing  a  high  density  MMM  structure  is  developing 
methods  for  removing  the  heat  from  the  cube.  One  way  to  extract  heat  convectively  from 
an  interconnected  MMM  package  is  to  use  an  electrically  insulating  super  high  thermal 
conductivity  substrate  material  (like  diamond)  to  conduct  the  heat  out  laterally  through 
the  substrate  material  to  heat  sinks  at  opposing  faces.  Artificial  diamond  substrates 
produced  by  state  of  the  art  chemical  vapor  deposition  (CVD)  techniques,  have  an 
extremely  high  thermal  conductivity  (700-1700  W/m°C)  and  exhibit  electrical  insulating 
characteristics.  Diamond  also  possesses  a  low  coefficient  of  thermal  expansion  (CTE) 
(1.5  ppm/°C)  which  is  close  to  the  CTE  value  of  Silicon  (2.6  ppm/°C).  Using  diamond  as 
substrates  to  conduct  the  heat  out  of  the  MMM  structure  frees  the  designer  to  fully 
utilize  all  of  the  space  between  modules  to  implement  a  high  area  density  array  of  vertical 
interconnects  between  all  adjacent  modules  in  the  stack  to  minimize  interconnect  delays. 

One  approach  that  SA  has  evaluated  uses  1mm  thick  synthetic  diamond  substrates 
with  VLSI  chips  bonded  to  their  upper  surface  and  a  high-density  multi-layer  Cu/BCB  or 
Cu/Polyimide  X-Y  interconnect  mat  on  the  bottom  side.  A  large  number  (on  the  order  of 
10,000)  of  4  mil  diameter  vias  are  laser-drilled  through  the  substrate  and  metallized  for 
passing  signals  from  the  chips  to  the  controlled  impedance  X-Y  interconnect  in  the 
bottom.  20  mil  diameter  fuzz  buttons  are  retained  in  the  spacer  boards  between  the 
diamond  substrates.  SA  has  determined  that  as  long  as  the  vertical,  Z  interconnects  and 
the  horizontal,  X-Y  MMM  interconnects  have  the  same  characteristic  impedance,  mixed 
vertical-horizontal  runs  should  cause  no  signal  integrity  problem  for  implementing  a  3-D 
diamond  cube  MMM. 
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Thermal  conduction  models  for  both  (MMM-L)  and  additive  MMMs  (MMM-D), 
using  diamond  as  the  thermal  management  substrate,  have  been  completed.  These  models 
assess  the  feasibility  of  conducting  heat  laterally  through  the  diamond  substrates  to  heat 
sinks  located  at  the  board  edges.  This  frees  the  systems  designer  to  fully  use  all  of  the 
active  area  of  the  populated  substrate  to  implement  a  high  density  array  of  vertical 
interconnects.  These  models  have  illustrated  that  up  to  500  watts  can  be  dissipated  from 
a  typical  10  cm  x  10  cm  diamond  MMM  (while  maintaining  a  viable  chip  junction 
temperature)  by  using  forced  liquid  convection  technology  or  even  an  advanced  phase- 
change  heat  exchange  approach  such  as  spray  cooling.  Demonstration  MMM 
brassboards  interfaced  with  these  selected  board-edge  coolers  are  being  developed  to 
experimentally  verify  these  thermal  modeling  results. 

Coupled  with  their  silicon-matched  bonding  properties  and  electrical  insulator 
characteristics,  artificial  diamond  substrates  are  an  enabling  material  technology  for 
implementing  extremely  high  density  microelectronics  packaging.  Signal  latency  can  be 
reduced  by  a  factor  of  3  over  conventional  2-D  MMMs  through  the  use  of  an  area  array 
interconnected  3-D  configuration.  Applications  vary  from  supercomputers  through 
mainframes  and  high  performance  workstations.  The  technical  benefits  of  this  technology 
include:  (1)  minimization  of  interconnect  delays  in  supercomputers  operating  above  500 
MHz,  (2)  enhancing  inter-processor  I/O  bandwidths  in  MPPs  and  scaleable  processors, 
and  (3)  providing  for  resource  additions  to  workstations  such  as  vector  and  graphics 
processors.  From  a  military  standpoint,  this  technology  can  lead  to  new  operational 
capabilities  in  signal  and  image  reconnaissance  systems.  From  a  commercial  viewpoint, 
this  technology  may  accelerate  supercomputer  and  various  other  computer  product 
developments  approaching  TERAOPS  performance  capability. 


3.8.  Task  Eight:  Analyze  the  system  house  responsibility  for  driving  MMM 
technology  in  advanced  electronic  systems. 


There  is  currently  ongoing  research  by  the  Government,  Universities  and  industry  in 
the  area  of  Mixed  Mode  Modules.  The  significant  advantages  that  MMM  technology 
offers  in  size,  cost  and  performance  lend  it  to  being  integrated  into  high  performance 
systems.  Mixed  Mode  systems  are  becoming  increasingly  important  in  the  electronics 
industry  with  increased  reliance  on  analog  sensors  combined  with  ADCs  and  digital  signal 
processing  or  control.  Electronic  packaging  will  have  to  be  available  to  support  these 
systems.  Lowest  assembly  cost,  lowest  size/weight  and  highest  performance  all  dictate 
use  of  mbced  mode  modules.  Examples  of  military  and  commercial  systems  include 
wireless  communication,  global  positioning  systems,  tactical  information  assistants, 
electronic  warfare  systems,  automotive  collision  avoidance  systems  and  smart  sensors  of 
various  kinds. 

Mixed  mode  system  packages  are  being  produced  presently,  but  it  is  not  the 
product  of  well  developed  design  methodology.  Unfortunately  the  lack  of  a  cohesive 
infrastructure  is  hindering  the  integration  of  MMMs  into  high  volume,  viable  systems. 
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There  is  a  need  for  coordination  among  the  developers  and  the  users  of  this  technology  to 
ensure  that  the  passive  components,  optimized  materials,  low  cost  manufacturing 
techniques  and  more  general  CAD  tools  being  developed  will  in  fact  meet  the  particular 
needs  of  mixed  mode  systems.  One  of  the  major  issues  for  developing  MMMs  will  be  a 
determination  of  who  is  responsible  for  integrating  all  of  the  parts  onto  one  substrate. 
The  device  designers  are  only  familiar  with  their  particular  technology,  and  the  substrate 
manufacturers  are  concerned  with  producing  only  that  part  of  the  MMM.  The  system 
houses  will  drive  the  rest  of  the  MMM  infrastructure.  System  designers  who  are 
responsible  for  integrating  all  of  the  devices  on  an  MMM  substrate  will  need  to  become 
familiar  with  this  technology  before  the  pieces  are  in  place.  This  technology  is  in  the 
early  stages  of  development,  but  one  of  the  key  considerations  in  its  progress  is  to  avoid 
the  problem  of  integrating  too  many  devices  together  on  one  MMM  host  substrate.  The 
more  integration,  the  narrower  the  market  for  the  MMM  module  due  to  higher  NRE  costs 
and  a  longer  learning  curve. 


3. 9.  Task  Nine:  Evaluate  various  architectures  for  3D  MMM  structures. 

The  following  is  a  summary  of  the  benefits  that  3D  packaging  would  bring  to 
military  systems  and  outlines  how  additional  effort  in  this  area  could  make  low-cost  3D 
electronics  a  practical  reality  for  both  military  and  commercial  applications.  The 
principal  problem  that  low-cost  3-D  electronics  packaging  will  address  is  the  critical  need 
in  a  wide  range  of  military  systems  to  greatly  reduce  the  size,  weight,  and  volume  that 
various  vehicles  and  platforms  must  carry.  This  need  is  particularly  obvious  in  the  case 
of  aircraft,  manned  and  unmanned,  but  applies  equally  well  to  platforms  and  vehicles  all 
across  the  military.  A  reasonable  goal  for  3D  electronics  by  the  year  1999  would  be  the 
capability  to  shrink  electronics  for  various  military  applications  by  at  least  20x  in  volume 
and  20x  in  weight  in  comparison  to  conventional  packaging  such  as  VME  cards.  The  size 
reduction  will  also  substantially  reduce  the  average  signal  wire  length  and  therefore  cause  a 
reduction  in  the  digital  output  drive  power  requirements  from  the  IC  chips. 

IC  devices  are  virtually  two-dimensional,  and  when  spread  in  a  single  plane  begin  to 
cover  substantial  areas  which  leads  to  longer  wire  length,  increased  latency  and  increased 
output  power  drive.  In  3-D  packaging  a  very  small  (~2mm)  vertical  stacking  pitch 
between  boards  is  achieved  by  mounting  the  bare  die  on  thin  MCM-type  substrates  and 
mounting  these  substrates  as  close  together  as  the  height  of  the  die  will  allow.  The  result 
is  a  “cube”  type  of  form  factor,  which  has  a  small  size  and  a  convenient  shape,  and 
because  of  the  low  overall  dimensions,  tends  to  offer  a  very  low  weight  overhead. 
Interconnections  between  the  boards  are  implemented  by  means  of  area  arrays  of  high 
density  vertical  interconnects.  Effective  use  of  these  vertical  interconnects  between 
boards  is  made  possible  by  high  densities  of  through  vias  from  the  top  to  the  bottom 
sides  of  the  various  MMM  boards  making  up  the  3D  stack.  The  results  of  this  3D 
electronics  packaging  approach  is  improvement  in  size,  weight,  volume,  power,  signal 
delay  and  bandwidths  which  translates  into  better  system  performance.  One  inherent 
problem  which  results  from  the  high  volume  density  of  chips  is  that  3-D  electronics 
packaging  represents  a  serious  challenge  in  thermal  management.  Though  some  promising 
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technologies  are  available  (e.g.,  diamond  substrates  as  a  means  of  distributing  the  heat), 
this  is  an  area  that  will  require  further  investigation. 

In  the  3-D  MMM  packaging  concept  that  SA  evaluated,  double  rows  of  IC  chips 
are  flip-chip  mounted  to  the  MMM  substrate.  In  between  these  double  rows  are  patches 
of  high  density  area  interconnect  pads  which  contact  the  inter-board  Z-interconnect 
media.  While  one  approach  to  thermal  management  is  to  use  diamond  substrate  MMMs, 
the  arrangement  of  the  IC  die,  and  vertical  interconnects  into  rows  allows  for  alternative 
methods  of  thermal  management.  Through-the-bulk  spray  cooling  down  the  chip  rows 
(between  Z-interconnect  areas)  could  be  used  as  could  thin  cold  plates  sandwiched 
between  the  backs  of  chips  and  the  adjacent  boards  with  cutouts  for  the  Z-interconnects. 
It  is  important  to  note  that  in  general,  the  spray  cooling  requirements  for  the  3-D  MMMs 
will  be  substantially  lower  than  for  conventional  packaging  cooling.  Because  of  the 
shorter  interconnect  lengths  in  3-D  MMM  packaging,  the  output  drive  power 
requirements  of  the  IC  chips  is  reduced  by  >2x  over  conventional  2-D  packaging.  The 
very  small  fluid  volumes  in  spray  cooling  lead  to  very  low  pump  power  requirements. 
Liquid  spray  cooling  can  reduce  the  cooling  weight  and  size  for  a  3-D  MMM  by  16:1. 

Another  critical  area,  besides  thermal  management,  is  the  development  of  high 
density,  low  cost  Z-interconnects.  What  is  required  is  a  combmation  of  high  compliance, 
with  uniform,  reliable  low  electrical  contact  resistances.  The  compliance  can  be  achieved 
by  making  the  conductors  out  of  polymer  fibers,  with  only  a  thin  metal  coating  around 
the  body  of  the  fibers.  The  fibers  would  be  suspended  in  a  sort  of  random  pattern  in  an 
elastomer  matrix.  The  fiber  matrix  could  be  sliced  (perpendicular  to  the  fiber  direction)  in 
~1  mm  thick  sheets.  The  protruding  ends  of  the  polymer  fibers  would  be  metallized  and 
nickel  plated.  In  order  to  achieve  reliable  low  metal-to-metal  contact  resistance,  even 
when  the  mating  force  between  the  fiber  ends  and  the  mating  metal  pads  on  the  boards  is 
low,  particle  interconnects  are  used  on  the  fiber  ends.  Very  tiny  diamond  particles  are 
incorporated  in  the  plated  metal  on  the  fiber  ends  to  achieve  the  penetration  of  the  surface 
oxide  layers  on  the  mating  metal  pads.  Because  of  the  large  densit}"  of  the  Z  interconnect 
fibers  within  the  matrix,  the  medium  would  not  have  to  be  aligned  to  the  mating  pads  on 
the  two  boards  being  mated,  and  therefore  there  is  no  need  to  customize  the  inter-layer 
interconnect  media  to  the  particular  pad  pattern  on  the  boards  being  mated.  This 
interconnection  approach  is  an  inexpensive,  high  volume  commodity. 


3.10.  Task  Ten:  Evaluate  various  technologies  and  processes  that  can  be  used  to 
integrate  passive  devices  into  Mixed  Mode  Modules. 


In  current  practice,  implementing  complete  mixed-mode  systems  or  sub-systems 
requires  mounting  individual  packaged  devices,  MCMs,  shielded  or  unshielded  analog 
modules,  and  numerous  passive  components  for  filters  and  RF  processing  on  a  variety  of 
substrate  types.  This  approach  constrains  the  size,  weight  and  performance  capabilities 
of  defense  applieations.  In  order  to  meet  military  requirements,  the  development  of 
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materials  and  processes  for  integrating  discrete  passives  (capacitors,  resistors,  and 
inductors)  directly  into  the  MMM  substrate,  is  important  since  embedded  passives 
technology  can  reduce  substrate  area  by  as  much  as  ten  times  while  providing  increased 
performance  and  reliability  at  1/4  of  the  cost.  In  addition,  many  secondary  issues 
associated  with  MMMs  such  as  cost  effectiveness,  yield,  testability,  and  reliability 
receive  little  attention  during  the  development  stage  of  a  technology,  though  their  impact 
is  as  important  as  first  level  research  questions. 

A  number  of  DARPA  programs  have  focused  on  technologies  for  the 
implementation  of  mixed-signal  MCMs.  Mixed-signal  MCM  technology  is  understood 
to  be  one  in  which  most  of  the  passive  components  can  be  fabricated  on  the  MCM 
substrate  itself.  The  advantages  include  reduced  cost  and  much  higher  densities  than 
conventional  surface  mount  technologies  can  provide.  The  current  focus  in  mixed-signal 
MCMs  is  for  mixed  RF/digital  circuitry.  The  first-order  focus  on  MCM  embedded 
resistor  and  capacitor  test  structures  is  on  evaluating  the  layer  properties  (sheet  resistance 
or  C/A),  since  it  is  easy  to  calculate  the  resistance  of  embedded  resistors  from  the  sheet 
resistance  and  the  capacitiance  of  embedded  capacitors  from  the  C/A  of  the  dielectric 
layer.  The  characteristics  of  inductors  is  deducible  fi-om  known  materials  characteristics 
and  known  physics  (Maxwell’s  equations),  but  figuring  out  the  inductance,  let  alone  Q(f) 
for  an  embedded  inductor  is  not  at  all  simple.  The  quantitative  details  of  how  a  bent  piece 
of  wire  turns  into  an  inductor  is  still  pretty  much  a  mystery.  Appendix  A  provides  more 
detail  on  how  to  design/analyze  MCM  embedded  inductors,  including  meeting  MCM 
inductor  design  goals  efficiently  and  quickly. 
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4.  PERSONNEL  SUPPORTED 

Associate  -  4955.5  hours 
Scientist/Engineer  -  3487.5  hours 
Sr.  Scient/Engin  II  -  627  hours 

5.  PUBLICATIONS:  None 

6.  interactions/transitions 

6.1.  Participation  at  Meetings,  Conferences,  Seminars 

DARPA  Electronic  Packaging  and  Interconnect  Review  Meeting,  February  28- 
March  3,  1995,  at  the  Crystal  Gateway  Marriot  in  Crystal  City 

DARPA  Electronic  Packaging  and  Interconnect  Review  Meeting,  March  13-14, 
1996,  at  the  Lansdowne  Convention  Center. 

6.2.  Consultative  and  advisory  functions: 

SA  has  provided  input  to  DARPA  Program  Managers  on  the  3D  Electronics 
information  that  is  provided  as  part  of  this  report.  This  information  was  used  by 
Dr.  James  Murphy,  Program  Manger,  DARPA/ETO,  in  a  briefing  to  the  Director 
of  DARPA,  Larry  Lynn  on  July  7,  1995. 

SA  provided  input  to  DARPA  Program  Managers  on  Embedded  Passives  for  a 
workshop  presentation  on  Embedded  Passives,  March  31,  1996 

6.3.  Transitions: 

As  a  result  of  research  performed  under  this  contract  on  3D  interconnects,  SA  was 
able  to  provide  input  to  DARPA  that  will  lead  to  future  technology  applications 
using  3D  electronics. 

SA  has  also  provided  input  to  DARPA  in  a  new  program  using  MMMs  for  all- 
digital  receivers  and  developing  embedded  passive  technologies  for  the  MMMs. 

7.  DISCOVERIES,  INVENTIONS,  PATENT  DISCLOSURES:  NONE 

8.  HONORS/AWARDS:  NONE 
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Importance  of  Inductors  in  Mixed-Signal  MCMs 

A  number  of  DARPA  programs  are  focused  on  technologies  for  the 
implementation  of  mixed-signal  MCMs.  The  MCM-L  coi\sortium  (Dr.  Frank  Patten 
and  Dr.  Robert  Parker)  has  as  its  initial  target  product  an  rf  pocket  pager  product 
which  mixes  rf  circuitry  with  digital  signal  processmg  chips.  Dr.  James  Murphy  has 
a  new  mixed-signal  foundry  program  focused  on  having  mixed-signal  MCM 
capabilities  commercially  available  in  the  hig^-density  MCM-D  technology  for  both 
military  and  commercial  applications.  Additionally,  Frank  Patten's  high-temperature 
superconductor  technologies  are  aimed  principally  at  high-density  rf/ microwave  or 
mixed-signal  applications. 

Mbced-signal  MCM  technology  is  tmderstood  to  be  one  in  which  most  of  the 
myriad  of  passive  components  (resistors,  capacitors  and  inductors)  required  in 
typical  analog/rf  designs  can  be  fabricated  on  the  MCM  substrate  itself.  The 
advantages  of  doing  so  (as  pointed  out  by  Dr.  Nick  Neclario  in  several  recent  talks) 
are  reduced  cost  and  much  higher  densities  than  conventional  "chip  shooter"  surface 
mormt  technology  (SMT)  can  provide. 

In  addition  to  the  capability  for  the  fabrication  of  the  passive  elements 
themselves,  the  most  serious  problems  with  mixed-signal  MCMs  is  crosstalk 
between  the  high-level  digital  circuits  and  sensitive  analog/rf  circuitry.  The  most 
serious  problem  involves  coupling  through  power  supply  planes,  which,  for  best 
margin  and  highest  density,  necessitates  having  bypass  capacitor  planes  to  decouple 
the  power  planes  (and  correspondingly  reduce  the  ac  impedance  of  the  power 
planes).  This  is  already  available  in  MCM-D  technology;  the  standard  nCHIP  MCM 
processes  have  available  50  nF/cm^  bypass  capacitor  planes,  usable  both  for  power 
supply  bypass  functions  and  for  implementing  the  embedded  capacitors.  (Since 
nCHIP  is  the  prime  contractor  on  Jim  Murphy’s  DARPA  mixed-signal  MCM 
foundry  program,  this  technical  issue  seems  under  control  there.)  There  are  a 
number  of  efforts  imder  way  to  bring  bypass  capacitor  planes  and  effective  (high 
capacitance  per  unit  area  [C/A])  embedded  capacitor  technologies  to  MCM-L. 
Within  Frank  Patten  and  Bob  Parker's  MCM-L  consortium,  there  is  an  embedded 
passives  effort  with  the  University  of  Arkansas  HiDEC  facility,  with  RPI  (Gene 
Rymaszewski)  developing  a  high  C/A  reactively-sputtered  capacitor  deposition 
technology  suitable  for  MCM-L.  Separately,  Jim  Murphy  has  a  3-M  program  for 
embedded  capacitors  in  MCM-Ls.  Hence,  while  a  high  capacitance  (high  C/ A) 
embedded  capacitor  technology  is  very  important  for  mbced-signal  MCMs,  this  issue 
is  being  treated  elsewhere,  an  will  not  be  discussed  further  here.  Similarly,  suitable 
embedded  resistor  technologies,  both  for  MCM-D  and  for  MCM-L,  are  being 
pursued  under  the  programs  mentioned  above,  and  elsewhere,  and  will  not  be 
discussed  further  here. 

As  noted  above,  the  focus  of  most  of  the  applications  interest,  both  military 
and  commercial,  in  mixed-signal  MCMs  is  for  mixed  rf/ digital  circuitry,  with  rf 
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frequencies  up  to  about  1.9  GHz  of  principal  interest.  A  typical  characteristic  of  rf 
circuit  design  (as  compared  to  baseband  analog  designs)  is  that  the  operating  rf  or  if 
frequencies,  fs,  src  typically  much  larger  than  the  signal  bandwidths,  Afs.  Whereas 
in  baseband  analog  designs  the  signal  frequencies  are  generally  limited  by  the  R-C 
time  constant  at  the  signal  nodes,  this  is  not  the  case  in  rf  designs.  For  example,  if  the 
capacitance  at  a  circuit  node  is  Cn,  and  the  real  part  of  the  circuit  impedance  at  the 
node  is  Rn,  then  the  3dB  signal  frequency  will  be  of  the  order  of  (fs)-3dB=l/2jdRn  Cn- 
With  typical  device  and  interconnect  node  capacitances,  baseband  signal  frequencies 
are  generally  limited  to  some  tens  of  MHz  with  typical  device  and  interconnect 
technologies.  In  rf  designs,  the  operating  frequencies,  fs,  are  typically  far  above  this 
range,  and  it  is  the  signal  bandwidth.  Afs,  not  fs  itself,  which  is  subject  to  the  R-C 
time  constant  limit  cited  above.  This  rf  design  "miracle"  of  having  only  the  signal 
bandwidth,  not  the  signal  frequency  itself,  is  made  possible  by  the  use  of  inductors 
as  energy  storage  devices  in  the  circuits. 

The  key  to  this  rf  "bandpass"  circuit  design  is  the  fact  that  while  capacitance  is 
ubiquitous,  both  in  the  active  devices  themselves  and  in  the  circuit  interconnects, 
and  the  capacitive  susceptance,  as  frequency  is  increased,  rapidly  dominates  the 
impedance  of  signal  nodes,  in  fact  the  capacitance  dissipates  no  energy  in  the  circuit 
and  generates  no  noise.  With  the  proper  use  of  inductors,  for  any  given  signal 
frequency,  fs,  it  is  possible  to  set  up  an  energy  exchange  between  lossless  inductors 
and  the  circuit  capacitance  such  that  the  capacitive  susceptance  effectively 
disappears,  with  respect  to  both  signal  amplitude  and  noise  considerations.  The 

reason  for  this  is  that  the  capacitive  susceptance  is  +jcoC,  as  compared  to  -j/ coL  for 
the  inductive  susceptance  (or  correspondingly,  the  impedances  are  -]/ coC  and  +jcoL 

respectively,  where  as  usual,  co=27rf).  Because  these  have  opposite  signs  because  of 
the  different  way  in  which  capacitors  and  inductors  store  energy,  they  can  cancel 
each  other  out".  Because  of  their  different  frequency  dependencies,  however,  they 
can  cancel  exactly  at  only  a  single  frequency,  fs,  or  "adequately"  within  a  range  of 
frequencies  within  the  bandwidth,  Afs,  around  this  center  frequency. 

Note  that  underlying  the  effectiveness  of  this  rf  bandpass  circuit  design 
concept  is  that  the  inductors  and  capacitors  are  energy  storage  devices,  as  opposed 
to  energy  dissipation  devices  (like  resistors).  If  both  inductors  and  capacitors  are 
lossless,  they  can  function  perfectly  in  these  bandpass  circuit  functions  (from  both 
the  frequency  response  and  circuit  noise  standpoints).  To  the  extent  that  they 
dissipate  a  significant  fraction  of  the  energy  they  are  supposed  to  be  just  storing,  the 
frequency  response  and  efficiency  of  the  rf  circuits  will  be  compromised  and  thermal 
noise  will  be  added  to  the  signals.  As  discussed  later  in  the  "Dependence  of 
Bandwidth,  Q  and  Noise  on  Rac  of  Inductor"  section  (Eqs.  33  -  46),  the  figure  of 
merit  for  inductors  and  capacitors  is  their  "quality  factor"  or  Q,  defined  (Eq.  33)  at  a 
frequency  f  as  Q(f)=7t  (maximum  energy  stored)/(energy  dissipated  per  cycle).  For 
nearly-ideal  rf  performance,  designers  want  Q  values  as  high  as  possible. 
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The  losses  in  both  the  energy  storage  materials  (e.g.,  dielectric  losses  in 
capacitors  or  hysteresis  losses  in  high  permeability  magnetic  core  materials)  and  in 
the  metal  conductors  associated  with  them  (IR  losses)  serve  to  limit  the  values  of  Q 
obtainable.  In  practice,  because  the  conductor  geometries  associated  with  capacitor 
structures  tend  to  be  highly  parallel  (or  very  low  resistance),  and  low  loss  tangent 
dielectric  materials  are  available  capable  of  giving  high  Q  values,  particularly  at 
modest  rf  frequencies.  The  situation  is  much  more  difficult  for  inductors.  Most  of  the 
inductors  used  in  rf  applications  are  "air  core",  meaning  that  the  active  magnetic 
"core"  region  in  which  most  of  the  magnetic  fields  reside  are  either  air  or  some  other 
dielectric  material  with  unity  relative  permeability  (|Xr=l).  While  these  materials 
have  no  hysteresis  losses,  rather  long  metal  thin  conductor  structures  are  required  to 
achieve  the  desired  inductance,  L,  values.  As  a  consequence,  the  ac  skin-effect 
resistance,  Rac/  values  for  the  inductors  tend  to  be  substantial,  and  hence  the  Q 

values  (Q=a)L/Rac  for  inductors)  tend  to  be  mediocre,  particularly  at  lower 
frequencies.  As  will  be  discussed  later,  this  problem  of  poor  inductor  Q  tends  to  be 
aggravated  as  the  size  of  the  inductor  structure  is  reduced.  Since  part  of  the  intent  of 
mixed-signal  MCM  technologies  is  to  be  able  to  be  able  to  increase  circuit 
densities/reduce  size  by  making  the  embedded  components  small  (keeping  the 
MCMs  small  is  also  a  part  of  the  cost  reduction  strategy),  making  MCM  embedded 
inductors  with  satisfactory  values  of  Q  (particularly  at  lower  rf  and  if  frequencies) 
within  size  constraints  is  a  challenging  undertaking.  However,  the  essential  role  that 
(reasonably)  high-Q  inductors  play  in  rf  circuit  design  makes  tackling  this  challenge 
very  important. 


Status  of  Embedded  Inductor  Design/Analysis  Capabilities 

I  have  attended  reviews  of  a  number  of  the  DARPA  mixed-signal  MCM 
programs  (e.g.,  Frank  Patten  and  Bob  Parker’s  MCM-L  embedded  passives  reviews, 
and  others)  and  have  noticed  a  curious  characteristic  which  is  striking  similar 
between  them.  That  characteristic  is  the  highly  empirical  approach  to  the  whole 
issue  of  MCM  embedded  inductor  design  and  analysis,  in  sharp  contrast  to  that  used 
for  embedded  resistors  or  capacitors.  In  embedded  resistors,  everybody  knows  that 
the  terminal  resistance  will  be  given  as  the  product  of  the  number  of  squares  in  the 
resistor  layout  tomes  the  sheet  resistance  of  the  resistor  layer.  If  different  sized  or 
shaped  structures  are  fabricated,  it  is  imderstood  that  it  is  for  the  purpose  of 
evaluating  second-order  fabrication  factors  like  contact  resistance,  edge  effects  in 
resistor  pattern  etching,  etc.,  not  to  imderstand  how  a  resistor  works.  The  same  is 
true  for  capacitors.  Once  the  capacitance  per  unit  area,  C/A=£oer/td/  bas  been 
evaluated,  everyone  knows  what  (apart  from  edge  fringing  effects)  what  the 
capacitance  of  a  given  structure  will  be.  When  different  sizes  or  shapes  are 
fabricated,  it  is  for  the  purpose  of  evaluating  second-order  fabrication  or  yield 
characteristics,  or  for  purposes  of  creating  structures  for  evaluation  of  operation  in 
different  regimes  of  frequency,  not  to  figure  out  what  the  capacitance  will  be. 
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While  the  first-order  focus  on  MCM  embedded  resistor  and  capacitor  test 
structures  is  on  evaluating  the  basic  layer  properties  (sheet  resistance  or  C/A),  that 
has  not  been  at  all  the  case  for  the  embedded  inductor  test  structures.  In  fact,  since 
the  basic  inductor  critical  materials  characteristics  are  weU  known  (lossless  p,  r=l  for 
the  core  and  reasonably  well  known  metal  resistivity  and  thickness  for  conductors), 
there  would  really  appear  to  be  no  need  for  such  first-order  inductor  test  structures 
at  all.  That  is,  the  characteristics  of  the  inductors  should  be  deducible  from  the 
known  materials  characteristics  and  the  known  physics  (Maxwell's  Equations).  In 
fact,  the  initial  inductor  test  structures  in  these  programs  were  done  with  minimal 
a  priori  knowledge  of  what  to  expect  for  even  the  inductance,  let  alone  Q,  of  the 
fabricated  embedded  inductors. 

The  reason  for  this  is  quite  simple.  Whereas  it  is  extremely  simple  to  calculate 
the  resistance  of  embedded  resistors  from  the  sheet  resistance  of  a  resistor  layer,  and 
capacitance  of  capacitors  from  the  C/A  of  the  dielectric  layer,  figuring  out  the 
inductance,  let  alone  Q(f)  for  an  embedded  inductor  is  not  at  all  simple.  While  most 
electrical  engineers  have  a  pretty  clear  quantitative  grasp  of  how  capacitors  and 
resistors  work,  the  quantitative  details  just  how  a  bent  piece  of  wire  turns  into 
inductance  remains  pretty  much  a  mystery  to  most  people.  The  analysis  work 
described  in  this  report  was  done,  and  the  report  written,  to  help  the  people  who 
need  to  know  just  how  to  design/analyze  MCM  embedded  inductors  to  avoid 
wasting  a  lot  of  time  building  test  patterns  (which,  in  the  end,  only  serve  to  verify 
Maxwell’s  Equations)  zero  in  on  MCM  inductor  design  goals  efficiently  and  qmckly. 
It  should  be  pointed  out  that,  because  of  the  substantial  complexity  of  the 
design /analysis  process  for  some  types  of  MCM  embedded  inductor  structures 
(particularly  planar  structures,  in  which  the  conductor  orientation  is  unfavorable,  or 
cases  when  other  conductors  [e.g.,  groimd  planes  or  other  circuitry]  lie  in  proximity 
to  the  inductor  structure),  it  is  good  practice  to  verify  the  inductance  and  Q(f)  (plus 
relevant  crosstalk,  etc.)  performance  of  the  final  designs  with  realistic  test  patterns 
which  match  the  intended  application  as  close  as  possible. 


Implementation  of  Inductors  in  Mixed-Signal  MCMs 

The  approach  adopted  for  this  work,  in  order  to  focus  this  effort  on  the  real 
problems  of  mixed-signal  MCM  design,  was  to  conduct  frequent  meetings  with  the 
engineers  involved  in  actual  mixed-signal  MCM  work.  Because  the  DARPA  mixed- 
signal  MCM  foundry  effort  is  just  starting,  I  focused  most  of  my  effort  on  meeting 
with  engineers  at  nCHIP  to  evaluate  both  the  customer  requirements  (passed 
through  from  inputs  from  their  system  customers)  and  the  economic /process 
constraints  on  mixed-signal  MCM-Ds.  While  some  of  the  detailed  focus  is  on 
MCM-D  technology  issues,  in  fact,  the  analyses  are  equally  applicable  to  MCM-L  or 
other  technologies. 

One  of  the  key  questioirs  that  came  up  in  early  discussions  the  choice  between 
of  the  implementation  of  embedded  inductors  in  planar  spiral  form  versus  nCHIP’s 
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wirGbond  implGniGntation  of  solanoidal  inductors.  SomG  aspects  of  this  question, 
such  as  the  much  greater  susceptibility  of  the  planar  spiral  inductors  to 
loss / degradation  of  Q  due  to  substrate  conductance  if  non-highly  insulating  were 
considered  (and  won't  be  discussed  here),  but  assuming  highly  insulating  substrates, 
issues  of  relative  area  requirements,  performance  capabilities  [e.g.,  Q  and  self- 
resonant  frequency],  crosstalk,  etc.  between  spirals  and  solenoids  are  of  importance. 
In  the  meetings,  it  was  pointed  out  that  since  all  planar  spirals  must  have  their  axes 
in  the  same  direction  (the  Z-direction,  assuming  the  MCM  substrate  is  the  X-Y  plane) 
the  crosstalk  between  them  will  be  maximum.  On  the  other  hand,  with  a  wireborid 
solenoid  embedded  inductor,  the  magnetic  dipole  axis  can  be  oriented  anywhere  in 
the  X-Y  plane,  so  it  is  possible  to  properly  orient  two,  even  fairly  closely  spaced, 
inductors  to  have  zero  crosstalk  (and  probably,  with  some  spacing  constramts,  to 
find  zero  crosstalk  orientations  for  three  inductors).  When  this  was  discussed  in  the 
meetings,  it  was  noted  that  for  this  to  be  useful,  we  have  to  have  the  means  to 
calculate  the  orientations  that  will  give  zero  crosstalk  values.  I  agreed  and  accepted 
the  task  of  developing  the  capability  to  calculate  the  magnetic  fields,  etc.  for  these 
structures  that  allow  such  questions  to  be  answered. 

Another  very  good  question  that  came  up  in  the  meetings  involves  the 
thickness  required  for  the  substrate  metallization  portion  of  the  loops  in  the 
wirebond  solenoids,  or  the  size  of  the  wirebond  wire  itself,  or  the  related  ipue  of 
how  to  calculate  the  ac  conductor  (skin  effect)  losses  in  these  structures.  The  issue  is 
that,  in  the  idealized  view  of  a  solenoid  inductor,  the  magnetic  field  is  strong  inside 
the  solenoid,  and  essentially  zero  outside.  If  that  is  the  case,  then  the  high-frequency 
surface  current  on  the  conductors  (since  the  surface  current  density  must  equal  the 
tangential  H-field  at  the  surface)  would  be  near  zero  on  the  outside  of  the  conductor, 
and  hence  the  conductor  loss  would  be  essentially  twice  as  large  as  you  would 
normally  think,  since  only  half  of  the  periphery  contributes  to  conduction.  Also,  tlus 
would  mean  that  the  substrate  metal  might  as  well  be  made  only  half  as  many  skin 
depths  thick  as  would  ordinarily  be  specified,  since  the  bottom  side  doesn't 
contribute  to  high-frequency  conduction  anyway.  (As  will  be  shown  later,  this  turns 
out  not  to  be  the  case.) 

Another  design  issue  relating  to  the  wirebond  solenoid  MCM  embedded 
inductors  is  what  tum-to-tum  spacing  (or  pitch,  S,  of  the  loops  making  up  the 
solenoid)  is  optimum.  If  the  pitch  is  very  tight  (S  small),  the  magnetic  coupling 
between  loops  will  be  tighter,  improving  the  inductance,  L,  to  number  of  turns,  N, 
ratio,  and  hence  the  Q  should  be  improved.  However,  small  S  also  implies  greater 
turn-to-turn  capacitance,  and  hence  lower  self-resonant  frequency.  Also,  as  I 
suggested  in  the  kick-off  meeting  for  the  program,  the  use  of  ribbon  bonds  (instead 
of  normal  round  wire  wirebonds)  offers  an  ideal  way  to  sharply  reduce  resistance  of 
the  inductor  loops  without  increasing  eddy  current  losses  in  the  conductors  (because 
the  wide  direction  of  the  conductor  is  parallel  to  the  H-field,  instead  of 
perpendicular  to  it  as  in  a  planar  spiral).  However,  there  will  be  some  sacrifice  in 
how  closely  the  turns  can  be  spaced  using  ribbon  bonds,  as  opposed  to  wirebonds. 
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and  hence  a  detailed  quantitative  analysis  of  the  effects  on  L  and  R  that  determine  Q 
will  be  necessary  to  optimize  Q  for  MCM  embedded  solenoid  inductors. 

In  light  of  these  numerous,  and  important,  design  questions  about  embedded 
inductors  in  mixed-signal  MCMs  requiring  detailed  quantitative  answers,  I  took 
upon  myself  the  task  of  coming  up  with  analytical  and  numerical  tools  of  analysis 
capable  of  dealing  with  both  planar  spiral  and  solenoidal  embedded  inductor 
structures.  This  has  turned  out  to  be  quite  a  bit  of  work,  but  very  valuable  in 
providing  exactly  the  kinds  of  insights  and  answers  needed  for  the  design  and 
analysis  of  MCM  embedded  inductor  structures. 


Stacked  or  Concentric  Loop  Approaches  to  Modeling  the  Magnetic 
Fields  and  Inductances  of  Solenoidal  or  Spiral  Embedded  Inductors 

There  exist  in  handbooks  innumerable  empirical  formulas  for  calculating  the 
approximate  inductance  of  various  spiral  or  solenoidal  inductor  structures. 
Unforhmately,  in  applying  such  formulas  the  designer  has  limited  knowledge  of 
whether  the  particular  formula  is  accurate  in  the  range  of  design  parameters  in 
which  he  is  interested,  and  they  provide  no  insight  at  all  into  the  type  of  design 
questions  that  were  cited  above.  Further,  they  only  serve  to  deepen  the  sense  of 
abstruseness  which  seems  to  surroimd  magnetic  fields  and  inductor  structures. 
While  no  designer  has  any  problem  converting  bulk  or  sheet  resistivity  and 
geometry  into  resistance,  or  converting  dielectric  constant  and  at  least  simple 
geometries  into  capacitance,  just  how  a  bent  wire  and  the  permeability  of  free  space 
gets  turned  into  magnetic  fields  and  inductance  is  pretty  much  a  mystery  to  most 
designers.  Since  the  principal  purpose  here  is  to  achieve  insight,  as  well  as 
quantitative  treatment  of  some  fairly  subtle  questions  about  MCM  embedded 
inductors  and  their  implementation,  I  decided  at  the  outset  that  the  only  approach 
which  could  accomplish  these  goals  would  be  to  start  from  first  principles  along  an 
analysis  path  that  is  both  imderstandable  and  mathematically  tractable.  I  found  the 
solution  to  both  of  these  requirements  in  the  recognition  that  both  the  planar  spiral 
inductors  and  the  solenoid  inductor  structures  may  both  be  viewed  simply  as  an 
assemblage  of  coaxial  circular  loops. 

Figure  1  shows  a  planar  spiral  inductor  structure  as  might  be  used  as  an 
MCM  embedded  inductor.  The  structure  is  drawn  in  Figure  1,  not  as  a  spiral,  but  as 
collection  of  nested  coaxial,  coplanar  circular  loops  of  decreasing  radii,  R.  It  is 
intuitively  obvious  that  this  structure  will  give  virtually  identical  inductance  and  H- 
field  characteristics  to  a  true  spiral,  and  would  probably  be  closer  to  the  kind  of 
layout  that  would  be  practically  done  in  MCMs  (actually,  octagons  are  typically 
used  to  approximate  the  circles).  Recognizing  this,  the  magnetic  field  produced  by 
an  N-tum  (N-loop)  spiral  can  be  obtained  as  the  sum  of  the  fields  from  the  current,  I, 
passing  through  the  collection  of  loops,  and  its  inductance  may  be  obtained  as  the 
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sum  of  the  self  inductances  of  the  N  loops  plus  all  of  the  combination  of  mutual 
inductances  between  the  different  loops  (N^  terms  in  all). 

Figure  2  shows  the  same  concept  applied  to  the  solenoid  inductor  case.  Here, 
instead  of  treating  the  actual  helical  conductor  geometry,  it  is  intuitively  obvious 
that  (since  the  conductor  spacings  and  radii  are  the  same)  the  treatment  as  a  series 
of  stacked  circular  coaxial  loops,  all  of  identical  radius,  R,  separated  from  one 
another  by  a  stacking  pitch,  S,  will  give  essentially  the  same  fields  and  inductance  as 
a  helix  of  radius,  R,  and  winding  pitch,  S.  As  for  the  spiral,  the  magnetic  field 
produced  by  an  N-tum  (N-loop)  solenoid  at  any  point  can  be  obtained  by  summing 
the  magnetic  fields  produced  by  the  current,  I,  passing  through  each  of  the  loops  in 
the  stack  (recognizing  that  each  will  be  at  a  different  height  relative  to  the 
measurement  point),  and  the  inductance  of  the  solenoid  may  be  obtained  as  the  sum 
of  the  self  inductances,  Lii,  of  the  N  loops  (which  will  be  identical  for  all  N  loops) 
plus  all  of  the  combination  of  mutual  inductances,  Mij,  between  the  different  loops 
(which  will  depend  on  the  vertical  separation  between  the  pairs  of  loops  in  the 
stack).  Again,  there  will  be  N  Lii  terms  and  N(N-1)  Mij  terms  (for  i^^j)  or  total 
terms  to  be  summed  to  give  the  inductance  for  the  N-loop  solenoid  (because 
Mij=Mji,  only  half  of  these  are  independent,  and  all  of  the  Lii  the  same  for  the 
solenoid  [but  not  for  the  spiral]).  While  summing  these  magnetic  field  or  mutual 
inductance  terms  does  represent  an  increasing  amoimt  of  arithmetic  as  the  number 
of  turns  or  loops,  N,  increases,  it  is  inherently  very  simple,  and  the  type  of 
computation  which,  with  proper  organization,  can  be  nicely  handled  in  a 
spreadsheet  (EXCEL  4.0  was  used  in  this  work). 

A  very  nice  feature  of  this  approach  of  modeling  MCM  embedded  spiral  or 
solenoidal  inductor  structures  as  a  collection  of  coaxial  loops  is  that  we  have  only 
one  electromagnetic  fields  problem  to  solve  (actually,  one  magnetic  fields  problem, 
since  we  are  analyzing  the  low-frequency  magnetic  fields  and  inductances  where  the 
dimensions  of  the  inductor  structure  are  small  in  comparison  to  a  wavelength). 
Ordinarily,  the  bane  of  E&M  is  that  every  time  you  want  to  examine  a  new  problem, 
you  are  back  to  solving  Maxwell's  Equations  for  the  new  geometry.  We  have  only 
one  general  magnetic  field  case  to  solve;  calculating  the  H-field  above  a  circular  loop 
of  radius  R,  at  a  measurement  point  at  some  height,  d,  above  the  plane  of  the  loop 
and  radius,  r  from  its  axis.  We  have  only  one  general  mutual  inductance  case  to 
solve;  the  case  of  the  mutual  inductance  between  two  coaxial  circular  loops  of  two 
different  radii  displaced  by  a  vertical  distance,  d.  Of  course,  for  the  spiral  inductor, 
all  of  the  loops  are  coplanar,  but  this  is  just  the  general  case  with  d=0,  so  we  still 
have  only  one  E&M  geometry  to  solve  (and  fortimately,  it  is  not  very  difficult). 
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Calculation  of  the  Low  Frequency  Magnetic  Field  Above  a  Loop 

I  had  considered  putting  the  details  of  the  analysis  for  the  calculation  of  the 
magnetic  fields  for  the  loop  an  appendix,  but  in  light  of  the  purpose  of  all  of  this  to 
dispel  the  mystery  surrounding  inductors,  that  would  be  counterproductive  to  our 
goal.  Further,  the  calculation  is  actually  very  simple,  with  very  few  "gory  details", 
and  the  simplicity  is  a  key  part  of  the  point,  so  let's  just  quickly  walk  through  the 
process  from  first  principles  to  getting  plots  of  the  magnetic  field  strength  vs.  radial 
position  at  various  heights  above  a  circular  loop. 

In  the  first  place,  because  we  are  interested  in  the  low-frequency,  "static" 
fields  (where  the  dimensions  of  the  inductor  structure  are  small  in  comparison  to  a 
wavelength),  we  don't  have  to  cope  with  the  sunultaneous  solution  of  NIaxwell  s 
Equations,  but  rather  can  go  directly  to  Ampere's  law  for  the  loop  geometry  defined 
in  Figure  3.  We  will  consider  only  situations  in  which  the  permeability,  \i,  of  the 
medium  is  constant  and  isotropic  (and  for  the  actual  calculations,  equal  to  the 
permeability  of  free  space,  {io=47t  x  10"^  henries/ meter),  so  at  any  point  r  the  B-field 
is  simply  given  as 


B(r)  =  pH(r)  Eq.l 

where  bold  face  type  is  used  to  indicate  vector  quantities.  Because  the  numbers  are 
more  convenient  for  H-fields,  the  field  plots  will  be  given  in  terms  of  H-field,  only 
converting  to  B-field  when  necessary  to  calculate  inductance,  from  the  B-field. 

Amperes  law  is  very  simple.  If  a  differential  segment  of  conductor,  dl,  of 
length,  dl,  located  at  some  point  in  space,  r',  carries  a  current,  I,  then  the  magnitude 
of  the  differential  H-field  produced  at  a  point  r  by  this  current  will  be  given  by 

dH(r)  =  [I  dl  sinb]  /[4j:  I  r  -  r’  1 2]  Eq.  2 

where  I  r-r'  1 2  is  just  the  square  of  the  distance  between  r  and  r’,  and  sincj)  is  just  the 
sine  of  the  angle  between  the  vector  (r-r')  and  dl.  The  direction  of  H(r)  is  mutually 
perpendicular  to  both  (r  -  r')  and  dl.  While  this  is  all  you  need  to  solve  the  loop  field 
problem,  it  turns  out  to  be  a  lot  simpler  to  keep  the  trigonometry  straight  if  you  use 
the  cross-product  vector  notation,  in  which  form  Eq.  2  is  expressed  as 

dH(r)  =  [I  dl  X  (r  -  r')]/[4K  I  r  -  r'  1 3]  Eq.  3 

The  advantage  of  this  notation  is  that  it  keeps  all  of  the  work  very  simple  and  easy. 
The  disadvantage  is  that  you  have  to  remember  the  definition  of  the  cross  product, 
which  is  simply 


A  X  B  =  ax(AyBz  -  AzBy)  +  ay(AzBx  -  AxBz)  +  az(AxBy-  AyBx)  Eq.  4 
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where  ax,  ay  and  az  are  the  unit  vectors  in  the  x,  y,  and  z  directions.  You  can  find  all 
of  this,  of  course,  in  any  good  E&M  text  (my  favorite  being  Ramo,  Whinnery  and 
Van  Duzer/  "Fields  and  Waves  in  Communication  Electronics  ,  in  which  Sections 
2.24  and  5.23  are  particularly  relevant  here),  but  the  point  in  going  through  this 
derivation  is  that  it  is  simple  enough  that  you  don't  need  to  "look  it  up".  (I  did  this 
the  first  time  through  on  the  plane  between  Dulles  and  California  and  got  it  right.) 

What  we  are  looking  for,  as  seen  in  Figure  3,  is  the  total  field,  H(r),  due  to  the 
current,  I,  flowing  in  the  loop  of  radius,  R,  in  the  X-Y  plane,  having  as  its  axis  the  Z 
axis.  To.  get  this  field,  H(r),  we  are  going  to  have  to  sum  all  of  the  differential 
currents,  dH(r)  (from  Eq.  3)  from  all  of  the  differential  conductor  segments,  dl, 
forming  the  loop.  A  convenient  way  to  do  this  is  to  specify  the  position  around  the 

loop  by  an  angle,  0,  taken  from  the  X-axis  (and  we  pick  our  measurement  point,  r,  to 
lie  in  the  X-Z  plane,  without  any  loss  in  generality,  since  the  loop  has  circular 
symmetry).  The  length  of  the  differential  conductor  segment  associated  with  a 

differential  angle,  d0,  around  the  conductor  loop  will  just  be  dl=Rd0,  and  the 


components  of  dl  will  be 

dlx  =  -Rsin0d0,  dly  =  Rcos0d0,  and  dlz  =  0  Eq.5 

The  location,  r',  of  the  field-generation  segment,  dl,  on  the  loop  will  be  given  by 

r'x  =  R  COS0,  r'y  =  R  sin0  and  r’z  0  E*!-  ^ 

If  the  measurement  point,  r,  is  at  height  z=d  in  the  X-Z  plane  (y=0),  ut  a  distance 
r  =  X  from  the  z  axis,  then  the  components  of  r  are 

rx  =  r,  ry  =0,  and  rz  =  d  Eq.  7 

and  the  components  of  (r-r')  are 

(r-r')x  =  r-Rcos0,  (r  -  r')y  =  -R  sinG ,  and  (r-r')z  =  d  Eq.  8 

which  gives,  for  Ir-r'l^ 

I  r  -  r'  1 2  =  r^  -  2rRcos0  +  R^  +d2 

=  R2[(r/R)2  -  2(r/R)cos0  + 1  +  (d/R)2]  Eq.  9 


where  the  latter  form,  normalizing  all  dimensions  to  R,  is  the  most  convenient.  Note 
that  for  substitution  as  the  denominator  in  Eq.  3,  Eq.  9  must  be  raised  to  the  3/2 
power. 

The  numerator  of  Eq.  3  is  obtained  from  dl  X  (r  -  r'),  the  axial  (Z)  component 
of  which  is  given  from  Eqs.  4, 5  and  8  as 
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[dl  X  (r-  r’)]z  =  dlx  (r  -  r')y  -  dly  (r  -  r')x 

=  (-R  sinO  d0)(-R  sinB)  -  (R  cos8  d0)(r  -  R  cos0) 

=  R2[l-(r/R)cos0]d0  Eq.  10 

While  the  axial,  Z-component  of  field  does  all  of  the  work  in  the  cases  of  interest,  for 
nulling  out  the  crosstalk,  we  also  would  like  to  know  the  radial  (X)  component  of 
field,  which  we  get  similarly  as 


[dl  X  (r-  r’)]x  =  dly  (r  -  r')z  -  dlz  (r  -  r')y 

=  (Rcos0  d0)(d)-(O)(-Rsin0) 

=  R2[(d/R)cos0]d0  Eq.  11 

While  there  is  a  local  differential  tangential  field  component,  dHy,  for  some  values 
of  q,  from  the  circular  S5anmetry  of  the  ring,  this  component  must  integrate  to  zero, 
and  so  we  need  not  bother  with  it  here. 

From  Eqs.  10  and  11,  we  obtain  the  differential  contributions  to  the  axial, 
Hz(r,d)  and  radial  Hx(r,d)  components  of  H(r),  dH(r),  at  height  d  and  radius  r,  as 

dHz(r,  d)  =(I/47tR)[l-(r/R)cos0]/[(r/R)2 . 2(r/R)cos0  +  1  4-  (d/R)2]3/2d0  Eq.  12 

We  get  Hz(r,  d)  by  integrating  dHz(r,  d)  around  the  loop  (from  0=  0  to  2tc),  which  by 
symmetry  we  can  simplify  by  just  multiplying  by  2  and  integrating  from  0  =  0  to  7t, 

Hz(r,  d)  =  (I/27tR)  J  [l-(r/R)cos0]d0/[(r/R)2 . 2(r/R)cos0  +  1  +  (d/R)2]3/2  Eq.  13 
0 

we  can  do  the  same  for  the  radial  field  component,  Hr(r,  d)  =  Hx(r,  d)  as 


Hr(r,  d)  =  (I/27iR)  (d/R)  J  cos0  d0/[(r/R)2  -  2(r/R)cos0  +  1  +  (d/R)2]3/2  Eq.  14 

0 

While  integrating  Eq.  13  for  Hz(r,  d)  analytically  looks  a  little  daunting,  for 

the  simpler  case  of  r=0  there  is  no  0  dependence  at  all,  so  we  obtain  the  axial  field 
strength,  Hz,  versus  height,  d  on  the  axis  of  the  loop  as 

Hz(r=0,  d)  =  (I/2R)/[1  +  (d/R)2]3/2  Eq.  15 

or  for  d=0  (in  the  plane  of  the  loop,  d=0,  at  its  center,  r=0). 


Page  11 


Hz(r=0,  d=0)  =  (I/2R) 


Eq.l6 


It  is  quite  possible  that  there  are  analytical  solutions  to  the  definite  integrals  of  Eqs. 
13  and  14.  (I  didn't  have  a  table  of  integrals  on  the  plane,  and  haven’t  really  tried 
hard  to  integrate  the  functions  subsequently.)  I  know  that  Hz(r,d)  can  be  expressed 
as  the  derivative  of  complete  elliptic  integr^  functions,  but  it  seems  a  lot  simpler  to 
just  numerically  integrate  these  to  get  the  H-field.  I  found  it  quite  straightforward  to 
just  do  this  in  an  EXCEL  4.0  spreadsheet  calculation.  The  only  concern  in  carrying 
out  this  integral  numerically  is  that  it  has  a  1/r  singularity  at  r=R  for  d=0  (the 

denominator  goes  to  zero  near  6=0).  This  is  of  course  correct;  the  H-field  would 
indeed  become  infinite  if  the  conductor  radius  were  reduced  to  zero.  Practically,  it 

means  that  the  integration  step  size,  A0,  must  be  made  quite  small  near  0=0  when 
r/R  is  near  unity  and  d/R  is  near  zero  (e.g.,  near  the  surface  of  the  conductor,  when 
the  radius  of  the  conductor  is  small). 

I  have  included  two  spreadsheet  H2(r,  d)  calculations  which  calculate,  at  one 
time,  the  magnetic  field  strengths  at  numerous  radial  distances,  r,  from  the  axis  at 
the  specified  height  above  the  loop.  It  is  important  to  note  that  these  calculations  are 
for  the  external  (outside  the  metal  conductor  itself)  magnetic  field  from  a  circular 
loop  fabricated  with  circular  cross-section  wire.  Within  the  wire,  at  low  frequencies, 
the  H-field  falls  linearly  to  zero  from  the  surface  of  the  conductor  to  its  center.  At 
higher  frequencies,  the  skin  effect  forces  a  more  rapid  fall-off  of  current  with 
distance  from  the  surface.  The  Macintosh  filenames  for  the  two  Hz(r,  d)  calculation 
EXCEL  4.0  spreadsheets  are  "New  Loop  B-Field  Calc  at  d"  and  "Expanded  Loop  B- 
Field  Calc  at  d"  (in  spite  of  the  names,  they  calculate  Hz  values  which  must  be 
multiplied  by  jx  if  you  want  B  z)-  The  larger  of  the  two  spreadsheets  ("Expanded...  ) 
calculates  Hz(r,  d)  for  75  values  of  radial  distance,  r,  at  the  specified  height,  d,  and 
uses  quite  a  small  integration  step  size  to  maintain  accuracy  quite  near  r=R  for  d=0, 
so  it  is  fairly  large  (and  takes  =30  seconds  to  recalculate  on  a  Macintosh  Ilfx  after  a 
parameter  is  changed).  Figures  4,  5  and  6  show  plots  of  Hz(r/R,  d/R)  [in  ampere- 
turns  per  meter]  for  20  different  heights,  d/R,  above  an  R=1  meter  loop  for  a  current 
of  1=1  ampere  in  the  loop  (to  scale  for  other  loop  sizes  or  currents,  just  divide  H 
values  by  the  actual  loop  radius,  R,  in  meters,  and  multiply  by  the  actual  current,  I, 
in  amperes.  As  seen  in  Figure  4,  the  H-field  at  the  center  of  the  loop  (r=0)  in  the 
plane  of  the  loop  is  indeed  0.5  ampere-tums/meter,  as  indicated  by  Eq.  16.  In  the 
plane  of  the  loop,  Hz(r,  d=0)  is  fairly  flat  near  the  center,  then  increases  as  1/  (R-r) 
near  the  conductor,  reversing  sign  for  r/R>l,  with  a  negative,  -1/ (r-R),  singularity 
and  drop-off  with  increasing  r/R.  As  the  height,  d/R,  is  increased,  the  field  inside 
the  loop  drops  off  somewhat,  the  dramatic  increase  near  r/R=l  becomes  much  less 
pronounced  or  disappears  (beyond  d/R=0.4  or  0.5),  with  the  radius  at  which  the 
field  direction  becomes  negative  increasing  to  r/R=1.85  at  d/R=0.4.  Figure  5  is 
plotted  with  a  ±5  ampere-tum/meter  scale  to  show  the  field  enhancement  near 
r/R=l  better,  while  Figure  6  plots  only  the  0.9^/R<l.l  region  on  a  ±65  a-T/m  scale. 
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Calculation  of  the  Mutual  Inductance  Between  Coaxial  Loops  at  d 

In  order  to  calculate  the  inductances  for  the  practical  inductor  structures 
shown  in  Figures  1  and  1,  we  need  to  be  able  to  calculate  the  self-inductance  of  a 
loop  of  radius,  R,  and  the  mutual  inductance  between  different  sized  coaxial  loops 
displaced  from  one  another  by  some  height.  As  discussed  in  Van  Duzer  and  other 
texts,  there  are  several  methods  for  the  calculation  of  the  mutual  inductance  between 
circuits.  Since  we  already  have  the  H-field  (and  hence  the  B-field  after  multiplying 
by  jx),  we  can  simply  use  the  siurface  integral  (Stoke'  theorem)  form  where  the 
mutual  inductance,  Ml2,  (which  is  the  ratio  of  voltage  induced  in  loop  1  to  the  time 
derivative  of  current  in  loop  2)  is  given  as 

Mi2  =  (1/I2)1  ^2  •  dSi  Eq.  17 

Si 

where  the  surface  integral  is  boimded  by  the  center  of  the  conductor  forming  loop  1. 
Consider  our  case,  as  in  Figure  3,  of  a  fixed  loop  (loop  2  in  Eq.  17)  of  radius,  R2/ 
lying  in  the  X-Y  plane  and  a  second  loop,  coaxial  with  the  first,  of  radius,  Ri,  lying 
above  it  at  a  height,  d.  The  surface  integral  will  be  conducted  over  the  area  of  the 
disc  of  radius,  Rl,  enclosed  by  loop  1.  Since  this  disc  is  horizontal,  dSl  is  vertical 
(pointing  in  the  Z  direction),  so  the  only  component  of  B-field  that  will  contribute  is 
Bz  =  p.  Hz.  Hence,  for  some  height  d,  aU  we  have  to  do  is  integrate  \i  times  the 
appropriate  Hz(r,  d)  curve  in  Figure  4  from  r=0  to  r=R  l  to  get  Mi2.  The  integral  is 
over  area,  dA=r  dr,  of  course,  so  we  have 


Mi2(d)  =  (p/l2) 


J  rHz(r,d)dr 


r=0 


Eq.  18 


As  another  check  on  the  accuracy  of  the  Hz(r,  d)  spreadsheet  calculations  described 
above,  I  carried  out  this  Eq.  18  mtegral  numerically  in  an  EXCEL  spreadsheet  and 
got  excellent  agreement  with  the  textbook  complete  elliptic  integral  solutions  for 
Mi2  of  Eq.  19,  integrating  over  Hz(r,  d)  results  pasted  in  from  the  spreadsheet  used 
to  calculate  the  curves  in  Figures  4, 5  and  6. 

Because  Eq.  18  represents  a  double  integral  (because  Hz(r,  d)  comes  from  the 
integral  in  Eq.  13)  calculating  mutual  inductances  from  Bz(r,  d)  in  a  spreadsheet 
using  Eqs.  13  and  18  is  a  bit  of  a  chore.  Fortimately  much  better  mathematicians  than 
I  have  been  working  on  this  problem  for  at  least  sixty  or  seventy  years  and  have  an 
analytical  solution  in  terms  of  the  complete  elliptic  integrals  of  the  first  kind,  K(k) 
and  the  second  kind,  E(k).  These  have  been  available  in  table  form  for  many  years 
(such  as  in  Jahnke  and  Emde,  "Tables  of  Functions",  since  at  least  the  1938  edition  [I 
use  the  1945  edition],  as  well  as  Abromowitz  and  Stegun  [but  be  careful  as  they 
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define  the  modulus,  k,  differently  from  J&E  and  most  E&M  texts]).  Fortimately,  they 
are  also  defined  functions  in  engineering/scientific  software  packages  such  as  the 
excellent  "HiQ"  package  available  for  the  Macintosh,  and  some  other,  computers. 
This  makes  it  quite  easy  to  generate  the  mutual  (and  self- )  inductance  data  needed  to 
analyze  the  practical  inductor  structures  of  interest.  The  analytic  expression  for  M 
between  two  loops  is 

M  =  p  [Rl  R2]1/2  {[(2/k)  -  k]  K(k)  -  (2/k)  E(k)}  Eq.  19 
where  the  modulus  of  the  complete  elliptic  iritegral  functions,  k,  is  given  from 

k2  =  (4  Ri  R2)/[d2  +  (Ri  +  R2)2]  Eq.  20 

Listing  1  shows  a  short  (=12  line)  listing  of  HiQ  "script"  code  that  calculates  M(R  2) 
for  any  specified  d  and  Ri  values  and  stores  the  results  in  an  array  for  convenient 
graphing  or  pasting  into  EXCEL  worksheets. 

The  functional  dependencies  are  easier  to  understand  if  we  rewrite  Eqs.  19 
and  20  in  terms  of  normalized  radii  and  heights,  as  in  Eqs.  12  -  15.  If  we  take  one  of 
the  loop  radii  as  fixed  and  call  its  radius,  R,  and  call  the  other  radius,  r,  and  write  the 
expressions  in  terms  of  the  normalized  dimensions,  r/R  and  d/R,  we  have  for 
Eq.  19, 


M  =  p  R  (r/R)l/2  {[(2/k)  -  k]  K(k)  -  (2/k)  E(k)}  Eq.  21 


and 


k2  =  4  (r/R)/[(d/R)2  +  (1  +  (r/R))2]  Eq.  22 

This  form  illustrates  clearly  how  inductance  is  proportional  to  R,  assuming  all  other 
geometric  factors  remaining  in  proportion.  Values  of  M(r/R)  for  various  heights, 
d/R,  are  shown  in  Figure  7.  As  is  obvious  in  Figure  4  from  the  flatness  of  the  B-field 
near  the  center  of  the  loop  for  small  values  of  d/R,  and  the  fact  that  the  integrated 
area  in  Eq.  18  varies  as  (r/R)2  [or  R22  in  Eq.  IS],  M  must  vary  at  least  quadratically 
with  r/R  near  r/R=0.  As  seen  in  Figure  7,  this  is  indeed  the  case  for  d/R  values 
under  0.5,  and  for  d/R  values  near  zero,  the  M(r/R)  curves  climb  into  a  logarithmic 
singularity  at  r/R=l,  as  expected  in  light  of  the  l/(l-r/R)  singularity  in  the  B-field. 
Beyond  r/R=l,  the  M(r)  values  drop  off  of  the  logarithmic  singularity  and  decrease 
relatively  slowly  at  larger  r/R  values.  Very  clearly,  mutual  inductance  is  maximized 
between  equally  sized  loops  spaced  as  closely  together  as  possible.  (That  is,  of 
course,  just  what  a  solenoidal  inductor  is!) 

In  the  calculation  of  the  inductance  of  the  spiral  inductor  structure  of 
Figure  1,  we  need  to  know  the  mutual  inductances  between  coplanar  (d=0)  coaxial 
loops  of  different  radii.  While  the  M(r/R)  curves  in  Figure  7  are  drawn  over  the 
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0  <  r/R  <  2  range,  for  d/R=0  this  is  really  not  necessary.  Since  M  12  =  M2I,  we  can 
always  reverse  the  definitions  of  the  "fixed  radius,  R,  and  the  variable  radius,  r 
whenever  r  >R,  as  we  do  not  have  to  worry  about  the  d/R  scaling.  (Of  course  we 
will  have  to  rescale  the  calculated  mutual  inductance  by  multiplying  the  value  from 
a  chart  based  on  R=1  meter  by  the  new,  larger  R  value.)  All  of  the  M(r/R) 
information  is  contained  within  the  0  <  r/R  <  1  range  with  the  simple  requirement 
for  scaling  with  the  magnitude  of  the  larger  radius,  R.  Figure  8  shows  the  curve  of 
M(r/R)  for  d=0  in  the  0  <  r/R  <  1  range,  along  with  a  lOxM  expanded  curve. 
Numerical  values  from  which  these  were  plotted  are  found  in  the  EXCEL  4.0  file 
"M(r/R)  for  Spiral  Inductors  v2C". 

Figure  7  can  also  be  used  to  find  the  self-inductances  of  loops  for  either  spiral 
or  solenoid  inductor  applications.  Because  of  the  l/(R-r)  B-field  singularity  for  d=0, 
if  the  radius,  a,  of  the  conductor  forming  the  loop  were  zero,  then  the  self  inductance 
of  the  loop  would  be  infinite  (although  only  logarithmically  infinite).  Hence,  it  is 
essential  to  specify  the  conductor  radiiis,  a,  for  calculation  of  loop  self-inductance. 
As  described  in  Van  Duzer  (section  5.24)  or  other  texts,  the  external  self-inductance 
of  a  loop  is  calculated  by  carrying  out  dae  surface  integral  in  Eq.  17  to  the  interior 
edge  of  the  conductor.  I.e.,  we  calculate  the  self-inductance  for  a  circular  loop  of 
radius,  R,  by  setting  d=0  in  Eq.  18  and  carrying  out  the  integral  from  r=0  to  r=R-a 
(instead  of  r=0  to  r=R  as  for  mutual  inductance).  This  also  means  we  can,  from 
Figure  8,  pick  off  the  external  self  inductance  value,  Le,  as  the  value  for  M 
corresponding  to  the  normalized  radius,  r/R  =  (1  -  a/R).  (Of  course,  if  the  loop  size 
is  different  from  R=1  meter,  then  the  M  value  must  be  multiplied  by  the  actual  value 
of  R  [in  meters].)  This  may  be  expressed  as, 

Le  =  M(r=R-a,  d=0)  External  Self-Inductance  Eq.  23 

This  external  self-inductance  value  calculated  in  this  way  is  not  quite  the  entire  self¬ 
inductance  of  the  loop.  As  mentioned  previously,  at  low  frequencies  (i.e.,  for  5»a, 
where  5  is  the  skin  depth  in  the  conductor  metal),  the  B-field  drops  linearly  from  the 
surface  value  to  zero  at  the  conductor  center.  This  gives  a  very  low-frequency 
internal  inductance  for  a  solid  wire  of  Li  =  [l/8k  henrys/meter  of  wire  length.  Since 
the  length  of  wire  around  a  loop  of  radius  R  is  1  =  27tR,  the  very  low  frequency 
internal  self  inductance  for  a  loop  is 

(Li)o  =  (p./4)  R  Internal  Self-Inductance  (dc)  Eq.  24 

at  higher  frequencies,  the  currents  will  flow  closer  to  the  surface  of  the  inductor  and 
the  internal  inductance  will  be  smaller.  Note  that  like  the  external  self  inductance 
(Eqs.  21  and  23),  the  internal  self-inductance  also  scales  with  R.  This  internal  self¬ 
inductance,  Li,  must  be  added  to  the  external  self  inductance,  Le,  from  Eq.  23  to  get 
the  total  self-inductance  of  the  loop. 


Ls  =  Le  +  Li  Total  Self-Inductance  Eq.  25 
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In  the  modeling  of  a  solenoidal  inductor  as  a  stack  of  loops  of  equal  radius,  R, 
shown  in  Figure  2,  we  need  not  only  this  loop  self-inductance  value,  Lq/  but  also  the 
mutual  inductances,  Mij,  between  loops  separated  by  height,  d,  values  integer 
multiples  of  the  loop  stacking  pitch,  S  (or  S/R  in  normalized  form).  Figure  9  shows  a 
plot  (data  from  file  "R=r=l  Mutual  Inductance  vs  d")  of  the  mutual  inductance,  M, 
between  equal  radius  (R=r=l  meter)  coaxial  loops  vs.  the  normalized  distance,  d/R, 
separating  them.  An  expanded  lOxM  curve  is  also  included.  Hence,  Figure  9, 
together  with  Figure  8  for  the  loop  self-inductance  values,  is  all  we  need  to  calculate 
the  inductance  of  multi-turn  solenoids. 


Calculation  of  the  Inductance  of  Solenoids  by  Stacked  Loop  Method 

Figure  2  shows  the  geometry  considered  for  the  modeling  of  a  solenoid 
inductor  by  a  series  of  stacked  loops.  As  noted  above,  with  the  self-  and  mutual 
inductance  data  illustrated  in  Figures  8  and  9,  all  of  the  E&M  part  of  the  work  in 
calculating  the  inductance  for  solenoids  is  done;  what  remains  is  a  little  circuit 
arithmetic.  This  is  easier  to  see  if  we  start  from  a  simple  case:  two  turns  or  two 
coupled  loops.  Since  the  same  current  flows  through  both  loops,  and  the  sense  of 
current  direction  aroxmd  the  loops  is  the  same  (which  means  that  the  sign  of  the 
mutual  inductance  terms  will  be  positive),  the  total  voltage,  V,  induced  by  a  rate  of 
change  of  current,  dl/dt,  will  be  given  by 

V  =  L  dl/dt  =  Lis  dli/dt  +  Ml2  dl2/dt  +  L2s  dli/dt  +  M2I  dli/dt  Eq.  26 

Since  M12  =  M21  (as  seen  from  Eqs.  19  and  20,  for  example)  and  the  fact  that 
12  =  II  =1,  we  must  have  that  the  total  inductance  for  the  two  loops,  L,  is  just 

L  =  Lis  +  Mi2  +  L2s  +  M21 

=  2  Ls  +  2  Mi2  for  Lis=l2s=Ls  and  M21=Mi2  Eq.  27 

Hence,  the  calculation  of  the  solenoid  inductance  is  simply  a  matter  of  summing  up 
some  loop  self-inductances  and  mutual  inductances  between  loops  to  get  the 
answer.  This  is  very  simple.  The  only  problem  arises  as  the  number  of  turns  or 
loops,  N.  increases,  the  number  of  L  and  M  terms  to  be  looked  up  and  added  is  N^, 
which  can  get  a  bit  tedious  for  large  values  of  N. 

Fortunately  spreadsheets  are  ideal  for  doing  this  type  of  calculation  where  the 
arithmetic  is  simple  but  the  organization  of  all  of  the  data  a  bit  troublesome.  I 
created  an  EXCEL  4.0  spreadsheet  (filename:  "Any  S/R  Solenoid  L  Calc  v3")  to  carry 
out  this  task.  How  this  works  is  most  easily  seen  by  reproducing  the  example  of 
Eq.  27  for  a  somewhat  bigger  case;  we  will  consider  five  turns  or  loops.  For  N=5,  we 
have 
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L  =  Lis  +  Mi2  +  Mi3  +  Mi4  +  Mis 
M21  +  L2s  +  M23  +  M24  +  M25 
M31  +  M32  +  L3s  +  M34  +  M35 
M41  +  M42  +  M43  +  L4s  +  M45 
M51  +  M52  +  M53  +  M54  +  Lss  Eq-  28 

Notice  that  the  upper  left  (2  x  2)  comer  of  this  square  array  is  the  same  as  the  sum 
for  the  N=2  case  in  Eq.  27.  Summing  the  upper  left  (3  x  3)  corner  gives  the 
inductance  for  an  N=3  loop  solenoid,  etc.,  just  as  long  as  we  order  the  loops 
sequentially  (i.e.,  it  is  assvuned  that  we  mean,  3  "adjacent"  loops  in  a  stack  when  we 
refer  to  a  3-tum  solenoid).  Eq.  28  illustrates  that  one  mutual  inductance  array  can  be 
used  to  calculate  the  inductance  of  a  solenoid  with  the  specified  dimensions  (radius, 
R,  loop  and  turns  pitch,  S,  [entered  as  S/R])  for  any  number  of  turns,  N,  from  N=1  to 
N=[dimension  of  array]  (N=5  for  Eq.  28).  In  the  Any  S/R  Solenoid  L  Calc  v3 
spreadsheet,  the  M  array  size  is  21  x  21,  so  the  inductance  of  up  to  N=21-turn 
solenoids  can  be  calculated.  (It  is  trivial  to  extend  the  size  of  the  arrays,  if  desired,  to 
accommodate  larger  numbers  of  turns.) 

The  way  the  spreadsheet  works  is  to  first  set  up  another  square  array  that 
calculates,  from  the  relative  positions  of  any  two  loops  in  the  stack,  the  distance 
between  them,  d/R.  This  separation  distance,  d/R  array  is  used  as  the  index  to  a 
mutual  inductance  lookup  table  (the  large  data  file  R=r=l  Mutual  Inductance  vs  d 
is  included  in  the  "Any  S/R  Solenoid  L  Calc  v3"  spreadsheet  to  look  up  the  M(d/R) 
values,  since  EXCEL  doesn't,  to  my  awareness,  have  the  complete  elliptic  integral 
functions  [K(k)  and  E(k)  in  Eqs.  19  and  21]  available  in  it  to  calculate  M(r/R).  There 
is  also  an  issue  of  calculation  of  the  self-inductance  of  the  loops.  Rather  than  include 
the  entire  "M(r/R)  for  Spiral  Inductors  v2C"  file  in  this  spreadsheet  as  a  lookup  table, 
you  have  the  option  of  either  pasting  in  an  exact  Le  value  from  this  file  for  the 
specified  conductor  radius,  a/R,  (Eq.  23),  or  you  can  omit  the  number  or  put  in  zero 
for  this  "Lext"  value,  and  the  approximate  Le  value, 

Le  «  p  R  [  ln(8/(a/R))  -  2]  for  a/R«l  Eq.  29 

will  be  used  (Van  Duzer,  Section  5.24).  The  "bottom  line  for  this  Any  S/R  Solenoid 
L  Calc  v3"  spreadsheet  is  a  table  of  inductance,  L,  vs.  number  of  turns,  N,  values  for 
integer  values  of  N  from  N=1  to  N=21  (as  the  arrays  are  sized  now). 

For  the  case  of  the  solenoid  inductor,  there  are  a  number  of  reasonably  good 
equations  which  give  at  least  pretty  good  approximations  to  the  inductance  to  an 
air-core  solenoid  over  a  fair  range  of  design  parameters.  The  simplest,  fairly 
accurate,  formula  (appearing  in  F.  E.  Terman,  "Radio  Engineering  Handbook"  for 
1943  and  most  other  handbooks  and  texts)  is 

L  =  (n  p  N2  R  )/[0.9  +  1/R]  =  (7C  p  N2  R  )/[0.9  +  N  (S/R)]  Eq.  30 
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where  in  the  second  version,  the  nominal  length  of  the  inductor,  1,  has  been  replaced 
by  1  =  N  S  ,  where  S  is  the  turns  pitch  as  shown  in  Figure  2.  A  much  more 
sophisticated  formula  is  given  in  Jahnke  and  Emde  (p.  87,  with  tables  on  pp.  88-89) 

L  =  (2/3)  p  R  N2  { [K(k)  +  ((d/l)2  -l)E(k)]/sina  -  (d/l)2}  Eq.  31 

where  d  =  2  R  is  the  diameter,  1  the  length,  and  a  is  defined  by  tana  =  d/1.  The 
fact  that  this  formula  for  the  inductance  of  thin-walled  solenoids  is  based  on  the  K(k) 
and  E(k)  complete  elliptic  integral  fimctions  (like  Eqs.  19-22)  gives  on  a  certain 
degree  of  confidence  in  its  rigorousness.  Listing  2  shows  a  short  listing  of  HiQ 
"script"  code  that  calculates  L  vs.  N  from  Eq.  31,  again  taking  1  =  N  S  for  the  length. 
Both  Eqs.  30  and  31  were  used  for  comparison  with  the  "Any  S/R  Solenoid  L  Calc 
v3"  spreadsheet  calculation.  As  seen  in  Figures  10-13,  in  general  the  results  were 
quite  close,  particularly  between  my  "stacked  loop"  spreadsheet  calculation  and  the 
J&E  formula  (Eq.  31).  For  S/R=0.2  in  Fig.  10,  the  inductance  numbers  from  both  the 
"Terman  Formula"  (Eq.  30)  and  the  "Jahnke  and  Emde"  formula  (Eq.  31)  are  nearly 
identical  with  the  results  from  this  "Stacked  Loop  Calc"  spreadsheet  calculation.  As 
seen  for  S/R=0.025  (an  extremely  fine  pitch;  S/R=0.2  with  R=0.0006  meters  is 
thought  to  be  about  the  limit  for  simple  wirebond  embedded  inductors  in  MCM’s)  in 
Figiire  12b,  the  percentage  differences  tend  to  be  larger  for  fewer  numbers  of  turns, 
and  the  agreement  is  better  with  the  J&E  formula  (Eq.  31).  This  is  clear  on  a  log  plot 
like  Figure  12b,  but  not  on  the  linear  plots  like  Figure  12a.  Note  that  for  a  fairly  wide 
turns  spacing,  S/R=0.4  in  Figure  13,  this  "Stacked  Loop  Calc"  spreadsheet 
calculation  gives  about  a  10%  higher  L  (for  a/R=0.03)  than  Eqs.  30  or  31.  Since  with 
larger  spacings  the  self-inductance  assiunes  higher  importance  (because  the  mutual 
inductances  are  relatively  smaller),  my  guess  would  be  that  the  stacked  loop 
calculation  is  probably  closer  to  reality. 

The  results  for  solenoid  inductance  calculations  using  this  "Any  S/R  Solenoid 
L  Calc  v3"  spreadsheet,  L  vs.  N  for  integer  values  of  N  from  1  to  21,  for  nine  different 
S/R  values  are  plotted  on  Figures  14  (linear  scale)  and  15  (log  scale),  with  numerical 
data  presented  in  Table  1.  It  should  be  pointed  out  that  there  is  much  more  data 
available  in  the  calculation  than  is  presented  in  this  summary.  In  particular,  the 
mutual  inductance  array  is  available  (e.g.,  Eq.  28),  such  that  the  details  of  the 
coupling  are  known  on  a  tum-by-tum  basis.  This  would  be  invaluable,  for  example, 
if  the  turns  were  split  off  into  two  or  more  different  circuits;  the  M  array  would 
provide  precise  information  for  the  transformer  design  and  characterization. 


Calculation  of  L  for  Planar  Spiral  Inductors  by  Stacked  Loop  Method 

Figure  1  shows  the  geometry  considered  for  the  modeling  of  a  planar  spiral 
inductor  by  a  series  of  nested  loops.  Because  all  of  the  loops  lie  in  the  same  plane,  all 
of  the  E&M  part  of  the  work  in  calculating  the  inductance  that  we  need  is 
summarized  in  Figure  8  alone.  Basically,  the  circuit  arithmetic  is  the  same  as 
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described  above  (Eqs.  26-28)  for  the  solenoid  case,  except  that  since  the  loop  radii  are 
all  different,  the  Ls  values  will  aU  be  different  and  we  have  to  remember  to  scale  the 
inductance  values  looked  up  out  of  an  R=1  meter  table  by  multiplying  them  by  the 
actual  R  [in  meters].  The  scaling  for  the  planar  spiral  inductor  is  based  on  the  radius 
of  the  center  of  the  conductor  of  the  outermost  (N=l)  loop.  The  loops  (turns)  are 
numbered  consecutively  from  the  outside  inward.  While  the  spreadsheet  calculation 
(fliename:  "Spiral  Inductor  L(N,S/R)  CalcVl")  is  capable  of  supporting  irregular 
turns  pitches,  in  the  calculation  the  turns  pitch,  S,  is  taken  as  constant  (entered  in 
normalized  form  as  S/R).  Hence,  the  N=1  loop  has  radius  R,  the  N— 2  turn  has  r— R-S, 
N=3  has  r=R-2S,  etc.  In  general  the  radius  of  the  Nth  turn  of  the  spiral  is 

r(N)  =  R[l-(N-l)(S/R)]  Eq.32 


Obviously,  since  r(N)  cannot  be  negative,  for  any  given  value  of  S/R,  there  is  only 
room  for  a  certain  number  of  turns  in  the  spiral.  It  turns  out  that  keeping  track  of 
that  fact  is  the  biggest  factor  in  making  the  "Spiral  Inductor  L(N,S/R)  CalcVl" 
spreadsheet  calculation  a  lot  more  complicated  than  for  the  solenoid  case.  However, 
the  approach  is  identical,  the  goal  being  to  build  up  the  inductance  array  like  Eq.  28. 
The  "M(r/R)  for  Spiral  Inductors  v2'*  data  used  to  plot  Figure  8  is  included  in  the 
spreadsheet  as  a  lookup  table  for  the  loop-to-loop  mutual  inductance  values  (and  for 
the  loop  self-inductance  values  as  weU,  using  r=R-a  [Eq.  23]). 


The  results,  including  both  the  L(N)  data  and  the  inductance  array  (M  array), 
from  this  "Spiral  Inductor  L(N,S/R)  CalcVl"  spreadsheet  calculation  for  S/R  values 
of  0.5,  0.4,  0.35,  0.3,  0.25,  0.2,  0.175,  0.15,  0.125,  0.11,  0.10,  0.09,  0.08,  0.07,  0.06,  0.05, 
0.04, 0.03, 0.02  and  S/R=0.01 ,  all  with  conductor  radius  a=0.4(S/R),  are  summarized 
in  the  file  "R=lm  Spiral  Inductor  L(N,S/R).S".  This  is  important  because  it  contains 
the  mutual  inductance  array  (e.g.,  Eq.  28)  for  each  case,  such  that  the  details  of  the 
coupling  are  known  on  a  tum-by-tum  basis.  This  would  be  invaluable,  for  example, 
if  the  turns  were  split  off  into  two  or  more  different  circuits;  the  M  array  would 
provide  precise  information  for  the  spiral  transformer  design  and  characterization. 


The  problems  with  planar  spiral  inductors  are  fairly  apparent  from  the  curves 
in  Figure  16.  Unless  the  S/R  value  is  extremely  small  (and  such  small  S/R  values 
lead  to  very  small  conductor  dimensions,  high  resistance  and  generally  degraded  Q 
values),  the  difference  in  radius  between  the  loops  leads  to  poor  coupling  (low 
mutual  inductance)  between  turns.  This,  and  the  smaller  sizes  of  the  loops  as  their 
radii  decrease  as  N  increases,  leads  to  the  rapid  saturation  of  the  increase  of  L  with 
N  in  spirals.  This  is  more  apparent  in  the  log-log  plot  in  Figure  17  (where  the  slope  is 
the  power  of  N  with  which  L(N)  is  increasing).  Contrast  the  curves  for  practical 
values  of  S/R  (e.g.,  S/R<0.2  for  low  resistance)  for  the  spirals  in  Figure  17  with  the 
equivalent  log-log  plot  for  the  solenoid  inductors  with  various  S/R  values  in 
Figure  18.  The  planar  spiral  is  convenient  for  fabrication  as  an  embedded  MCM 
inductor.  Unfortunately  that  is  about  all  the  good  that  can  be  said  of  it,  since  it  is 
electrically  substantially  inferior  to  the  solenoid  inductor  structure. 
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Dependence  of  Bandwidth,  Q  and  Noise  on  Rac  of  Inductor 

The  prime  specification  for  an  embedded  MCM  inductor,  next  to  that  of  its 
inductance,  L,  is  its  Q  at  a  given  frequency,  f,  which  is  defined  as  n  times  the 
maximum  energy  stored  divided  by  the  energy  dissipated  per  half  cycle. 

Q(f)  =  K  (maximum  energy  stored)/ (energy  dissipated  per  cycle)  Eq.  33 

With  a  sinewave  current, 

i(t)  =  ip  sin(27ift)  Eq*  34 

through  an  inductor  having  an  inductance,  L,  and  an  ac  series  resistance,  Rac/ 

maximum  energy  stored  =  (1/2)  L  ip^  Eq-  35 

and  the  average  power  dissipation  in  the  series  resistance,  Rac,  is 

average  power  dissipation  =  Rac  Irms^  =  Bac  (ip/ '^^9'  -  (1/^)  E^ac  ip^  Eq.  36 

so  that  the  energy  dissipated  per  half  cycle  (i.e.,  in  a  period  of  time  t=l/2f),  is  just 

energy  dissipated  per  half  cycle  =  average  power  dissipation/ (2f) 

=  1(1/2)  Rac  ip^]/  (2f)  Eq.  37 

so  that  from  Eqs.  33, 35  and  37,  the  inductor  Q  at  the  signal  frequency,  f,  is  given  as 

Q(f)  =  7C  [(1/2)  L  ip2]/{[(l/2)  Rac  ip2]/(2f)  =  27:f  L/Rac  Eq.  38 

Since  the  angular  frequency  (radians  per  second)  is  just  co=27tf,  this  is  frequently 
written  equivalently  as 


Q  =  2Kf  L/Rac  “  ®  L/Rac 

The  factor  of  it  in  the  definition  of  the  Q  of  an  inductor  (Eq.  33)  comes  from 
the  original  definition  of  Q  in  terms  of  the  bandwidth  of  a  simple  resonant  circuit.  If 
the  inductor,  L,  is  connected  in  parallel  with  a  capacitance,  C,  for  Rac=0  (or  Q=~), 

the  circriit  will  resonate  at  a  frequency  (Oq  =  27ifo  given  by 

©o  =  2jifo  =  l/V(LC)  Eq.40 

For  finite  loss  (finite  Q),  Rac^iO,  the  series  RLC  circuit  resonates  at  a  frequency,  cor 
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cor  =  27cfr  =  coq  [1  -  (1/2Q)2]0.5  Eq.  41 

where  the  finite  Q  resonant  frequency,  ©r,  is  very  slightly  lower  than  the  lossless 

(Q=oo)  resonant  frequency,  ©o  (s-S*/  0.0078%  lower  for  Q=40  or  0.125%  lower  for 
Q=10).  If  this  RLC  resonant  circuit  is  driven  by  an  ac  current  source,  the  voltage  will 
reach  a  peak  value  at  the  frequency  fr/  and  drop  to  a  voltage  of  1/V2=0.7071  times 
this  peak  value  at  frequencies  Af  above  and  below  this  (the  half-power  points), 
where,  for  reasonably  high  values  of  Q,  Af  is  given  by 

Af  =  fr/(2Q)  Eq.42 

Using  the  usual  definition  of  bandwidth  as  the  frequency  difference  between  the 
half-power  points  (or  about  -3dB  voltage  points),  or  BW=2AF  in  this  case,  we  have, 

Q  =  fr/2Af  =  (resonance  frequency) /(bandwidth)3db  Eq-  43 

which  was  the  original  definition  of  Q  used  in  rf  circuit  work.  Hence,  the  factor  of  7i 
in  the  energy  definition  of  Q  (Eq.  33)  is  required  for  compatibility  with  the  earlier  rf 
circuits  defiriition  (Eq,  43). 

In  this  "reasonably  low  loss"  (e.g.,  Q  above  5  or  10)  case,  the  parallel  resonant 
circuit  (capacitor  C  in  parallel  with  inductor  L,  with  its  series  resistance,  Rac) 
resonance  looks  essentially  like  a  pure  resistance,  Rp,  whose  value  is  given  by 

Rp  =  Q2Rac  atf  =  fr  Eq.44 

If  Eq.  44  at  first  glance  males  it  appear  that  increasing  Rac  would  increase  Rp,  this  is 
wrong,  as  increasing  Rac  decreases  Q  (Eq.  39),  so,  due  to  the  Q2  term  in  Eq.  44, 
Rp  decreases  with  increasing  Rac-  This  may  be  clearer  when  the  expression  for  Rp  at 
resonance  is  written  differently, 

Rp»l/[(©C)2Rac]=  IZcl2/Rac  at©=©r  Eq.45 

where  IZcl=l/©C  is  the  magnitude  of  the  capacitor  impedance  at  ©=C0r,  or 
equivalently  as 


Rp «  (©L)2/Rac  =  \Z\\ 2/Rac  at ©=©r  Eq. 46 

where  I  Zc !  =  ©E  is  the  magnitude  of  the  inductor  impedance  at  ©=©r.  Both  of  these 
forms  show  that  Rp  indeed  is  proportional  to  the  inverse  of  Rac,  so  that  for  highest 
impedance  in  the  resonant  circuit,  the  ac  resistance,  Rac  must  be  kept  as  small  as 
possible.  Also,  it  can  be  shown  that  the  equivalent  thermal  noise  current  (Johnson 
noise)  for  the  resonant  circuit  is  the  same  as  that  calculated  for  a  resistor  of  value  Rp. 
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Magnetic  Field  Intensity/Current  Distribution  on  Solenoid  Conductor 

As  seen  in  the  previous  section,  the  performance  quality  of  an  inductor,  in 
terms  of  attainable  impedance,  bandwidth,  Q,  noise,  efficiency,  etc.  is  critically 
related  to  just  how  low  the  effective  series  resistance  of  the  inductor,  Rac/  can  be 
made.  As  discussed  in  my  report  of  6/1/96,  external  factors  such  as  dielectric  losses 
in  substrate  materials  can  lower  the  Q  values  of  MCM  embedded  inductors  if 
substrates  like  silicon  are  used,  normally  the  limiting  factor  for  Rac  is  the  skin  effect 
losses  in  the  conductors.  (For  ferrite  or  other  ferromagnetic  core  inductors  with 
substantial  relative  permeabilities  [jir  »  1]/  tbe  core  material  losses  can  dominate, 
but  we  will  be  considering  principally  the  case  of  non-magnetic  [air  or  other  pr  =1] 
core  inductors  here.) 

Calculation  of  the  dc  resistance,  Rdc/  of  a  metal  conductor  having  resistivity, 
p,  is  quite  simple,  since  for  most  geometries  the  resistance  per  unit  length,  Rdc/l  is 
given  from  the  cross-sectional  area,  Ac,  as 

Rdc/l  =  P /Ac  Eq.47 


(The  exception  is  geometries  for  which  the  direction  of  current  flow  or  cross- 
sectional  area  of  the  conductor  changes  very  rapidly  along  the  current  path,  such  as 
at  sharp  comers  in  conductors  or  abrupt  changes  in  conductor  size,  but  these  are  not 
serious  issues  in  most  inductor  geometries.)  For  example,  the  dc  resistance  of  a 
circular  wire  of  radius  a  or  diameter  dw=2a  will  be  given  by 

Rdc/l  =  p/in  a2)  =  4p/(7t  dw^)  for  Circular  Wire  Eq.  48 

Unfortunately,  things  are  not  so  simple  for  the  case  of  the  ac  skin-effect 
resistance  of  conductors  in  inductors.  There  are  in  fact  expressions  for  the  ac  skin- 
effect  resistance,  Rac/  of  certain  specific  conductor  geometries,  such  as  that  of  a 
circular  wire  used  as  the  center  conductor  of  a  coaxial  cable,  but  if  the  same  wire  is 
wrapped  with  very  close  turns  spacing  around  a  core  to  make  a  solenoid  inductor, 
this  coax  line  expression  for  Rac  wiU  not  give  the  correct  value  of  the  ac  resistance  of 
the  wire.  The  problem  is  that  of  "current  crowding",  or  non-uniform  distribution  of 
ac  current  aroimd  the  conductor  periphery.  Because  of  the  circular  symmetry  of  a 
coaxial  line,  the  current  flows  uniformly  around  the  periphery  of  the  center 
conductor  wire,  even  as  the  skin  effect  forces  the  current  to  flow  nearer  and  nearer  to 
the  surface  as  frequency  is  increased.  This  uniformity  of  the  skin-effect  current 
around  the  center  conductor  of  a  coaxial  line  is  a  direct  reflection  of  the  uniformity 
of  the  tangential  magnetic  field  around  the  surface  of  the  conductor.  (This  is  just  a 
specific  instance  of  the  general  magnetics  principal  that  the  change  in  tangential 
magnetic  field,  Ht  (in  amperes  per  meter)  at  the  surface  of  a  conductor  equals  the 
surface  current  density,  Js,  flowing  in  the  conductor. 


Page  22 


In  the  case  of  a  solenoid,  the  magnetic  field  at  the  surface  of  the  conductors  is 
strongly  influenced  by  the  overall  form  of  the  inductor,  not  just  by  the  local  shape  of 
the  conductors  as  in  a  coaxial  line.  This  could  have  important  implications.  For 
example,  in  the  idealized  view  of  a  solenoid  inductor,  the  magnetic  field  is  strong 
inside  the  solenoid,  and  essentially  zero  outside.  If  that  is  the  case,  then  the  high- 
frequency  surface  current  on  the  conductors  (since  the  surface  current  density  must 
equal  the  tangential  H-field  at  the  surface)  would  be  near  zero  on  the  outside  of  the 
conductor,  and  hence  the  conductor  loss  would  be  essentially  twice  as  large  as  you 
might  otherwise  think,  since  with  this  "current  crowding",  only  half  of  the  periphery 
contributes  to  conduction.  Also,  this  would  mean  that  the  substrate  metal  return 
path  in  a  wirebond  inductor  might  as  well  be  made  only  half  as  many  skin  depths 
thick  as  would  ordinarily  be  specified,  since  the  bottom  side  wouldn’t  contribute  to 
high-frequency  conduction  anyway.  (As  I  will  show  later,  this  turns  out  not  to  be  the 
case  for  geometries  of  interest.)  In  any  event,  in  order  to  calculate  Rac  for  the 
conductors  in  a  solenoid,  or  other  inductor  structures,  we  must  first  examine  the 
magnetic  field  distribution  around  the  conductors,  in  order  to  determine  what 
expression  for  Rac  should  be  suitable  for  use  in  that  geometry. 

The  geometry  for  the  calculation  of  the  solenoid  magnetic  field  was  shown 
previously  in  Figure  2,  and  the  stacked-loop  calculational  approach  in  which  the 
magnetic  field  at  any  point  in  the  solenoid  is  obtained  as  the  sum  of  the  fields  from 
each  individual  loop  at  that  point  (Figures  3-6).  Note  that  in  these  H(r/R)  plots,  what 
is  plotted  is  the  z-component  of  magnetic  field,  H7(r/R),  versus  radial  distance  from 
the  axis  of  the  solenoid,  measured  in  the  plane  of  symmetry  normal  to  the  axis 
passing  through  the  center  of  the  center  loop  of  the  solenoid.  This  is  convenient  in 
that  in  this  plane  of  symmetry,  H(r/R)=Hz(r/R)  since  Hr=0  by  symmetry.  This 
means,  of  course,  that  only  odd  numbers  of  loops  (turns)  in  the  solenoid  can  be 
treated,  but  this  plane  gives  the  worst-case  field  differences  between  the  inside  and 
outside  of  the  conductors.  Figure  19a  shows  Hz(r/R)  for  a  turns  pitch,  S,  to  loop 
radius,  R,  ratio  of  S/R=0.2  (a  comfortable  value  even  for  modest-sized  wirebond 
loops)  for  r/R=0  (axis  of  solenoid)  to  r/R=2  (outside  of  the  solenoid  loops  by  one 
loop  radius)  for  solenoids  of  N=l,  3,  5,  7,  9, 11  and  13  turns.  Note  that  initially  the 
field  on  the  axis,  Hz(r=0),  is  nearly  proportional  to  the  number  of  turns,  indicating 
reasonably  strong  coupling  between  1st  and  2nd  nearest  neighbor  loops,  but  the 
increase  in  Hz(r=0)  is  small  between,  for  example,  N=ll  and  13  turns  because  of  the 
substantial  axial  distance  to  these  end  turns  (1.2R).  On  the  other  hand,  as  the  total 
length  of  the  solenoid  becomes  large  (e.g.,  2.4R  for  N=13  turns  with  S/R=0.2),  the 
field  distribution  within  the  solenoid  (r/R<l)  becomes  much  flatter,  and  much  larger 
than  the  field  outside.  For  example,  for  N=13  turns,  Hz=3.97  ampere-turns /meter 
(A/m)  at  r=0,  increasing  to  4.21  A/m  at  r=0.8R,  while  outside  the  solenoid  at  r=1.2R 
the  field  is  only  Hz=-0.6  A/m. 

While  the  gross  solenoid  field  in  Figure  19a,  particularly  for  larger  numbers  of 
turns,  is  much  higher  inside  the  loops  than  it  is  outside,  to  understand  the 
implications  of  the  field  distribution  on  the  ac  skin-effect  resistance  of  the 
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conductors,  we  must  look  in  detail  at  the  fields  immediately  adjacent  to  the  surface 
of  the  wires.  The  field  in  this  region  (near  r=R)  is  shown  in  detail  in  Figure  19b. 
Remember  that  the  stacked-loop  calculation  of  Hz(r)  for  the  solenoids  in  Figures  19- 
25  assume  an  extremely  small  conductor  radius  for  the  loops,  which  leads  to  the 
near-singularity  in  Hz(r)  at  r=R.  This  is  just  fine,  since  it  is  understood  that  for  a 
given  actual  conductor  radius,  a,  the  external  Hz(r)  field  shown  in  the  Figures  is 
correct  for  the  region  outside  the  conductor  (i.e.,  for  r  <  (R-a)  and  r  >  (R+a)).  At  dc, 
the  Hz-field  within  the  metal  conductors  drops  linearly  from  the  value  at  the  surface 
(r=(R-a)  or  r=(R+a))  to  zero  at  the  center  (r=R)  (this  is  for  a  straight  circular  wire;  it 
will  be  slightly  off-center  for  a  solenoid).  For  higher  frequencies,  the  dropoff  of  the 
H-field  as  one  goes  away  from  the  surface  is  more  rapid,  determined  by  the  skin 

depth,  5.  It  is  the  details  of  this  dropoff  in  H-field  within  the  conductor  (since  it 
mirrors  the  surface  current  density  within  the  conductor)  that  determines  how  the  ac 
resistance  is  to  be  calculated.  Note  that  because  of  the  l/(R-r)  dependence  of  the  H- 
field  near  the  conductor,  the  vertical  scale  factor  has  been  compressed  from  ±8  A/ m 
in  Figure  19a  to  ±20  A/m  in  Figure  19b.  At  this  scale  factor,  it  is  clear  in  Figure  19b 
that  at  small  distances  from  the  center  of  the  loop  conductor  (r/R  near  1),  the 
dominant  effect  is  that  of  the  l/(R-r)  behavior  of  the  H-field  around  circular  wires, 
as  opposed  to  the  solenoid  shape  effects.  For  example,  even  for  a  fairly  large 
diameter  conductor,  a/R=0.03,  the  surface  H-field  (and  correspondingly,  ihe  surface 
current  on  the  conductor,  Js)  at  its  highest  point  (r=R-a  or  r/R=0.97  in  Figure  19b) 
varies  from  5.76  to  7.57  A/m  for  N=1  to  13  turns,  while  the  field  at  the  lowest  point 
on  the  outside  surface  varies  from  -4.87  to  -3.78  A/m,  so  that  the  variation  (Eq.  52)  is 
from  ±8.4%  at  N=1  to  ±33%  at  N=13,  with  a  variation  of  ±21%  for  N=5,  a  more 
typical  number  of  turns.  Since  these  numbers  represent  the  extremes  in  surface 
current  density  over  the  conductor  wire  surfaces,  simply  using  the  same  expressions 
for  Rac/  the  ac  skin  effect  resistance,  as  for  a  straight  wire  (such  as  the  center 
conductor  of  a  coaxial  line)  should  be  quite  accurate. 

For  larger  values  of  the  solenoid  turns  pitch,  S,  such  as  the  S/R=0,4  case 
shown  in  Figure  20a  and  20b,  this  trend  is  even  clearer.  The  gross  Hz(r)  field  within 
the  solenoid  is  quite  uniform  (Figure  20a  for  r/R<0.8)  for  larger  numbers  of  turns 
because  of  the  greater  length  of  the  solenoid  (4.8R  for  N=13  turns),  but  the  field 
strength  at  the  center  is  lower  than  in  Figure  19a  because  of  the  reduced  coupling 
between  the  more  widely  spaced  turns.  This  also  means  that  the  fields  at  the  surface 
of  the  conductors  is  even  more  uniform  for  S/R=0.4  in  Figure  20b  than  for  S/R=0.2 
in  Figure  19b,  ±23%  at  N=13  or  ±17.5%  at  N=5  (at  N=1  it  is  the  same  as  in  Figure  19b, 
of  course).  Hence  the  use  of  a  straight-wire  formula  for  Rac  is  even  better  justified 
for  this  S/R=0.4  case. 

In  the  opposite  direction,  to  solenoids  with  tighter  (finer)  wiring  pitches, 
nCHIP  has  fabricated  wirebond  solenoids  (details  of  which  will  be  discussed  later) 
which  have  loop  areas  equivalent  to  that  of  an  R=0.00107  meter  radius  circle  with  a 
turns  pitch  of  S=98.2  pm  using  1  mil  diameter  wire,  and  it  appears  reasonable  to 
achieve  a  pitch  of  S=125  pm  with  1.5  mil  diameter  wire  (radius,  a=19.05  pm),  which 
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corresponds  to  S/R=0.1169  and  a/R=0.0178.  Hence,  the  S/R=0.1  case  illustrated  in 
Figures  21a  and  21b  is  probably  representative  of  the  practical  limit  for  fine  pitch  in 
MCM  wirebond  solenoids.  (nCHIP’s  1  mil  wire  example  has  S/R=0.0918,  but  the  Q 
is  poorer  because  of  the  sm^er  wire  size,  and  the  difference  in  surface  fields,  which 
depends  principally  on  a/S  is  worse  for  the  1.5  mil  wire  on  125  ^im  pitch  case  than 
for  the  1.0  mil  wire  on  98.2  \Lm  pitch.)  We  note  in  Figure  21a  the  tight  coupling 
between  turns  as  evidenced  by  the  nearly  linear  increase  of  Hz(r=0)  with  N  up  to  7 
to  9  turns  (curves  for  N=l,  3,  5,  7,  9,  11,  13, 15, 17,  19  and  21  turns  are  shown  in 
Figures  21a  and  21b).  As  would  be  expected,  the  inside/outside  field  variation  is 
larger  for  this  fine  turns  pitch  case.  However  for  a  wirebond  conductor  radius  of 
a/R=0.0178,  at  N=5  turns  the  field  uniformity  is  ±15%,  and  even  for  N=21  turns,  it  is 
±34%. 


Inasmuch  as  a  ±34%  variation  in  surface  current  density  over  the  periphery  of 
the  conductor  wire  is  significant,  it  is  worth  calculating  how  much  such  a  non 
uniformity  in  current  density  would  be  expected  to  perturb  the  accuracy  of  a 
calculation  of  the  wire  high-frequency  resistance,  Rac-  If  we  make  the  assumption 
that  the  surface  current  density  aroimd  the  conductor  is  sinusoidal,  that  is, 

Js(<l))=Ja  +  Jpsin((j))  Eq.49 

where  Ja  is  the  average  value  of  surface  current  around  the  wire  and  Jp/Ja 
represents  the  field  uniformity  figure  as  quoted  above,  while  ({)  denotes  the  angular 
position  aroxmd  the  wire  (with  (j)=0  directed  parallel  to  the  solenoid  axis,  (j)=7i/2  the 

peak  field  point  on  the  inside  of  the  loop  and  ^-3n/2  the  minimum  field  point  on 
the  outside  if  the  loop).  The  use  of  a  straight-wire  value  for  Rac  would  be  assuming 
that  the  value  of  Js  was  everywhere  the  same  as  its  average  value,  Ja.  Because  power 
dissipation  varies  as  Js^/  the  non  imiform  current  distribution  in  Eq.  49  will  give  a 
higher  power  dissipation,  and  hence  a  higher  value  for  Rac/  than  the  xmiform 
straight-wire  case.  The  difference  between  the  actual  solenoid  Rac  value,  (Rac)sol/ 
and  the  straight-wire  calculated  Rac  value,  (Rac)coax/  will  be  given  by 

(Rac)sol/(Rac)coax  =  [Average0s(«2)]/[Averageas(c|)))]2 

=  1(1/211)  J  [Js(«)2  di).  ]/l  Ja  P  Eq.  50 

0 

which,  substituting  Eq.  49  for  Js((|))  and  carrying  out  the  integral  gives 

(Rac)sol/ (Rac)coax  =  [  Ja^  +  0-5  Jp^]/  Qa^) 

=  1+0.5  Gp/Ja)^  Circular  Conductor  Eq.  51 
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where  Jp/Ja  is  the  same  as  the  field  uniformity  values  cited  above.  Specifically,  for  a 
wire  radius,  a,  and  solenoid  loop  radius,  R,  the  value  of  Jp/Ja  is  given  by 

Jp/Ja  =  [Hz(R-a)  -  Hz(r+a)]/[[Hz(R-a)  +  Hz(r+a)]  Eq.  52 

(where  magnitudes,  rather  than  signed  quantities,  have  been  assumed  for  the 
Hz(R-a)  and  Hz(r+a)  values).  For  the  S/R=0.1,  a/R=0.0178  wirebond  solenoid 
example  cited  above  (Figure  21b),  for  N=5  we  foirnd  Jp/Ja  =  15%,  which,  from  Eq.  51 
gives  (Rac)sol/(Rac)coax  =  l-OHS  (which  means  only  a  1.15%  error  assuming  a 
straight-wire  loss  expression),  or  for  the  extreme  N=22  turns  case,  for  which  we 
found  Jp/Ja  =  33.7%,  Eq.  51  gives  (Rac)5ol/(Rac)coax  =  1-0567,  or  a  5.67%  error  in 
using  a  straight-wire  calculation  of  Rac  lor  ^  solenoid.  This  is  a  very  reasonable  error 
for  an  inductor  Q  calculation,  and  the  correction  factor  from  Eq.  51  can  always  be 
used  to  correct  the  solenoid  Rac  values  calculated  using  a  straight  circular  conductor 

expression  for  Rac- 

Figures  22a  and  22b  (for  S/R=0.05)  and  23a  and  23b  (for  S/R=0.03)  are 
illustrative  of  what  more  conventional  (tightly  woimd,  S=2a,  insulated  wire  single¬ 
layer)  solenoids  look  like.  Inasmuch  as  there  appears  no  practical  way  to  fabricate 
such  fine  pitches  in  a  wirebond  solenoid  for  an  MCM,  these  are  included  here  for 
reference  only.  The  S/R=0.03  case  of  Figures  23a  and  23b  does,  however  serve  as  an 
excellent  starting  point  to  examine  another  important  problem;  the  use  of  stripe  or 
ribbon  elements  in  fabricating  mixed-signal  MCM  embedded  inductors.  One 
example  of  this  would  be  the  use  of  ribbon  bonds  in  place  of  ordinary  round  wire 
wirebond  loop  conductors.  Another  is  the  need  to  analyze  the  Rac  for  metal 
microstrip  return  paths  for  wirebond  inductors  (in  which  the  wirebonds  form  the 
top  half  of  the  solenoid  loops  while  substrate  metallization  stripes  form  the  bottom 
return  path.  We  would  like  to  have  some  insight  into  the  fields  associated  with  these 
flat  conductor  structures  similar  to  what  is  shown  in  Figures  19-23  for  circular 
conductors  in  solenoids.  Such  an  example  was  carried  out  using  the  stacked-loop 
approach.  The  geometry  for  this  calculation  is  shown  in  Figure  24.  Here,  a  ribbon 
conductor  of  dimensions  0.03R  x  0.15R  is  simulated  with  a  stack  of  5  circular  cross- 
section  conductors,  each  carrying  a  current  of  0.2  A  (to  give  1  A  total  conductor 
current,  as  in  the  previous  examples).  These  ribbon  conductors  (e.g.,  7  loops  as 
shown  in  Figure  24)  are  stacked  on  a  pitch  of  S/R=0.21  to  complete  the  solenoid. 
Figures  25a  and  25b  show  the  Hz(r/R)  plots  for  N=1  to  N=7  turns  (ribbon  conductor 
loops)  of  the  0.03R  x  0.15R  conductors.  Note  that  the  gross  fields  for  the  ribbon 
solenoid  in  Figure  25a  look  virtually  the  same  as  for  the  S/R=0.2  circular  conductor 
case  in  Figure  19a.  The  more  detailed  field  plots  near  the  conductors  show  much 
lower  fields  near  the  ribbon  conductor  surface  in  Figure  25b  than  for  the  circular 
conductor  in  Figure  19b.  This  is  expected,  of  course,  because  of  the  width  of  the 
ribbon,  which  avoids  the  strong  l/(R-r)  dependence  of  Hz(r)  near  the  conductor). 
For  example,  at  the  inner  surface  of  the  ribbon  conductor,  r=0.985  in  Figure  25a,  for 
N=1  turn,  Hz(0.985)=3.596  A/m,  as  compared  to  Hz(0.985)=ll.ll  A/m  for  the  round 
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conductor  in  Fig.  19b  (3. lx  higher),  while  for  N— 7  turns,  H2(0.985)— 4.72  A/m  for  the 
ribbon  solenoid  vs.  Hz(0.985)=14.36  A/m  for  the  circular  wire.  The  field  at  the 
outside  solenoid  surface  is  Hz(r/R=1.055)=-2.85  A/m  for  N=1  turn  and  -2.06  A/m 
for  N=7  turns  for  the  ribbon  solenoid  (as  compared  to  -10.1  A/m  for  N=1  and  -9.3 
A/m  for  N=7  with  circular  wire).  This  gives  a  current  uniformity  from  Eq.  52  for  the 
ribbon  solenoid  of  Jp/Ja=±11.6%  for  N=l,  increasing  to  J  p/Ja=±39.2%  for  N=7.  (The 
circular  wire  was  more  uniform  at  4.71%  at  N=1  and  21.4%  for  N=7  turns.) 

If  the  field  distribution  around  the  periphery  of  the  ribbon  conductor  varied 
sinusoidally  as  assumed  for  the  circular  conductors  in  Eqs.  49-51,  then  the  correction 
factor  for  the  ac  resistance  would  be  (Rac)sol/ (Rac)coax  =  1.0067  for  N=1  for  the 
ribbon  loop  and  (Rac)sol/(Rac)coax  =  1*0769  for  the  N=7  turn  solenoid.  If  the  surface 
current  distribution  were  flat  (constant  high  value,  Js=Ja+Jp/  on  the  inside  and  a 
constant  lower  value,  Js=Ja-Jp/  in  the  outside),  then  the  integral  averaging  process 
would  give,  in  place  of  Eq.  51, 

(Rac)sol/ (Rac^late  =  [  Ja^  +  Jp^]/  0^^) 

=  1  +  (Jp/Ja)^  Flat  Conductor  Eq.  53 

Using  this  uniform  field  (worst  case)  correction  factor,  we  would  have  the  correction 
factor  for  the  ac  resistance  as  (Rac)sol/ (Rac)coax  =  1.0134  for  N=1  for  the  ribbon  loop 
and  (Rac)sol/(Rac)coax  =  1*154  for  the  N=7  turn  ribbon  solenoid.  This  is  still  a 
reasonably  small  correction  factor,  but  it  should  not  be  omitted  from  wirebond 
solenoid  Q  calculations. 


Skin-Effect  Surface  Impedance  of  Thick  Metal  Conductors 

The  dominant  effect  determining  the  high-frequency  losses  or  Q  in  inductors 
is  the  "skin  effect",  or  reduction  of  the  depth  of  penetration  of  ac  fields  into 
conductors  as  frequency  increases.  The  depth  of  penetration,  or  skin  depth,  5,  is 
given  in  any  E&M  text  (e.g.,  Ramo,  Whinnery  and  Van  Duzer,  "Fields  and  Waves  in 
Communications  Electronics",  sections  4.12  and  5.14)  for  a  frequency,  f,  as 

5  =  [p/(7i|if)]*^-^  Eq*  54 

where  for  most  metals  of  interest  as  conductors  (non-ferromagnetic,  iJ.r=l/  metals) 
the  magnetic  permeability,  p=po  Pn  is  just  p=p  o=47t  x  10-7  henries/ meter.  The  values 

for  20°C  metal  resistivity,  p,  in  Eq.  54  are  taken  here  as  p=1.7241xl0-8  Q-m  for 
copper,  p=2.440xl0-8  fl-m  for  gold  and  p=2.620xl0-8  il-m  for  aluminum.  The  skin 

depth,  5,  as  given  by  Eq.  54,  is  defined  by  the  exponential  attenuation  of  the  electric 
field  strength  of  a  normal-incidence  plane  wave  onto  a  flat  plate  of  the  conductor,  of 
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thickriBss  large  in  comparison  to  5,  where  the  magnitude  of  the  field  strength  at  a 

depth  X  into  the  conductor  decays  as  e"^/^.  This  leads  to  a  surface  impedance 
(internal  impedance  for  a  unit  length  and  width  of  the  conductor),  Zs,  for  the 

angular  frequency  co=27cf,  given  as 

Zs  =  Rs  j®Ri 

where  the  real  part  (the  surface  resistance,  Rs,  from  which  Rac  is  obtained)  is 


Rs  =  p/6  =  [7tf|ip]0-5  Eq.  56 

and  the  imaginary  part  (internal  inductance  in  surface  region  of  the  conductor)  is 

Li  =p/(co5)  =  Rs/co  Eq.  57 

Figures  26,  27  and  28  show  convenient  plots  of  the  skin  depth,  5,  and  surface 
resistance,  Rs  versus  frequency  for  copper,  gold  and  aluminum  metal  conductors 
respectively. 

As  an  example  of  the  use  of  Rs  to  calculate  the  ac  resistance  of  a  conductor, 
consider  the  case  of  a  thick  (metal  thickness,  tm  »  25)  straight  metal  plate 
conductor  of  width,  w  (an  example  might  be  a  stripline  signal  conductor.  While  the 
dc  resistance  of  this  plate,  from  Eq.  47  is  given  as 

R<ic/l  =  p/Ac  =  p/(lmw)  Eq.  58 

the  ac  resistance,  Rac/  is  obtained,  assuming  uniform  field  around  the  periphery, 
from  the  surface  resistance,  Rs,  and  the  periphery,  2w  (assuming  here  that  tm  «  w, 
so  that  the  actual  periphery  [=2w+2tm]  is  approximately  equal  to  2w),  is  given  by 

Rac/1  =  Rs/ [periphery]  =  Rs/(2w)  Eq.  59 

In  addition  to  the  requirements  for  the  validity  of  Eq.  59  that  the  conductor  be  much 

thicker  than  25  and  that  the  magnetic  field  be  reasonably  uniform  around  the 
periphery  of  the  conductor,  there  is  an  additional  requirement  that  the  radius  of 

curvature  of  the  surface  of  the  conductor  be  much  larger  than  the  skin  depth,  5.  For 
example,  for  a  circular  wire  with  radius,  a,  much  greater  than  5,  Rac  is  given  by 

Rac/1  =  Rs/(27ia)  Circular  Wire  for  a  »  6  Eq.  60 

When  a  «  5,  Eq.  48  can  be  used.  In  between,  things  are  much  more  complicated. 
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Calculation  of  the  ac  Resistance,  Rac,  for  Circular  Cross-Section  Wire 

The  low-frequency  (a  «  5,  Eq.  48)  and  high-frequency  (a  >>  5,  Eq.  60) 
expressions  for  the  resistance  of  circular  cross-section  wires  are  very  simple,  but  in 
between  these  asymptotic  regions  the  expressions  for  Rac  are  rather  messy  (see 
Ramo,  Whinnery  and  Van  Duzer  sections  5.16, 5.17  and  5.18).  Defining  q  =  (aV2)/  5, 
we  have  for  Rac/1  (in  ohms/meter)  for  a  circular  wire  of  radius  a, 

Rac/1  =  [Rs/(7taV2)][Ber(q)  Bei’(q)  -  Bei(q)Ber'(q)]/[(Ber’(q))2  +  (Bei'(q))2]  Eq.  61 
while  the  internal  inductance  of  the  circular  wire  is  given  by 

coLi/1  =  [Rs/(7caV2)][Ber(q)  Ber’(q)  -  Bei(q)Bei’(q)]/[(Ber'(q))2  +  (Bei'(q))2]  Eq.  62 

In  these  expressions,  Ber(q)  and  Bei(q)  refer  to  the  real  and  imaginary  parts  of  the  Jo 
Bessel  function  of  a  complex  argument  (referred  to  as  Kelvin  functions)  which 
fortunately  are  available  in  tables  and  are  defined  functions  in  advanced  math 
packages  such  as  "HiQ"  (and  presumably  Mathematica).  Ber'(q)  and  Bei'(q)  refer  to 
the  first  derivatives  of  these  functions,  which  are  not  tabulated  or  defined  functions 
in  math  packages,  but  which  may  be  obtained  by  performing  numerical 
differentiation  operations  on  the  Ber  and  Bei  functions. 

The  "HiQ"  script  used  to  calculate  the  ac  resistance  per  meter,  Rac/1,  and  the 
internal  inductance,  Li/1  (referred  to  as  Lint  or  Lint/1  in  the  calculations  and  plots)  is 
shown  in  Figure  29.  This  code  was  originally  developed  as  part  of  a  calculation  of 
the  transient  response  of  small  (hence  lossy)  coaxial  lines,  in  which  the  transient 
response  is  obtained  by  taking  an  8192-point  IFFT  (inverse  Fourier  transform)  of  the 
complex  frequency-domain  response.  Because  of  the  substantial  number  of 
frequency  points  involved,  the  code  was  streamlined  to  approximate  the  expressions 
in  Eqs.  61  and  62  where  possible  (at  higher  frequencies)  with  analytical 
approximations  I  developed.  These  are  far  more  accurate  than  Eq.  60  (which  gives 
6%  to  7%  errors  for  1.5  mil  bond  wires  even  at  f=lGHz),  and  can  avoid  the  need  for 
getting  involved  in  the  messy  Ber/Bei  calculations  in  some  frequency  ranges. 

Figures  30a  and  30b  show  the  results  of  calculations  of  Rac/1  (in  Q/m)  and 
Li/1  (plotted  in  nH/m  so  that  it  can  be  plotted  on  the  same  scale  as  Rs)  for  1.5  mil 
(0.0015"  diameter,  or  a=0.00001905  meter  radius)  straight  circular  wires  made  of 
aluminum,  gold  and  copper.  Figures  31a  and  31b  show  the  same  for  1.0  mil  diameter 
(a=0.00001270  meter  radius)  wires.  As  expected  from  Eq.  48,  a  comparison  of  Figures 
31a  and  31b  with  30a  and  30b  shows  that  the  calculated  Rdc/1  is  2.25x  larger  for  the 
1.0  mil  diameter  wires  as  compared  to  1.5  mil  wires  of  the  same  material,  while  at 
200  MHz,  Rac/1  for  1.0  mil  wires  is  about  1.58x  larger  than  for  1.5  mil,  or  at  1  GHz, 
about  xl.54  larger  (as  compares  to  the  factor  of  1.50  suggested  from  Eq.  60,  with  the 
absolute  values  of  Rac  from  Eq.  60  in  error  by  a  substantial  factor). 
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Calculation  of  the  Q  of  a  Circular  Solenoid  (Circular  Wire  Coil) 

Since  we  have  the  capability  for  calculating  both  the  inductance,  L,  (through 
the  stacked-loop  approach,  as  used  for  Figures  10  -  15  and  18  and  Table  1,  as 
executed  in  the  EXCEL  spreadsheet  calculation,  "Any  S/R  Solenoid  L  Calc  vS")  and 
the  ac  resistance,  Rac/  of  the  conductors  (Figures  29-31),  we  can  calculate  the  Q  of  a 
wire  solenoid  inductor  from  Eq.  39.  EXCEL  spreadsheet  calculations  (filenames:  "1.5 
mil  A1  Solenoid  L,  Q(f)  V2"  and  "1.0  mil  A1  Solenoid  L,  Q(f)  VI"  [the  latest  versions  of 
the  calculation,  using  lookup  table  interpolation],  and  1.5  mil  Au  Solenoid  L,  Q(f) 
V4"  and  "1.0  mil  Au  Solenoid  L,  Q(f)  V3"  [earlier  versions])  were  coded  which 
calculate  values  of  inductance,  L  (external  inductance,  Lext  well  as  total 
inductance,  Lt=Lext+Lint(f)/  at  frequencies  of  f=25MHz  and  200Mh2  are  shown  in 
the  spreadsheet;  Lt  is  calculated  for  all  frequencies  for  the  Q(f)  calculations)  and  the 
Q(f)  values  (at  100  frequency  points  from  f=25MHz  to  2.5GHz)  for  numbers  of  turns 
in  the  circular  solenoid  from  N=1  to  N=24  turns  (easily  extendible  to  greater 
numbers  of  turns).  Because  EXCEL  doesn't  support  the  exotic  Kelvin  and  elliptic 
integral  functions,  lookup  tables  are  required  in  the  spreadsheets  for  M(d/R)  for  r=R 
[Eqs.  21  and  22]  and  for  Rac(f)  and  Lint(f)  [Eqs.  61  and  62],  with  the  values  for  these 
tables  originally  generated  from  "HiQ"  calculations  (Listing  2  and  Figure  29).  The 
M(d/R)  for  r=R  lookup  table,  "R=r=l  Mutual  Inductance  vs  d",  is  the  same  for  and 
wire  diameter  or  metal  composition  as,  unlike  the  self-inductance,  the  mutual 
inductance  between  loops  is  not  a  fimction  of  wire  diameter.  The  spreadsheet  "1.5 
mil  A1  Solenoid  L,  Q(f)  V2"  is  the  latest  version  of  the  calculation,  using  interpolation 
between  lookup  table  points  instead  of  the  standard  EXCEL  lookup  fimction  for 
better  accuracy.  Rac(f)  hint(f)  are  functions  of  both  wire  diameter  and  skin 

depth,  5(f),  and  hence  these  lookup  tables  (filenames:  "Rac/1  for  l.Omil  Al,Au&Cu 
Wires"  and  "Rac/1  for  1.5mil  A1  &  Gold  Wire")  are  specific  to  the  particular  wire 
diameter  and  metal  resistivity  chosen  for  the  solenoid  inductor  design. 

Because  of  the  size  of  these  lookup  tables,  and  the  fact  that  the  100  x  24 
calculated  Q(f,  N)  data  array  is  already  substantial  in  size,  in  order  to  keep  the 
spreadsheet  reasonably  small,  the  stacked-loop  calculation  of  the  external 
inductance,  Lext/  of  the  solenoid  is  streamlined  in  comparison  to  the  full  square  Mij 
array  used  in  the  "Any  S/R  Solenoid  L  Calc  v3"  spreadsheet.  The  approach  is,  given 
the  loop  radius,  R,  and  the  turns  pitch,  S,  to  look  up  from  the  "R=r=l  Mutual 
Inductance  vs  d"  lookup  table  the  mutual  inductance  values  corresponding  to  the 
various  separations  between  turns  (i.e.,  Mi2  for  d/R=S/R  [nearest  neighbors],  M13 
for  d/R=2S/R  [2nd  nearest  neighbors],  M14  for  d/R=3S/R,  M15  for  d/R=4S/R, 
etc.).  For  an  N-tum  solenoid  with  constant  loop  spacing,  S/R,  there  are  just  N-1  of 
these  discrete  mutual  inductance  values  involves  in  the  calculation  (e.g.,  from  Eq.  28, 
all  nearest-neighbor  mutual  inductances  are  the  same  [Mi2=M21=M23=M32=M34= 
M43=M45=M54],  all  2nd-nearest  neighbor  M's  are  the  same,  etc.).  In  addition  to 
these  N-1  mutual  inductance  terms,  scaled  from  the  "R=r=l  Mutual  Inductance  vs  d 
lookup  table  by  multiplying  the  R=1  meter  table  values  by  the  actual  loop  radius  (in 
meters),  tiie  external  self-inductance  of  the  loop,  Ls  (referred  to  as  Le  in  Eqs.  23,  25 
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and  29,  but  we  will  term  Ls  to  differentiate  it  from  the  external  inductance  of  the 
whole  solenoid,  Lext)  is  calculated  from  a/R  using  Eq.  29  (or  a  more  precise  value 
may  be  pasted  in  from  the  "M(r/R)  for  Spiral  Inductors  v2C  file  for  r=R-a  (where 
the  values  were  calculated  from  Eqs.  21  and  22). 

Given  the  Ls=Lls=L2s=L3s...  self-inductance  value  and  the  N-1  discrete 
mutual  inductance  values,  the  calculation  of  the  solenoid  inductance,  Lext/  is 
straightforward.  For  example,  simply  examining  Eq.  28  and  counting  the  equivalent 
terms  gives  for  this  N=5  case 

L  =  5Ls  +  2[  4Mi2  +  3Ml3  +  2Ml4  +  M15  ]  for  N=5  Eq.  63 

Clearly,  for  the  general  case  of  N  turns, 

L(N)  =  NLs  +  2[  (N-1)Mi2  +  (N-2)Mi3  +  (N-3)Mi4  + ...  +  MiN  ]  Eq-  64 

In  the  spreadsheet  calculations,  since  the  full  sequence  (N=l,  2, 3  ...  24)  of  number  of 
turns  was  desired,  the  recursive  form  of  Eq.  64  was  used, 

L(N)  =  L(N-1)  +  Ls  +  2[  Mi2  +  M13  +  M14  +  ...  MiN  ]  Eq.  65 

In  the  Q(f)  calculations,  the  total  inductance  is  used,  which  is  just  the  external 
inductance  from  Eq.  65  plus  the  length  of  the  wire  (l=27tRN)  times  the  Li/1  internal 
inductance  per  meter  (Eq.  62;  taken  from  the  wire-dependent  lookup  table  in  the 
spreadsheet),  with  the  Rac  value  obtained  from  the  Rac/1  column  of  the  same 
lookup  table  in  the  same  way. 

Figure  32  shows  a  plot  of  the  Q(f)  for  frequencies  up  to  f=2GHz  for  a  circular 
cylindrical  coil,  about  84  mils  in  diameter,  of  1.5  mil  diameter  aluminum  wire 
wound  with  a  5  mil  pitch,  for  N=1  to  10  turns  in  the  solenoid  coil  inductor.  (The 
choice  of  the  wire  size  and  wire  pitch  are  for  compatibility  with  the  wire  sizes  and 
pitches  achievable  with  MCM  wirebonding,  and  the  diameter  gives  the  same  area  as 
demonstrated  for  wirebond  loops.)  At  low  frequency  (up  to  50  MHz  or  so),  Rac  is 
nearly  constant  (see  Figure  30a),  and  Lext  (which  is  the  major  portion  of  L)  is  always 
independent  of  frequency,  so  Q=2KfL/R  ac  is  proportional  to  frequency  up  to  50Mhz 
or  so.  At  high  frequencies,  Rac  (Eq.  60),  because  of  the  Vf  dependence  of  Rs  (Eq.  56), 
varies  approximately  as  Vf,  which  leaves  Q  varying  more  or  less  as  Vf,  as  is  clear  in 
Figure  32.  It  is  clear  in  Figure  32  that  the  Q  increases  as  the  number  of  turns  in  the 
solenoid  increases.  If  there  were  no  mutual  inductance  between  turns,  the  Q  would 
be  independent  of  the  number  of  turns,  N,  and  be  equal  to  the  N=1  curve  in  Figure 
31,  because  with  all  Mij=0  in  Eq.  64  you  would  have  L=NLs/  and  wire  length  is 
proportional  to  N,  so  Rac=27iRN(Rac/l)/  which  gives  Q=27tfL/Rac=^fLs/[(Rac/l)E^]/ 
which  is  independent  of  N.  With  the  mutual  inductance  coupling  between  turns 
added,  L  in  Eq.  64  increases,  with  no  increase  in  Rac  (because  the  wire  length  is 
unchanged),  so  the  Q  increases  in  proportion  to  L(N)/NLs, 
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Q(N)  =  Q(N=1)[L(N)/NLs] 


Eq.  66 

where  Ls=L(N=l).  If  the  coupling  between  turns  were  infinitely  tight  (coefficient  of 
coupling,  k=l),  that  is,  Mi2=Ml3=Mi4=Mi5...=Ls,  then  from  Eq.  64  (using  the  fact 
that  the  sum  of  the  first  n  integers=(n/2)(n+l)  with  n=N-l  to  evaluate  the  mutual 
inductance  portion  of  Eq.  64)  we  would  have 

L(N)  =  N^Ls  for  Mi2=Mi3=Mi4...=MiN=Ls  Eq.  67 

For  this  case  of  infinitely  tight  turns  coupling  (completely  unrealistic,  except  with 
the  use  of  high-permeability  cores),  we  would  have  from  Eqs.  66  and  67, 

Q(N)/Q(1)  =  L(N)/NLs  =  N  for  Ml2=Mi3=Mi4...=MiN=Ls  Eq.  68 

Because  of  the  finite  turns  pitch,  S  (or  more  precisely,  the  pitch  to  radius  ratio,  S/R), 
the  mutual  inductance  between  more  distant  turns  in  air-core  solenoids  is  much 
smaller  than  between  nearest-neighbor  turns  (e.g.,  in  the  example  of  Figure  32,  while 
Ls=5.4575nH,  with  S/R=0.116912,  Ml2=3.004nH,  Mi3=2.103nH,  Mi4=1.600nH, 
Mi5=1.264nH,  Mi6=1.201nH,  Mi7=0.694nH,  Mi8=0.581nH,  Mi9=0.490nH  and 
Mio=0.416nH).  Hence  the  increase  of  inductance,  L,  with  increasing  numbers  of 
turns,  N,  from  Eq.  64  falls  far  below  the  N  ^  dependence  of  Eq.  67  as  the  length  to 
radius  ratio  of  the  solenoid,  NS/R,  becomes  large.  Correspondingly,  the  increase  in 
Q  with  increasing  numbers  of  turns,  N,  slows  as  NS/R,  becomes  large,  as  is  quite 
evident  in  Figure  32. 

Along  with  the  "free  Q"  that  the  mutual  inductance  offers,  it  also  has  the 
advantage  of  reduced  frequency  dependence.  As  the  mutual  inductance  contribution 
to  L  increases,  the  relative  contribution  of  the  frequency-dependent  internal 
inductance  of  the  wire  conductor  becomes  relatively  smaller.  For  example,  for  the 
coil  of  Figure  32,  with  a  single  turn  (N=l),  the  external  inductance  is  5.4575nH,  while 
L(200Mhz)=5.6568nH  (+3.65%)  and  L(25MHz)=5.7871nH  (+6.04%  greater  than  Lext)- 
With  N=4  turns,  the  external  inductance  is  9.43x  larger  than  for  N=l,  51.467nH, 
while  L(200Mhz)=52.264nH  (+1.55%)  and  L(25Mh2)=52.785nH  (+2.56%  greater  than 
Lext)-  In  addition  to  this  x2.36  reduction  in  the  change  in  L  with  frequency,  there  is  a 
corresponding  increase  in  Q.  For  example,  at  f=200  MHz,  Q=23.9  for  N=l,  while  for 
N=4  turns,  Q=55.2  at  200MHz,  x2.31  higher.  (This  factor  of  x2.31  for  Q(4)/Q(l)  is 
exactly  that  predicted  from  Eq.  66.) 

Figure  33  (derived  from  the  "1.0  mil  A1  Solenoid  L,  Q(f)  V2"  spreadsheet  calc) 
shows  the  effect  of  using  smaller  diameter  wire  in  the  coil;  1.0  mil  diameter 
aluminum  wire  instead  of  1.5  mil  diameter  as  in  Figure  32.  While  the  turns  pitch  can 
be  reduced  (S=98.2  ^im  is  assumed  in  Figure  33,  as  compared  to  125  pm  in  Figure  32), 
improving  the  mutual  inductances  due  to  the  smaller  S/R,  and  Ls  is  higher  because 
a/R  is  smaller,  the  increased  L  is  offset  by  the  larger  Rac/1  values  so  the  result  is 
lower  Q  (e.g.  for  N=4  turns,  L=58.57nH  and  Q(200MH2)=39.04)  for  1.0  mil  wire. 
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Calculation  of  the  ac  Resistance,  Rac,  for  Substrate  Conductor  Stripes 

In  the  fabrication  of  MCM  wirebond  solenoids,  while  the  top  side  is  created 
from  the  wirebond  loops,  the  bottom  side  of  the  solenoid  is  most  conveniently 
implemented  with  stripes  of  metallization  on  the  MCM  substrate.  In  order  to 
calculate  the  Q  for  such  embedded  MCM  wirebond  inductors,  we  must  calculate  the 
ac  resistance,  Rac/l  of  these  metal  stripe  return  paths  forming  the  bottom  half  of  the 
inductors  (Eqs.  61  and  62  are  used  for  the  wirebond  loops  on  the  top  half).  As  was 
discussed  previously,  for  low  frequencies,  6»tin  (skin  depth  »  metal  thickness), 

Eq.  58  is  used,  while  at  high  frequencies,  5«tin,  for  uniform  fields  Eq.  59  is  used. 
Unfortunately,  for  the  frequency  range  of  interest  and  the  metal  thicknesses,  tm, 

typical  for  MCMs,  usually  5  -tm,  so  neither  of  these  simple  approximations  works. 
Fortunately,  while  this  situation  is  more  complicated  than  the  simple  asymptotic 
expressions  (Eqs.  58  and  59),  it  is  not  as  messy  as  for  the  round  wire  case  (Eqs.  61 
and  62).  The  case  of  a  thin  conductor  sheet  (or  layer  on  another  material)  with  plane 
wave  incident  on  fine  side  is  treated  in  E&M  texts  (e.g.,  Ramo,  Whinnery,  and 
Van  Duzer  sections  5.19  and  5.20),  wherein  the  surface  impedance,  R  spl/  of  a  metal 
plate  of  thickness,  tm,  is  found  to  be,  for  single-sided  conduction, 

Rspl  =  f^s  [smh(2tin/5)  +  sin(2tin/5)]/[cosh(2tin/5)  -  cos(2tin/S)]  1-Side  Eq.  69 

where  Rs  (Eq.  56)  is  the  surface  resistance  of  the  bulk  metal.  (Since  EXCEL  does  not 
have  the  sinh(x)=(eX-e-x)  and  cosh(x)=eX+e-x)  functions,  these  exponential 
definitions  of  the  functions  are  used  in  the  spreadsheet  calculations.)  In  fact,  as 
discussed  in  conjunction  with  Figures  24  and  25,  and  Eq.  53,  the  substrate  stripe 
conductors  in  the  wirebond  solenoid  inductors  have  substantial  magnetic  field 
strengths  on  both  sides  of  the  metal  "plate".  If  these  fields  are  equal,  then  the 
conduction  on  both  sides  will  be  equal  and  the  center  line  should  be  a  plane  of 
synunetry,  so  we  should  be  able  to  consider  the  plate  as  two  plates  of  thickness  tin/2 
in  parallel,  each  with  single-sided  conduction  as  m  Eq.  69.  Hence,  for  equal  2-sided 
conduction  of  a  plate  of  thickness,  tna,  we  should  have 

Rsp2  =  (Rs/2)  [sinh(tm/5)  +  sin(tm/5)]/[cosh(tm/5)  -  cos(tm/5)]  2-Sided  Eq.  70 

For  a  metal  stripe  of  width,  w,  then  the  ac  resistance  per  meter  is  given  as 

Rac/l  =  Rsp2/w  for  Equal  2-Sided  Conduction  Eq.  71 

When  the  conduction  (H-field)  on  the  two  sides  of  the  conductor  are  somewhat 
different,  the  Rac  correction  factor  from  Eq.  53  may  be  used  with  Eq.  71  as  a  good 
approximation,  but  the  "resonance  effect",  wherein  Eq.  69  sometimes  gives  an  Rspl 
value  somewhat  lower  than  the  bulk  surface  resistance,  Rs,  wll  presumably  not  be 
accurately  modeled  in  this  asymmetric  conduction  case. 
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The  Rspl  and  Rsp2  thin  plate  sheet  resistance  values  (Eqs.  69  and  70)  are 
plotted  versus  frequency  for  aluminum  in  Figure  34,  as  derived  from  the 
spreadsheet  calculation  "Thin  Plate  R(Freq)  for  A1  V2".  Note  that  the  "resonance 
effect"  never  exceeds  8.29%  (that  is,  the  single-sided  Rspl  value  never  drops  below 
91.71%  of  Rs,  or  Rsp2  is  never  less  than  91.71%  of  Rs/2),  so  that  the  fact  that  this 
effect  may  not  be  modeled  well  by  Eq.  70,  particularly  for  substantially  asymmetric 
current  distributions,  is  not  a  serioxis  problem.  It  is  important,  however,  to  apply  the 
correction  factor  from  Eq.  53  when  the  H-fields  on  the  two  sides  are  substantially 
different.  Of  course,  if  virtually  all  of  the  H-field  is  on  one  side  of  the  conductor,  then 
Rspl  (Eq.  69)  should  be  used  in  Eq.  72  instead  of  Rsp2  (i-e.,  for  the  given  metal 
thickness,  tm,  the  upper  curve  at  right  side  in  Figure  34  should  be  used  instead  of 
the  lower  curve)  to  calculate  Rac- 


Calculation  of  the  Q  of  an  MCM  Wirebond  Solenoid  Inductor 

The  goal  of  all  of  this  is  to  develop  the  capability  to  predict  the  performance 
of  MCM  embedded  inductors.  In  particular,  we  are  interested  in  solenoid  inductors 
fabricated  using  wirebond  loop  conductors,  since  these  offer  much  higher 
performance  and  lower  crosstalk  than  planar  spiral  inductors,  and  there  are  no  CAE 
tools  which  handle  this  structure  (as  there  are  for  spiral  inductors).  Figure  35  shows 
the  general  structure  of  an  MCM  embedded  wirebond  inductor  using  substrate 
metallization  returns.  (An  alternative  structure,  also  investigated  by  nCHIP,  uses 
low-profile  wirebonds  for  the  return  path,  but  iinless  the  substrate  metallization  is 
very  thin,  this  structure  will  have  lower  Q,  3S  well  as  lower  inductance,  than  that  of 
Figure  35.)  The  wire  radius,  a,  substrate  metallization  stripe  width,  w,  length,  Im/ 
and  thickness,  tm,  as  well  as  the  turns  pitch,  S,  symbols  are  in  agreement  with  the 
usage  in  the  various  foregoing  formulas  and  technical  discussions. 

A  calculation  of  the  inductance  and  Q  of  a  wirebond  solenoid  inductor 
structure  requires,  of  course,  information  about  the  detailed  shape  of  the  wirebond 
loops  themselves,  since  it  is  this  shape  which  determines  the  loop  area,  Aloop,  and 
effective  radius,  Reff  (Reff  is  necessary  for  the  calculation  of  the  self-inductance  of 
the  loops,  and,  through  S/Reff,  the  mutual  inductance  terms  as  well),  as  well  as  the 
wire  length  per  turn,  Iwt  (necessary  for  the  calculation  of  the  internal  inductance 
and  ac  resistance  terms.  Lint  and  Rac).  Figure  36  shows  the  e.xtended  elliptical  shape 
assumed  for  the  wirebond  inductors  for  the  L(N)  and  Q(f,N)  calculations  shown  in 
Figures  37  and  38.  The  wirebond  base,  lm=0.330  cm  is  that  taken  for  the  nCHIP  first- 
generation  test  patterns,  and  the  loop  height,  58.5  mils,  is  based  on  an  average  of 
loop  height  measurements  on  1  mil  aluminum  wirebonds  on  a  test  coupon  sent  by 
nCHIP  to  Mayo.  The  calculations  of  Aloop  (=0.03591  cm^)  and  Iwt  (=0.4805  cm)  are 
shown  on  Figure  36. 
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The  determination  of  the  effective  loop  radius,  Reff,  for  this  more-or-less 
semicircular  wirebond  loop  shape  in  Figure  36  is  not  straightforward.  There  are,  in 
fact,  potentially  two  values  needed  for  Reff:  (Reff)Ls  (Reff)Mij*  (Reff)Ls  would 
be  the  value  of  R,  which,  substituted  in  the  circular  loop  expressions  (Eqs.  21-23  or 
the  "M(r/R)  for  Spiral  Inductors  v2C’’  file  [Figure  8]),  would  give  a  value  for  the 
calculated  external  self-inductance  of  the  loops  equal  to  the  actual  self-inductance  (as 
needed  for  use  in  Ls  in  Eqs.  25-28  and  63-68).  Similarly,  the  value  of  (R€ff)Mij  would 
be  the  value  of  R,  which,  taken  with  the  actual  turns  pitch,  S,  gives  a  value  of 
d/R=S/R=S/Reff  which  gives  the  correct  values  for  the  inter-loop  mutual 
inductances,  Mij,  in  Eqs.  21-22  or  the  "R=r=l  Mutual  Inductance  vs  d”  lookup  table 
file  [Figure  9].  (These  Mij  values  are  needed  in  Eqs.  26-28  and  63-68  for  the 
calculation  of  the  solenoid  inductance.)  Precise  calculation  of  these  (Reff)Ls 
(I^ff)Mij  values  for  the  "odd-shaped"  (=half-circle)  wirebond  loop  of  Figure  36 
would  represent  a  fairly  "computationally  challenging  problem.  Fortunately,  by  far 
the  strongest  factor  determining  the  inductance  of  solenoids  is  the  area  of  the  loop, 
Aioop/  not  details  of  the  shape.  Hence,  for  the  purposes  of  this  calculation,  we  have 
taken  for  Refh  the  value  which  gives  the  correct  loop  area,  Aioop;  is,  we  assume 

.  (Reff)Ls  =  (R€ff)Mij  =  Reff  =  Sqrt(Aloop/7t)  Eq*  72 

At  some  point  in  the  future,  we  will  try  to  obtain  precise  expressions  for  (Reff)Ls  ^nd 
(Reff)Mij  for  some  specific  wirebond  loop  geometry,  but  the  use  of  Eq.  72  as  an 
approximation  for  Reff  (=0.10692  cm  in  Figure  36)  should  be  fairly  accurate  (slightly 
overestimating  the  L  and  Q  values  for  the  wirebond  solenoids). 

While  the  nCHIP  first-generation  wirebond  solenoid  inductor  samples  used 
1  mil  diameter  aluminum  wire  on  an  S=98.2  jim  pitch,  as  noted  in  conjunction  with 
Figures  32  and  33,  the  use  of  a  larger,  1.5  iml  diameter,  wire  will  give  improved  Q 
(even  though  there  is  some  sacrifice  necessary  in  increased  turns  pitch,  S).  For  the 
calculations  of  Figures  37  and  38,  S=125  pm  was  assumed  wth  dw  =  2a  =  1.5  mil 
diameter  A1  wirebonds  forming  the  loops.  With  5=0.0125  cm  and  Reff=0.1069  cm, 
the  pitch  to  radius  ratio  is  S/Reff=0.1169,  the  same  as  shown  for  the  coil  in  Figure  32, 
so  the  L(N)  values  are  essentially  the  same  as  shown  there  (Lint  is  slightly  different). 
An  EXCEL  4.0  worksheet,  "1.5  mil  A1  Wirebond  L,  Q(f)  V3",  was  coded  for  the 
calculation  of  L(N)  and  Q(f,  N)  for  these  MCM  embedded  wnrebond  inductors.  The 
external  inductance  portion  is  the  same  as  for  the  coil  ("1.5  mil  A1  Solenoid  L,  Q(f) 
V2")  case  used  for  Figure  32,  only  with  the  ac  resistance,  Rac/  and  internal 
inductance.  Lint/  terms  modified  to  reflect  the  mixed  wire  and  substrate  stripe 
contributions.  For  an  N-tum  wirebond  inductor,  Rac  is  given  as 

Rac  =  N  [  Iwt  (Rac/l)wire  +  Im  Rsp2/ w  ]  Eq-  73 

where,  as  shown  with  Figures  35  and  36,  Iwt  is  the  length  of  the  bond  wire  in  one 
loop  and  (Rac)wire  is  from  Eq.  61  or  Figures  30a  and  30b.  In  Eq.  73,  Im  is  the  metal 
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stripG  length,  and  w  its  width,  and  Rsp2  is  the  2-sided  conduction  sheet  resistance 
from  Eq.  70  or  Figure  34.  The  internal  inductance.  Lint  is  handled  the  same  way, 

Rint  =  N  [  Iwt  (Li/l)wire  +  im  Rspl/  (27tf  w)  ]  Eq.  74 

only  using  Eq.  62  (or  Figure  30a)  for  (Li/l)wire/  and  we  are  using  Li  =  Rac/ ®  to 
obtain  the  internal  inductance  for  the  substrate  metal  stripe. 

Figure  37  shows  the  inductance  at  f=200  MHz  (black  solid  curve  and  data 
points)  and  also  the  external  inductance  (gray  curve)  vs.  number  of  turns,  N,  for  N=1 
to  24  turns  in  the  MCM  wirebond  solenoid  inductor  of  Figures  35  and  36.  Note  that 
N=4  turns  is  adequate  to  achieve  an  inductance  of  52.68  nH.  Figure  38  shows  the 
Q(f)  for  this  wirebond  inductor  structure  for  a  substrate  metal  thickness  of 
tjY\=25  [iTTi  for  numbers  of  turns  from  N=1  to  10.  The  behavior  of  ihe  Q(F,  N)  curves 
in  Figure  38  is  generally  the  same  as  for  the  wire  coil,  as  discussed  in  association 
with  Figure  32  (Eqs.  66-68),  except  that  the  Q  values  are  higher.  For  example,  for 
N=4,  the  wirebond  solenoid  gives  Q(200  MHz)=58.18  or  Q(1  Ghz)=135.8,  whereas 
the  circular  coil  of  the  same  1.5  mil  A1  wire  gave,  for  N=4,  Q(200  MHz)=55.19  or 
Q(1  Ghz)=133.17.  The  higher  Q  values  are  because  the  ac  resistance  of  the  w=100  pm 
substrate  aluminum  stripe  is  lower  than  that  of  the  dw=2a=38.1  pm  diameter  A1 
wire,  more  than  offsetting  the  longer  conductor  path  length  that  the  wirebond 
inductor  has  for  the  same  l^ff.  Note  that  due  to  the  approximations  associated  with 
Eq.  72,  and  the  fact  that  we  have  not  included  the  small  current  non-imiformity 
correction  for  the  substrate  stripes  (Eqs.  52  and  53),  the  values  of  L  and  Q  calculated 
here  will  be  slightly  high,  and  even  more  so  if  there  are  metal  planes  (package  lids, 
surroxmding  substrate  ground  planes,  etc.)  in  proximity  to  the  inductor  (as  these 
reduce  inductance,  and  may  increase  losses  somewhat).  In  addition,  sometimes  real 
metals,  due  to  surface  roughness  and  non-ideality  of  the  deposition  approach,  give 
skin-effect  losses  somewhat  greater  than  calculated  from  the  nominal  resistivity  of 
the  metal. 

With  all  of  these  effects,  it  would  be  wise  to  margin  the  calculated  Q  value  by 
at  least  20%  in  comparison  to  the  requirements  of  the  real  circuit  application.  For 
example,  for  the  customer  with  the  requirement  for  Qmin=‘=0  f=200  MHz  in  a 
47  nH  inductor,  a  calculated  Q=50  or  higher  would  be  a  reasonable  target.  Since  the 
MCM  embedded  wirebond  structure  gives  a  calculated  Q(200  MHz)=55.2  for  N=4 
turns  (L=52.68  nH)  this  should  be  more  than  adequate  to  meet  the  customer's 
requirements.  In  fact,  nbte  that  we  actually  have  substantially  more  Q  than  required 
to  meet  the  customer's  needs.  Further,  of  the  200  MHz  ac  resistance  of  Rac=l•141^1, 
the  majority,  (Rac)wire=0-851f^,  comes  from  the  wire,  and  only  (Rac)stripe=0.290n 
comes  from  the  tm=25  pm  thick  by  w=100  pm  wide  substrate  stripe.  Hence,  there 
would  appear  to  be  ample  opportunity  to  save  money  in  the  MCM  process  by 
reducing  the  substrate  aluminum  metallization  thickness  from  25  pm  to  a  value 
more  compatible  with  standard  MCM  processing. 
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Standard  MCM-D  fab  lines  typically  use  fairly  thin  metal  planes  for  their 
interconnects  (e.g.,  standard  nCHIP  processes  use  aluminum  thicknesses  varying 
from  1  pm  to  3  pm).  A  number  of  these  layers  could  be  used  in  parallel  to  increase 
the  effective  thickness  of  sputtered  aluminum,  but  going  to  a  25  pm  metal  thickness 
would  require  fairly  major  additions  to  the  process  which  would  presumably  add 

expense.  Because  of  the  frequency  dependence  of  the  skin  depth,  5,  the  selection  of 
an  adequate  value  for  metal  thickness,  tm/  wiU  be  strongly  dependent  on  operating 
frequency.  What  we  would  like  is  a  quantitative  evaluation  of  what  metal  layer 
thickness,  tni/  would  be  adequate  for  the  substrate  return  path  stripes  in  wirebond 
solenoid  inductors  to  achieve  a  given  Q  value  at  a  given  frequency.  The  spreadsheet 
calculation  described  above  is  capable  of  providing  this  information.  As  an 
illustration  of  this  capability,  the  "1.5  mil  A1  Wirebond  L,  Q(f)  V3"  spreadsheet  was 
duplicated  using  different  metal  thicknesses  for  the  metal  plate  impedance  lookup 
table  [obtained  from  the  "Thin  Plate  R(f),  L(f)  for  A1 V2"  spreadsheet].  The  results  for 
the  wirebond  solenoid  structure  of  Figures  35  and  36  for  substrate  aluminum  layer 
thicknesses  of  tm=1.5  pm,  3.0  pm,  6  pm,  9  pm,  12  pm,  15  pm  and  25  pm  are 
summarized  for  f=200  MHz  in  the  "200MHz  Q(N)  vs  Stripe  t.Sum"  spreadsheet. 
These  results  are  plotted  as  Q(200MHz)  vs.  N  in  Figure  39.  This  same  data,  only  for 
the  case  of  f=950  MHz,  is  summarized  in  the  "950MHz  Q(N)  vs  Stripe  t.Sum", 
spreadsheet,  and  is  plotted  as  Q(950MHz)  vs.  N  in  Fig.  40. 

We  can  see  from  Figure  39  that,  while  at  200  Mhz,  tra=15  pm  gives  optimum 
Q,  with  tm=12  pm  only  slightly  poorer  (<3%  lower  Q  than  15  pm),  even  thinner 
metals  would  suffice  to  meet  the  Q(200  MHz)=40  for  L=47  nH  customer  specification 
cited  earlier.  Assuming  we  would  need  to  target  a  calculated  Q=50  to  be  confident  of 
an  actual  delivered  Q=40,  for  N=4  turns,  tin=6  pm  (calculated  Q=46.39)  would 
appear  to  be  a  bit  too  thin,  but  tm=9  pm  (calculated  Q=53.2)  should  be  more  than 
enough  metal  thickness  to  meet  the  specified  200  MHz  Q  value  with  this  wirebond 
inductor  geometry. 

Figure  40  shows  the  Q(N)  at  f=950  MHz  for  this  same  set  of  metal  thicknesses. 
Here  we  see  that  a  metal  thickness  of  tin=9  pm  gives  optimum  Q  (even  better  than 
for  tin=15  pm  or  25  pm  due  to  the  "resonance  effect"  in  Eqs.  69  and  70,  discussed  in 
association  with  Figure  34),  but  a  metal  thickness  of  tin=6  pm  gives  calculated  Q 
results  almost  as  good  as  tin=9  pm.  Because  of  the  questionable  accuracy  of 
modeling  of  the  resonance  effect  by  Eq.  70  for  asymmetric  fields,  these  tin=6  pm  Q 
values  should  probably  be  margined  by  5%  or  so  beyond  the  usual  design  margin. 
In  any  event,  at  these  higher  frequencies,  there  is  no  need  to  have  particularly  thick 
substrate  metal  layers. 

While  the  substrate  metallization  thicknesses  do  not,  for  200  MHz  and  higher 
frequency  operation,  need  to  be  very  thick  (e.g.,  25  pm  is  not  required),  there  is  no 
margin  for  reducing  the  diameter  of  the  wirebond  wire  in  the  solenoid  inductors. 
For  example.  Figure  41  shows  the  Q(f)  curves  for  various  values  of  the  number  of 
turns,  N,  from  N=1  to  10,  for  the  case  of  a  1.0  mil  aluminxim  wire  diameter  (all 
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curves  for  tin=25  ^im).  Note  in  Figure  41,  for  example,  that  for  N=4  turns  (L=58.96 
nH),  the  calculated  Q(200  MHz)  value  for  the  1.0  mil  bondwire  is  43.31,  less  than 
75%  of  the  Q(200  MHz)  for  the  1.5  mil  diameter  wire  case  (Figure  38),  and  not  nearly 
enough  to  provide  suitable  design  margin  for  meeting  the  customer’s  specification 
(calculated  Q=50  required  to  insure  meeting  Q=40  specification).  Further,  for  any 
given  Q,  thinner  substrate  metal  can  be  used  with  the  thicker  wirebond  wire,  for  a 
net  cost  savings,  so  going  with  the  thicker  bond  wire  is  best.  (The  use  of  ribbon 
bonds  in  place  of  round  wire  wirebond  loops  could  also  prove  advantageous  in 
some  cases). 


Accounting  for  Capacitance  Effects  on  MCM  Embedded  Inductors 

Along  with  the  desired  inductance  characteristics  of  MCM  wirebond 
inductors,  nature  tmavoidably  brings  undesired  capacitance  effects.  I  have  carried 
out  extensive  calculations  of  the  tum-to-tum  capacitance  for  the  wirebond  solenoid 
embedded  MCM  structure  of  Figure  35.  This  data,  taken  along  with  the  tum-by-to 
mutual  inductance  and  self-inductance  parameters  calculated  in  the  various 
spreadsheets  discussed  previously,  form  the  basis  for  quite  a  detailed  ac  model  for 
the  embedded  inductors  (e.g.,  SPICE  models).  While  these  detailed  models  should 
be  quite  accurate  in  predicting  the  (complex)  ac  impedance,  Z(f),  of  the  embedded 
inductor  structures,  it  is  important  to  begin  with  a  clear  rmderstanding  of  just  what 
these  models  imply  with  regard  to  the  usual  parameters  quoted  for  such  devices. 
The  parameters  most  generally  quoted  are  the  inductance,  L,  the  self-resonant 
frequency,  fo  (or  SRF),  the  quality  factor,  Q(f),  and  the  internal  capacitance,  Cq,  of 
the  inductor. 

It  is  important  to  imderstand  that  when  a  single  number  such  as  the  internal 
capacitance,  Co,  or  self-resonant  frequency,  fo,  is  cited,  it  is  necessarily  within  the 
context  of  a  simplified  conceptual  model  of  the  inductor.  The  actual  frequency- 
dependent  impedance  characteristics  can  be  quite  complex  (particularly  for  many- 
turn  inductors),  containing  much  more  detailed  information  (e.g.,  exhibiting 
multiple  resonances,  etc.)  than  can  meaningfully  be  conveyed  by  only  a  couple  of 
parameters.  Figure  42  shows  the  simplified  model  for  inductors  that  is  usually 
assumed  when  these  key  inductor  parameters  are  defined.  Note  that  this  is  a  two- 
terminal  model  for  the  inductor,  appropriate  when  the  inductor  structure  is  well 
isolated  from  ground  planes  or  surroimding  circuitry,  such  that  the  capacitive  (or 
inductive,  for  that  matter)  coupling  to  other  electrodes  may  be  ignored.  Note  that 
this  is  not  always  the  case,  and  the  more  complex  models  we  will  be  generating  for 
SPICE  or  other  circuit  simulations  can,  at  least  in  principle,  take  these  couplings  into 
account.  However  for  these  cases  with  heavy  external  coupling,  the  simple 
parameters  like  self-resonant  frequency,  internal  capacitance,  etc.  are  in  general  not 
uniquely  defined  (they  tend  to  be  dependent  on  how  the  inductor  structure  is 
connected  in  the  circuit). 
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The  standard  simplified  model  for  the  isolated  inductor  in  Figure  42  is  that  of 
the  inductor,  L,  and  its  ac  series  resistance,  Rac/  shunted  at  the  inductor  terminals 
with  a  lossless  internal  capacitance,  Cq.  As  noted  in  conjunction  with  Eq.  40,  for 
reasonably  high  values  of  Q  (or  correspondingly,  low  values  of  Rac)/  with  the 
inductor  shunted  at  its  terminals  by  an  external  lossless  capacitor  of  capacitance 
value,  Cext/  the  resonant  frequency,  fr,  will  be  given  by 

fr  =  1/{27C  Sqrt[L  (Co  +  Cext)]}  Eq.  75 

or  the  reciprocal  of  the  resonant  frequency  by 

1/fr^  =  4k^  L  (Co  +  Cext)  Eq.  76 


This  straight-line  relationship,  illustrated  in  Figure  42,  forms  the  basis  for  the 
functional  definitions  for  the  basic  inductor  parameters.  (In  fact,  this  is  how  these 
parameters  are  defined/measured  on  unknown  physical  inductors  using  the 
HP  4342A  "Q  Meter").  In  general,  inductors  are  intended  for  use  at  frequencies 
substantially  lower  than  their  self-resonant  frequency  (SRF),  fo,  given  by  Eq.  75  with 
Cext=0  as 


fo  =  l/2WLCo  (SRF)  Eq.77 

In  fact,  the  simplified  inductor  model  of  Figure  42  may  be  thought  of  as  a  second- 
order  approximation  to  the  low-frequency  behavior  of  the  real  physical  inductor 
(where  the  first-order  approximation  would  be  Z(f)=Rac+j2JifL)- 

While  this  simplified  model  works  very  well  for  (reasonably  well  isolated) 
inductors  at  low-to  medium  frequencies,  at  very  high  frequencies  (near,  or  above, 
the  SRF),  all  of  the  complexities  of  multiple  resonances,  etc.,  can  make  the  behavior 
of  real  physical  inductors  more  complex  than  predicted  by  the  simplified  model.  In 
most  cases  this  is  not  a  serious  problem  in  the  application  of  the  inductors,  since  they 
are  not  normally  used  in  this  frequency  regime,  but  it  can  prove  a  problem  in 
characterizing  the  inductors  in  terms  of  the  parameters  of  the  simplified  model.  For 
example,  if  one  were  to  simply  open-circuit  the  inductor  and  attempt  to  measure  the 
SRF  in  order  to  obtain  Co  from  Eq.  77,  one  might  well  end  up  with  parameter  values 
that  rather  poorly  characterize  the  inductor  behavior  in  the  frequency  range  of  its 
intended  use.  What  this  means  in  terms  of  the  1/ vs.  Cext  plot  in  Figure  42  is  that 
while  the  upper  right  portion  of  the  straight  line  (e.g.,  for  Cext^Ci  in  Figure  42)  is 
indeed  straight,  at  much  higher  frequencies,  near  the  Cext=0  there  could  be 
very  substantial  departures  form  this  "straight-line"  relationship.  Experimentally, 
this  problem  is  avoided  by  carrying  out  all  characterization  measurements  with 
lossless  external  capacitance,  Cext  values  (e.g..  Cl  and  Cz  in  Hgure  42)  that  are  large 
enough  so  that  the  measurement  frequencies  are  kept  in  the  well-behaved  range  of 
interest.  As  seen  in  Figure  42,  the  use  of  two  (or  more)  external  capacitance  values 


Page  39 


allows  Co  to  be  obtained  directly  as  the  negative  intercept  of  the  1/ i-p-  vs.  Cext  plot 
(and  the  true  inductance,  L,  value  to  be  obtained  using  the  formula  shown  in  the 
figure),  with  the  SRF  then  calculated  from  Eq.  77. 

Note  that  these  techniques,  intended  for  evaluating  the  parameters  of 
unknown  physical  inductors,  work  equally  well  for  the  purpose  of  extracting  these 
basic  simplified  model  parameters  (e.g.,  L,  Cq,  Q(f)  and  SRF)  from  much  more 
detailed  composite  SPICE  models  for  the  devices.  A  very  common  problem  is  that  a 
very  detailed  device  model  gets  generated  (e.g.,  containing  all  of  the  distributed 
capacitance,  mutual  inductance,  self  inductance,  etc.  information  on  a  turn-to-tum 
basis),  but  one  is  hard  pressed  to  imderstand  from  the  elegantly  detailed  model  just 
what  the  L,  Co,  SRF  or  Q(f)  values  for  the  inductor  are.  One  simply  runs  the  inductor 
model  two  or  more  times  in  SPICE  over  frequency  with  fixed  capacitor  values  (e.g.. 
Cl,  C2...  in  Figure  42)  shimting  the  terminals,  using  a  unit  ac  current  source  drive  to 
find  the  resonant  frequencies,  fi,  f2,  etc.  The  values  of  Cl,  C2,  etc.  should  be  selected 
to  keep  the  resonant  frequency  values,  fi,  f2— /  (and  3dB  bandwidth  values  for  Q(f)), 
within  the  range  of  the  intended  application  frequency  for  the  device.  The  relations 
described  in  Figure  42  and  Eqs.  75  -  77  are  then  used  to  extract  the  parameters  for  the 
simplified  inductor  model.  It  is  also  useful  as  a  comparison,  to  run  the  simplified 
inductor  model  itself  to  compare  the  predicted  Z(f)  with  that  simulated  from  the 
more  detailed  model.  Quite  frequently  it  will  be  foimd  that  for  well-isolated 
inductors,  the  simplified  inductor  model  will,  with  properly  derived  parameters, 
give  adequate  accuracy  for  Z(f)  over  the  frequency  range  of  interest,  such  that  the 
use  of  more  computationally  intensive  and  "messy"  detailed  models  can  be  avoided 
in  practical  circuit  simulations. 

At  this  point  it  should  be  pointed  out  that  these  definitions  of  inductor 
parameters  and  the  simplified  inductor  model  into  which  they  fit  (and  in  terms  of 
which  they  are  defined)  represent  a  soimd  second-order  model  for  the  frequency 
dependence  of  inductors  in  the  low-to-medium  frequenc}^  regime.  Within  the 
framework  of  this  simplified  inductor  model  in  Figure  42,  the  inductor  Q  is  given  by 

Q(f)  =  2KfL/Rac  Eq- 78 

Note  that  this  definition  of  the  Q  of  an  inductor  having  significant  internal 
capacitance  is  the  same  as  that  of  Eq.  38,  wherein  no  mention  was  made  of  internal 
capacitance  effects,  etc.  The  reason  for  fius  is  that  the  basic  definitions  of  Q  as  noted 
in  Eqs.  33  -  46  are  all  related  to  energy  storage/energy  dissipation  or  resonant 
bandwidth  factors  that  are  uninfluenced  by  the  presence  of  lossless  shimting 
capacitors. 

The  fact  that  the  definition  of  Q  should  not  depend  on  the  presence  of  lossless 
capacitance  elements  shunting  the  inductor  is  clear  from  a  simple  circuit  example. 
Consider  a  comparison  between  two  inductors  having  the  same  true  inductance,  L, 
and  series  resistance  values,  Rac/  but  with  one  having  no  internal  shunt  capacitance 
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(Co=0)  and  one  having  a  finite  lossless  capacitance,  Co=Ciin  parallel  with  it.  If  the 
intended  application  of  the  inductor  is  to  function  in  a  parallel-resonant  (  tank  ) 
circuit  at  a  resonant  frequency  fr=f2  h\  Figure  42,  then  to  use  the  first  (Co=0) 
inductor  in  the  circuit,  a  capacitor  of  value  Cext=C2  is  placed  in  parallel  with  the 
terminals  of  the  inductor.  To  use  the  second  (Co=Cl)  inductor  for  this  same  tank 
circuit  function  at  the  same  frequency,  it  would  be  placed  in  parallel  with  a  smaller 
capacitor,  of  value  Cext=Cl-C2.  In  either  case,  the  total  capacitance  (combined 
internal  +  external)  across  the  inductor  will  be  the  same  (Cl),  and  hence,  assuming 
all  the  capacitors  to  be  essentially  lossless,  the  resulting  circuit  bandwidth,  and 
hence  (by  Eq.  43)  the  Q,  must  be  the  same,  just  as  indicated  by  Eqs.  39  and  78. 

While  the  internal  capacitance,  Co,  of  an  inductor  does  not  directly  lead  to 
problems  with  Q  in  many  applications,  the  value  of  the  internal  capacitance  or  SRF 
can  dictate  whether  a  given  inductor  is  usable,  due  to  circuit  constraints,  in  a  given 
application.  An  obvious  case  would  be  the  tank  circuit  example  above  if  the 
intended  application  required  a  resonant  frequency  higher  than  fr=fl  (the  self¬ 
resonant  frequency  of  the  inductor  having  Co=Cl  in  the  example;  see  Figure  42). 
While  most  bandpass  filter,  and  many  other,  applications  of  inductors  use  them  in 
parallel  with  capacitors  of  substantial  value,  this  is  not  always  the  case,  and  even 
when  it  is,  the  inductor  would  be  usable  only  when  the  value  of  circuit  capacitance 
required  exceeds  Go-  Further,  for  variable-frequency  tuned  circuit  applications 
where  the  inductor  is  operated  in  parallel  with  a  variable-capacitance  element  such 
as  a  varactor  diode,  the  presence  of  self-capacitance  will  reduce  the  tuning  range 
available  for  a  given  varactor  capacitance  ratio.  Hence,  in  general,  it  is  desirable  to 
have  low  Cq  (or  correspondingly,  high  SRF)  for  circuit  reasons  having  little  to  do 

with  Q. 
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Definition  of  Geometry  for  Circular  Loop  Magnetic  Field  Calculation 


so  B^Field  is_given 
as  B(r,d)=nH(r,d) 


Z-Component  of  Magnetic  Field  vs.  Radial  Position 

at  Height,  d,  Above  Circular  Loop  of  Radius,  R,  Carrying  1=1  Ampere 
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Z-Component  of  Magnetic  Field  vs.  Radial  Position 

at  Height,  d,  Above  Circular  Loop  of  Radius,  R,  Carrying  1=1  Ampere 


Z-Component  of  Magnetic  Field  vs.  Radial  Position 

at  Height,  d,  Above  Circular  Loop  of  Radius,  R,  Carrying  1=1  Ampere 
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//  Circular  Loop  Inductor  Mutual  Inductance  Calc  R.  C.  Eden  July  7,  1996 
II  General  calculation  set  up  with  d=0  for  use  with  spiral  inductors 


//  Physical  Constants:  Units  are  meters,  seconds,  farads,  Henrys,  Hz,  etc. 
pi=3. 14159265353979323846; 

mu=(4E-7)*pi;  /*  Permeability  of  Free  Space  in  Henries/meter  */ 

clt=2.99792458E8;  /*  Speed  of  Light  in  m/s  */ 

Eo=l/ (mu*clt*clt) ;  /*  Permitivity  of  free  space  */ 


Calculation  for  Spiral  Inductors;  d=0,  vary  radius  b 
Radius  of  1st  loop  in  meters  */ 

Height  of  2nd  loop  above  1st  loop  in  meters  */ 

Step  size  on  b  in  meters  */ 

Initial  Value  of  radius  b  in  meters  */ 

for  i=l  to  9  do 
b=0.0+i*bstep; 


II  Mutual 

Inductance 

a=l; 

/* 

d=0.0; 

/* 

bstep=0.1; 

/* 

b=0; 

/* 

ksq={4*a*b)  /  (d''2+(  (a+b)  ''2) )  ; 

k=Sqrt  (ksq) ;  /*In  HiQ  argument  of  fn  is  called  k  instead  of  k''2  as  in  text: 
Kfn  =  Cornell (ksq);  /*  Conplete  Elliptic  Integral  of  the  First  Kind  */ 

Efn  =  comel2 (ksa) ;  /*  Complete  Elliptic  Integral  of  the  Second  Kind  ’ 

M  =  mu*Sqrt (a*bf* ( ( { (2/k) -k) *Kfn) - { (2/k) *Efn) ) ; 

//  Barr[i]=b; 

//  Marr[i)=M; 

BMarr[i,l]=b; 

BMarr [ i , 2 ] =M; 


end  for; 

//*  Note  that  HiQ  uses  Abroraiwitz  and  Stegtin  definition  cf  the  K{k)  &  E(k) 

//  Complete  Elliptic  Integral  functions  in  which  k''2  is  the  modulus 

//or  argument  of  the  function,  instead  of  k,as  in  the  older  Jahnke  and  Emde 

//  "Tables  of  Functions"  forms  for  K(k)  and  E(k),  (as  used  in  Ramo,  Whinnery, 

//  and  Van  Duzer  and  most  other  E&M  texts) ,  so  in  HiQ  we  must  square 

//  the  argument  (modulus)  of  the  function  first,  before  calling  the  function: 

//  This  HiQ  Script  file  Calculates  the  Mutual  Inductance  Between  Two  Circulai 
//  Loops  which  are  Co-Axial,  but  may  be  displaced  by  a  height,  d,  between  ths 
//  The  Radii  of  the  two  loops  are  "a"  and  "b"  (dimensions  in  meters) 

//  For  the  case  of  a  spiral  inductor,  the  loops  are  coplanar,  so  d=0 


Mutual  Inductance  Between  Coaxial  Circular  Loops  vs.  Radius  and  Height 

cl/R=0.0 
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Mutual  Inductance  Between  Coaxial,  Coplanar  (d=0),  Circular  Loops  vs.  Radius 
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r/R,  Radius  of  Smalier  Loop,  r,  Normalized  to  Radius  of  Outer  Loop,  R 


Mutual  Inductance  vs.  Separation,  d,  Between  Two  Coaxial  R=r=1  Meter  Circular  Loops 
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//  Inductance  of  Thin-Walled  Cylindrical  Solenoid  R.  C.  Eden  July  8,  1996 
//  Using  Complete  Elliptic  Integral  formula  from  Jahnke  and  Emde  page  87. 

//  L=N^2*  (d)  (1/3)  *  [  { (K{k)  +  { (d/l)  ''2-1)  *E(k) )  /sin (alpha) }-  (d/1)  ^2)  ] 

//  Where  tan (alpha) =d/l=2R/l,  1  is  length,  and  k=sin(alpha)  is  modulus 


//  Physical  Constants:  Units  are  meters,  seconds,  farads,  Henrys,  Hz,  etc. 
pi=3.14159265358979323846; 

mu=(4E-7) *pi;  /*  Permeability  of  Free  Space  in  Henries/meter  */ 

clt=2.99792458E8;  /*  Speed  of  Light  in  m/s  */ 

Eo=l/ (mu*clt*clt) ;  /*  Permitivity  of  free  space  */ 


//Mutual  Inductance  Calculation 

Rr:l;  /*  Radius  of  Solenoid  in  m.eters  */ 

d=2*R;  /*  Diameter  of  Solenoid  in  meters  */ 

SoverR=0.1;  /*  Ratio  of  turns  Spacing,  S,  to  Radius,  R  */ 

Nstep=l;  /*  Step  size  on  Number  of  Turns,  N  */ 

N=l;  /*  Starting  Number  of  Turns  */ 


for  i=l  to  50  do 

N=i*Nstep;  /*  Number  of  Turns,  N  */ 

l=N*SoverR;  /*  Length  of  Solenoid,  1  (meters)  */ 
alpha=arctan(d/l) ; 
k= sin (alpha) ; 

ksq=k*k;  /*  Note  that  HiQ  uses  Abromiwitz  and  Stegun  definition  of  the  */ 
//  Corrplete  Elliptic  Integral  functions  in  which  k''2  is  the  m.oduli 
//  or  argument  of  the  function,  instead  of  k,  so  must  square  k  fii 
tanSq=  (d/1)  ''2; 

Kfn  =  Cornell  (ksq) ;  /*  Complete  Elliptic  Integral  of  the  First  Kind  ’ 

Efn  =  comel2 (ksq) ;  /*  Complete  Elliptic  Integral  of  the  Second  Kind 

L= (N*N*d*mu/3 ) * ( ( {Kfn+ (tanSq-1) *Efn) /k) -tanSq) ;  /*  Inductance  in  Henries  ^ 

NLarr[N, 1)=N;  /*  Store  Number  of  Turns,  N,  in  First  Column  */ 

NLarr [N, 2]=L;  /*  Store  Inductance,  L,  in  Second  Column  */ 

NLarr [N, 3]=1;  /*  Store  Length  of  Solenoid,  1,  in  Third  Column  */ 

end  for; 


//  This  HiQ  Script  file  Calculates  the  Inductance  of  a  Solenoid  of  Radius  R 
//  and  Turns  Spacing  S  (entered  as  S/R  ratio  called  "SoverR"  in  HiQ  Script 
//  for  50  values  of  the  number  of  turns,  N,  (N=l  to  50)  in  the  solenoid 
//  Filename  =  "SoleniodlnductanceCalcJEvlpl" 


h  I 


!/! 


z 


Inductance,  L  [Henries]  for  R=1  meter  Solenoid  having  N  Turns  (Loops)  with  Pitch  S 
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Inductance,  L  [Henries]  for  R=:1  meter  Solenoid  having  N  Turns  (Loops)  with  Pitch  S 
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Inductance,  L  [Henries]  for  R=1  meter  Solenoid  having  N  Turns  (Loops)  with  Pitch  S 


Inductance  vs.  Number  of  Turns  for  Solenoid 
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Inductance,  L  [Henries]  for  R=1  meter  Solenoid  having  N  Turns  (Loops)  with  Pitch  S 
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Inductance,  L  [Henries],  for  R=1  meter  N-Turn  Solenoids  having  Various  Pitches 


Inductance  vs.  Number  of  Turns  for  R=1  meter  Solenoids  with  Various  Pitches,  S/R 
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Inductance,  L  [Henries],  for  R=1  meter  N-Turn  Solenoids  having  Various  Pitches, 


Inductance  vs.  Number  of  Turns  for  R=1  meter  Solenoids  with  Various  Pitches,  S/R 
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Summary  of  R=1  meter  Solenoid  L  vs  N  Data,  Including  Very  Low  Frequency  Internal  Self-Inductance 

Richard  C.  Eden  7/14/96 
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Note  that  S/R  values  of  0.2  or  greater  appear  within  range  of  realizability  of  MCM 
embedded  wirebond  inductor  impiementation  [at  R(effective)=0.0006  meters] 


Inductance,  L  [Henries],  for  R=1  meter  N-Turn  Spirals  having  Various  Pitches, 


CO 


Inductance  vs.  Number  of  Turns  for  R=1  meter  Spirals  with  Various  Pitches,  S/R 


Summary  of  R=1  meter  Spiral  Inductor  L  vs.  N  Data  for  Various  Turn  Pitches,  S/R 
L  Values  Include  Very  Low  Frequency  Internal  Self-Inductance;  Conductor  Radius  a/R=0.4S/R 

Richard  C.  Eden  7/14/96  R  is  Radius  of  Center  of  Conductor  of  Outer  (N=1)  Loop;  Turns  are  Numbered  from  Outside  In 


CD 

in 

in 

in 

in 

in 

in 

in 

in 

m 

UJ 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

CvJ 

TT 

CO 

CM 

CM 

T— 

CM 

1 

o 

OO 

CO 

CO 

CD 

CD 

CO 

r>- 

o 

CO 

s 

00 

CD 

OO 

o 

u 

o 

CO 

CM 

CO 

CO 

CD 

CO 

•M- 

CO 

r- 

CD 

CM 

CM 

CO 

in 

CD 

CD 

CD 

in 

in 

in 

in 

in 

in 

in 

in 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

T“ 

a> 

CO 

CO 

CD 

CO 

CM 

m 

00 

T“ 

o 

tT 

r> 

CO 

CO 

CD 

CO 

in 

• 

O 

o 

TT 

CO 

00 

"O- 

o 

CO 

o 

d 

ir- 

C7> 

in 

CM 

o 

CO 

r- 

CO 

CM 

co 

CD 

8 

o 

CO* 

CM* 

CO* 

rr* 

in 

in 

CD 

CD* 

CD 

in 

in 

in 

in 

in 

in 

in 

m 

in 

UJ 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

CM 

CO 

CM 

CD 

ymm 

CO 

CD 

o 

CD 

o 

o 

CO 

CD 

o 

T" 

N. 

CO 

s 

o 

r-" 

CD 

o 

o 

• 

d 

CO 

r^ 

CM 

o 

in 

o 

T— 

o 

00 

CO 

CD 

K 

CO 

-r- 

CM* 

CO* 

CO 

CD 

in 

in 

in 

in 

in 

lO 

LU 

LU 

LU 

LU 

LU 

LU 

CO 

CD 

CD 

TT 

TT 

CO 

CM 

ymm 

o 

TT 

CO 

O 

in 

O 

CD 

CM 

•o- 

in 

CD 

O 

d 

'T 

CD 

CD 

CD 

CO 

in 

in 

s 

CD 

CM 

CO 

1-^ 

T-^ 

CM 

CO 

CO 

CO 

ID 


CO 

CD 

in 

in 

in 

m 

LU 

LU 

LU 

LU 

LU 

N 

N 

CO 

in 

00 

o 

8 

CD 

CD 

O 

CD 

in 

CO 

« 

d 

00 

CD 

CO 

CM 

O 

CM 

CD 

tT 

CO 

CD 

CM 

CM 

CD 

CO 

in 

in 

in 

§ 

LU 

LU 

LU 

LU 

LU 

CM 

CD 

CM 

CO 

CO 

in 

CM 

in 

s 

o  d 

CD 

T” 

OO 

TT 

CO 

CD 

O 

CD 

55 

CO 

CO 

ir^ 

CD 

CD 

in 

in 

lO 

LU 

LU 

LU 

LU 

in 

CM 

CO 

o 

CM 

Tr“ 

CO 

t— 

o 

o 

CD 

CM 

o 

CO 

y 

in 

in 

y- 

CM 

CM* 

y 

y 

CD 

in 

in 

in 

T 

'C- 

"T 

T 

TT 

T 

XT 

CO 

CO 

CO 

CO 

cq 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

UJ 

LU 

LU 

UJ 

LU 

LU 

LU 

LU 

T“ 

CM 

CD 

CD 

CO 

CD 

in 

00 

CD 

CO 

CO 

00 

in 

m 

'M- 

CO 

CD 

TJ* 

CD 

o 

o 

TT 

in 

CD 

s 

CD 

o 

CO 

CO 

CD 

CO 

r;:: 

o 

CD 

CO 

o 

TT 

rr 

CO 

CO 

CO 

CO 

CO 

CD 

CO 

CO 

T~ 

o 

o 

CM 

in 

O 

00 

d 

* 

CM 

CO 

CM 

in 

CM 

in 

r- 

in 

CO 

CO 

CD 

CM 

f:- 

CO 

OO 

o 

O 

in 

CO 

o 

cq 

o 

1^ 

CM 

o 

CD 

OO 

CO 

CM 

cq 

K 

CM 

in 

CD 

y^ 

y 

CM 

CO 

CO* 

in 

CD 

CO 

00* 

CD 

y 

y 

cq 

in 

m 

in 

T 

? 

T 

'M- 

CM 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

UJ 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

UJ 

CM 

CM 

CD 

00 

CD 

O 

CO 

O 

CM 

in 

in 

CO 

y“ 

T“ 

CO 

TT 

o 

CD 

CD 

o 

O 

CD 

00 

CD 

o 

CM 

CD 

Jo 

o 

in 

CD 

in 

CO 

o 

CD 

TT 

in 

CO 

o 

CM 

CM 

CD 

TT 

CM 

o 

CO 

Jo 

CD 

CD 

00 

CO 

CD 

CD 

CD 

d 

* 

CO 

O 

in 

rr 

O 

o 

■y 

CO 

CO 

CD 

r- 

in 

tT 

CM 

o 

CO 

in 

o 

CM 

in 

in 

in 

CD 

O 

CD 

in 

o 

CO 

CM 

00 

o 

in 

y 

CD* 

CM 

tT* 

N-* 

y^ 

y 

CM 

CM* 

CO 

CO 

in 

in 

CD 

CO* 

r^* 

00* 

cd 

CD 

CD 

CD 

in 

in 

in 

in 

'O- 

'rr 

'T 

'M* 

'M' 

Tf 

T 

TT 

T 

"T 

T 

CO 

rvi 

LU 

LU 

UJ 

LU 

LU 

UJ 

LU 

UJ 

tu 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

UJ 

LU 

LU 

CD 

CM 

CD 

CO 

CM 

CM 

CD 

OO 

CO 

in 

O 

O 

00 

y 

o 

CO 

xr 

o 

o 

CD 

CD 

in 

CD 

CD 

CD 

CD 

O 

CD 

TT 

CM 

• 

o 

CD 

o 

CD 

5> 

CO 

TT 

CO 

CM 

CO 

CO 

CD 

h- 

'M' 

OO 

O 

O 

CD 

o 

■  * 

CO 

o 

CO 

CO 

o 

CO 

r-. 

CD 

in 

CO 

OO 

CO 

CD 

CD 

CO 

CO 

o 

To 

CD 

CD 

CD 

in 

CM 

CD 

CD 

CO 

T— 

in 

CD 

CO 

CD 

o 

cq 

CD 

CD 

CM 

in 

in 

y 

CD* 

CD 

1-^ 

CM 

CM 

CM 

cd 

cd 

cd 

in* 

in 

in 

in 

CD 

CD 

CO 

in 

in 

m 

in 

*T 

xt 

■M- 

T 

TT 

•M- 

"O" 

CD 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

UJ 

LU 

LU 

LU 

LU 

LU 

UJ 

O 

CM 

CD 

o 

CM 

CD 

CM 

O 

t— 

O 

CD 

'M* 

CD 

00 

CO 

O 

in 

CO 

CO 

o 

o 

00 

CO 

r- 

in 

CD 

CO 

CD 

'M' 

O 

r- 

y“ 

O 

CD 

y- 

in 

CD 

o 

o 

TT 

CD 

CO 

m 

in 

CO 

CO 

CD 

CO 

CO 

o 

CM 

O 

CO 

d 

TT 

CM 

CM 

CD 

CM 

o 

O 

o 

00 

CD 

CO 

CO 

CM 

CO 

CO 

O 

in 

o 

CM 

cq 

CD 

CD 

'M- 

CD 

CM 

in 

CO 

o 

CO 

in 

CD 

p 

y~^ 

CM 

p 

in 

1-^ 

cd 

in 

cd 

y 

1-^ 

1-^ 

y 

CM* 

CM 

CM 

cd 

cd 

cd 

cd 

cd 

XT 

’’sr* 

CO 

in 

in 

in 

in 

in 

"M* 

"T 

XT 

T 

•M- 

-O' 

-M- 

■M- 

TT 

T 

LO 

lU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

CM 

CO 

CD 

CM 

r>. 

TT 

CD 

CO 

CM 

CO 

CD 

TT 

CD 

in 

CM 

CM 

▼“ 

O 

"y 

o 

O 

CO 

CD 

CO 

00 

CD 

CD 

CD 

CM 

CM 

O 

O 

CD 

CD 

in 

O 

CD 

CO 

CD 

CO 

CO 

CD 

CO 

O 

O 

in 

CO 

CD 

CO 

•y 

'M- 

y” 

yy 

CD 

Q 

d 

in 

O 

CO 

CO 

CD 

CM 

CM 

o 

CO 

CD 

CD 

CO 

CM 

in 

CO 

o 

•M* 

CD 

o 

p 

CO 

CO 

in 

G> 

CM 

p 

cq 

p 

P 

'M' 

iq 

p 

CO 

p 

p 

p 

y 

CO* 

in 

r--* 

CD 

1-^ 

y 

y 

y 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

<q 

m 

iq 

in 

in 

in 

T 

t 

*T 

*T 

■M" 

1 

CO 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

1 

CD 

CD 

•M- 

in 

CO 

CD 

CM 

in 

CO 

y 

r- 

in 

CO 

CO 

o 

CM 

in 

CO 

CD 

o 

tT 

CM 

CD 

CO 

CO 

in 

y“ 

CD 

CO 

in 

• 

O 

CM 

in 

CD 

r^ 

in 

CO 

O 

CO 

OO 

CD 

CM 

OO 

CO 

CD 

o 

* 

O 

Tj- 

CO 

CD 

CO 

CO 

CM 

o 

y^ 

CD 

xj- 

OO 

CD 

O 

p 

y 

p 

CO 

CO 

p 

p 

iq 

r-. 

p 

p 

p 

p 

p 

cd 

CD 

cd 

y 

y 

y 

y 

t— 

T-’ 

y 

r 


p 

p 

p 

CM 

UJ 

Ljj 

LU 

CO 

TT- 

"M- 

T“ 

CD 

CM 

d 

d 

CO 

CD 

O 

▼— 

in 

CO 

cvi  CD  CO 


in 

p 

LU 

p 

LU 

p 

LU 

•M- 

CO 

CO 

T— 

O 

1^ 

CD 

9 

TT 

in 

o 

o 

CO 

CM* 

in 

in 

o 

p 

CD* 

CO 

p 

LU 

p 

LU 

CO 

CO 

CO 

d 

d 

CM 

•M- 

CO 

CO 

CO 

CO 

LU 

LU 

in 

CM 

CO 

CD 

CM 

d 

d 

CD 

O 

CM 

CM 

00 

CO 

-r-  CO 


p 

in 

p 

in 

p 

in 

p 

T 

T 

T 

Tt- 

00 

lU 

ill 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

rr 

CO 

o 

CD 

CO 

GO 

CD 

r-- 

o 

CM 

CD 

CO 

in 

CO 

O 

m 

in 

O 

in 

in 

O 

o 

CO 

m 

5 

CO 

CO 

o 

o 

in 

CO 

d 

* 

•y 

'y 

CO 

TT 

CO 

CD 

CO 

CO 

00 

CO 

CO 

o 

iq 

p 

CD 

iq 

CM 

CO 

p 

p 

cq 

cq 

y 

CM 

^* 

CO* 

h-* 

CD* 

ir^ 

y 

p 

p 

p 

in 

in 

p 

p 

p 

*T 

00 

CM 

LU 

LU 

LU 

LU 

LU 

UJ 

LU 

LU 

LU 

LU 

LU 

UJ 

CO 

CO 

CM 

o 

y- 

CO 

CM 

f^ 

XT 

CO 

o 

CO 

CD 

CD 

O 

r*». 

CO 

CM 

r-* 

O 

CO 

00 

in 

O 

O 

CO 

in 

CO 

pi 

in 

o 

CO 

d 

CO 

TT 

in 

CO 

CO 

▼— 

CM 

oo 

CM 

CO 

o 

cq 

P 

CD 

O 

cq 

p 

p 

ir- 

CM 

in 

r^* 

cd 

d 

y 

■r-* 

p 

p 

in 

p 

p 

in 

p 

p 

in 

in 

p 

cn 

CD 

LU 

LU 

LU 

LU 

LU 

LU 

UJ 

lU 

LU 

LU 

UJ 

in 

CO 

CO 

CM 

CO 

CO 

CD 

CD 

y“ 

in 

o 

CO 

CD 

r- 

CO 

in 

CM 

CD 

CD 

r^ 

y 

o 

CO 

CO 

CD 

CO 

r- 

in 

CM 

o 

CM 

CM 

Q 

h". 

CD 

O 

CD 

CO 

in 

CM 

CD 

cq 

in 

CD 

T-* 

CO 

CO 

t— 

CD 

CD 

CD 

y 

CM 

cd 

in 

cd 

i^* 

cd 

O* 

cd 

CD 

p 

in 

p 

in 

p 

in 

p 

m 

in 

LU 

LU 

LU 

UJ 

LU 

LU 

LU 

LU 

LU 

LU 

CM 

■M* 

CO 

CM 

CM 

CM 

CO 

o 

CO 

CO 

OO 

CO 

r^ 

CD 

CO 

CM 

o 

CO 

CD 

CO 

CD 

CO 

o 

CO 

d 

CO 

■M- 

CM 

CO 

CO 

in 

p 

p 

'Cf 

p 

p 

CM 

CO 

CM 

cd 

in 

cd 

cd 

"  M 

oc  cc  4i 

CO  ^  ^ 


T-C'4CO^m<DNCOO> 


II 

oc 

</) 


II 

OC 

79 


T-CMCO^IOIONCOO 


O  T-  CO  ^ 


intDNooooi- 

1-  T-  r-  T-  T-  CM  CM 


Inductance,  L  [Henries],  for  R=1  meter  N-Turn  Spirals  having  Various  Pitches, 


Inductance  vs.  Number  of  Turns  for  R=1  meter  Spirals  with  Various  Pitches,  S/R 
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Number  of  Turns  (Loops)  in  Planar  Spiral  Inductor,  N 
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Inductance,  L  [Henries],  for  R=1  meter  N-Turn  Solenoids  having  Various  Pitches, 


Inductance  vs.  Number  of  Turns  for  R=1  meter  Solenoids  with  Various  Pitches,  S/R 
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Solenoid  Magnetic  Field  Intensity  vs.  Radial  Position  for  S/R=0.2 

H(r)=:Hz(r)  in  plane  of  center  loop  of  N-Turn  Solenoid  of  Radius,  R,  Carrying  1=1  amp  of  Current 


r/R,  Radius  from  Axis  of  Solenoid  at  which  H-Field  Is  Measured  (normalized  to  radius  of  solenoid,  R) 


Solenoid  Magnetic  Field  Intensity  vs.  Radial  Position  for  S/R=0.2 

H(r)=Hz(r)  In  plane  of  center  loop  of  N-Turn  Solenoid  of  Radius,  R,  Carrying  1=1  amp  of  Current 
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Solenoid  Magnetic  Field  Intensity  vs.  Radial  Position  for  S/R=0.4 

H(r)=Hz(r)  In  plane  of  center  loop  of  N-Turn  Solenoid  of  Radius,  R,  Carrying  1=1  amp  of  Current 
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r/R,  Radius  from  Axis  of  Solenoid  at  which  H-Field  is  Measured  (normalized  to  radius  of  solenoid,  R) 


Solenoid  Magnetic  Field  Intensity  vs.  Radial  Position  for  S/R=0.4 

H(r)=Hz(r)  In  plane  of  center  loop  of  N-Turn  Solenoid  of  Radius,  R,  Carrying  1=1  amp  of  Current 


Solenoid  Magnetic  Field  Intensity  vs.  Radial  Position  for  S/R=0.1 

H(r)=Hz(r)  in  plane  of  center  loop  of  N-Turn  Solenoid  of  Radius,  R,  Carrying  l=1  amp  of  Current 
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r/R,  Radius  from  Axis  of  Solenoid  at  which  H-Field  is  Measured  (normalized  to  radius  of  solenoid,  R) 


Solenoid  Magnetic  Field  Intensity  vs.  Radial  Position  for  S/R=0.1 

H(r)=H2(r)  In  plane  of  center  loop  of  N-Turn  Solenoid  of  Radius,  R,  Carrying  1=1  amp  of  Current 
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Solenoid  Magnetic  Field  Intensity  vs.  Radial  Position  for  S/R=0.05 

H(r)=H2(r)  In  plane  of  center  loop  of  N-Turn  Solenoid  of  Radius,  R,  Carrying  1=1  amp  of  Current 
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Solenoid  Magnetic  Field  Intensity  vs.  Radiai  Position  for  S/R=0.05 

H(r)=H2(r)  In  plane  of  center  loop  of  N-Turn  Solenoid  of  Radius,  R,  Carrying  1=1  amp  of  Current 
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Solenoid  Magnetic  Field  Intensity  vs.  Radial  Position  for  S/R=0.03 

H(r)=Hz(r)  in  plane  of  center  loop  of  N-Turn  Solenoid  of  Radius,  R,  Carrying  1=1  amp  of  Current 
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Solenoid  Magnetic  Field  Intensity  vs.  Radial  Position  for  S/R=0.03 

H(r)=H2(r)  in  plane  of  center  loop  of  N-Turn  Solenoid  of  Radius,  R,  Carrying  1=1  amp  of  Current 


CVOOOCD'^CVJOCOCO'^CVJOOOCO'^CVJOCVjTrcOCOOCVJ'^COOOOCVj'^l'CDCOOCM 
cococmcmcmcvjcvjt-i-'t-t-t-  •  »  * 

P!OU9|os  ujni-N  lo  dooj  ja;uao  p  eu8|d  u]  ‘(j)zh  ‘Pl^y  ogauSei'^  p  ;uauodiuoo-z 


o> 

o 


CO 

o 


o 


CO 

o 


LO 

o 


o 


CO 

o 


CVJ 

o 


q: 

Ts" 

o 

c 

o 

o 

« 

o 

i/) 

2 

’■5 

a 


■o 

o 

“ra 

E 

k. 

o 

■o 

o 

L, 

3 

Vi 

cz 

o 

IE 


ay 

05 


2 

o 

E 

I 

X 


o  sz 


CO 

ay 


ay 

d 


CD 

ay 


in 

ay 


o  3 


C3 

2 

o 

c 

o 

o 

(/> 

o 

W 

*x 

< 

o 


ay 

d 

CO 

ay 

d 


w 

3 

a 

X 


Osl 

(J) 


O) 

d 


ay 

d 


w\ 


V) 


^  o  X 
T3  0  O  ^ 

C  c!)  C 

S  w  0.5 

O  O 

O  ^  ^  OQ 

\Q: 


r 


N 

X 


S 

c 

o 

O 

o 

© 


§.^ 


o 

CO 

±i  Q. 

DC  05 


jQ 


> 

'N 

\ii 


0.15Rx0.03R  Ribbon  Solenoid  Magnetic  Field  Intensity  vs.  Radial  Position  for  S/R=0.21 

H(r)=Hz(r)  in  piano  of  cenlor  0.15R  x  0.03R  ribbon  conductor  of  N-Turn  Solenoid  of  Radius,  R,  Carrying  1=1  amp  of  Current 
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0.15Rx0.03R  Ribbon  Solenoid  Magnetic  Field  Intensity  vs.  Radial  Position  for  S/R=0.21 

H(r)=Hz(r)  In  plane  of  center  0.15R  x  0.03R  ribbon  conductor  of  N-Turn  Solenoid  of  Radius,  R,  Carrying  1=1  amp  of  Current 
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Skin  Depth  and  Surface  Resistivity  vs.  Frequency  for  Copper 
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Skin  Depth  and  Surface  Resistivity  vs.  Frequency  for  Aluminum 
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Signal  Frequency,  F  (Hz) 


//AIRounclWireRsLinl_1p5mil  1,5  mil  Aluminum  Round  Wire  R.C.Eden  8/16/96 

project  X,  FRLArr,  RsArr,  LintArr,  j,  Zo,  mu,  C,  Rsij,  Lext,  Lint,  R,  Z2,  Z1oc,  freq,  En 
project  Ra,  Ro,  gam,  ZoDC,  cit,  Eo,  ceff,  Rs,  Ldc,  LlNTdc,  RCENTdc,  Rdc; 
local  i,  y,  SqtEr,  Rsel,  w,  Sqrtw,  zzz,  VZt; 


//  Calculational  Parameters 

FreqStep=0.25E8;  /*  FreqStep  for  calculation*/ 

NS=1 01 :  /*  Number  of  frequency  steps  in  calculation  */ 

Inf=1e14; 

//  Physical  Constants:  All  units  are  meters,  seconds,  ohms,  farads,  henrys.  Hertz,  etc. 

pi=3.1 41 59265358979323846; 

j=(0,1): 

rSqrt2=1/sqrt(2): 

mu=(4E-7)*pi; 

clt=2.99792458E8:  /*  Speed  of  Ught  in  m/s  */ 

Eo=1/{mu*clt*clt);  /*  Permitrvity  of  free  space  */ 

Er=3.90:  r  Relative  dielectric  const  of  coax  insulator  */ 


SqtEr=Sqrt(Er); 

ceff=clt/SqtEr;  /*  Speed  of  light  in  the  dielectric  */ 

rho=2.62E-8;  /*  Resistivity  of  Aluminum  at  20°C  in  n-m  */ 


//  Line  Dimensions  and  dc  impedance,  ZoDC  (ac  Zo  is  complex  for  lossy  line) 
Ra=1.27E-5*1.5;  /*  Center  Conductor,  1.500  mil  dia  */ 

RCENTdc=rho/(pi*Ra*Ra):  t*  Coax  Center  Conductor  resistance/meter  */ 

Rdc=(rho/pi)*((1/(Ra*Ra))):  /*  Actual  dc  resistance  as  there  is  no  shield  */ 
LINTdc=(miy(8’pi)):  /*  Lo'  value  for  Gardiol's  Eq.  4.71  incl  shield  */ 

x=1.0000;  /*  x=1  meter  interconnect  line  length  */ 


//  Function  Definitions 

hh=0.02;  /*  Derivative  step  size  used  when  these  are  called  */ 

function  Fctl(x) 

return  bei(0,x): 

end  function; 
function  FctR(x) 

return  ber(0,x); 

end  function; 


//  Begin  Actual  Frequency  Domain  Calculation: 
i=1 ;  r  freq=0  dc  case  prior  to  frequency  sweep  ’/ 

freq[i]=FreqStep*(i-1); 
w=2*pi*freq[i]; 

RsArr[i]=Rdc; 

LintArr(i]=LINTdc; 

FRLArr{i,1]=FreqSfep*(i-1): 

FRLArr[i.21=Rdc; 

FRUrrti,31=LINTdc; 


//  Frequency  Sweep: 
for  i=2  to  NS  do 

freqp]=FreqStep*(i-1);  /*  Save  frequency  values  array  for  use  with  MagY(0  plots*/ 
w=2*pi*freq[i]; 

//  Skin  effect  depth  delta=sqrt((2*rho)/(w*mu));  and  zeta=(Ra*sqrt(2.0))/delta; 
zeta=Ra*sqrt((w*mu)/rho):  r  if  zeta>10,  the  high  frequency  approx  is  good  */ 
if  zeta>10then 

Rac=Rdc*(zeta/2)*(rSqrt2+(0.501 1275/(zeta-0.475)));  /*  Approx  RacBer  Fn  */ 
zzz=0.51 56/(zeta-0.667);  /*  temporary  var  in  numerical  Approx  to  LintBer  */ 
Lint=LINTdc*(4/zeta)*(rSqrt2-(zzz*2zz));  /*  Approx  to  UntBer  Fn  */ 

end  if; 

if  zela<=10then 

BUR=ber{0,zeta); 

Bl=bei(0,zefa); 

BURP=derivative(FctR,3,hh,zefa); 

BIP=derivafive(Fctl,3,hh,zeta); 

Rac=Rdc*(zefa/2)*((BUR*BlP)-(BI*BURP))/{(BIP*BIP)+(BURP*BURP)); 

Lint=LINTdc*(4/zefa)*((BI*BIP)+(BUR*BURP))/((BIP*BIP)+(BURP*BURP)); 

end  if; 

LintArr(i]=Lint; 

RsArrti]=Rac; 

FRLAnii,1]=FreqSfep*{i-1); 

FRLArrti,2]=Rac; 

FRLArr{i,3]=Lint; 
end  for; 


O.OE+0  5.0E+7  1.0E+8  1.5E+8  2.0E'f8  2.5E+8  3.0E+8  3.5E+8  4.0E+8  4.5E+8  5.0E+8  5.5E+8  6.0E+8  6.5E+8  7.0E+8  7.5E+8  8.0E+8  8.5E+8  9.0E+8  9.5E+8  1.0E+9 

Signal  Frequency,  f  (Hz) 


ac  Resistance,  Rac/I,  and  Internal  Inductance,  Lint/I,  for  1.5  mil  Diameter  Al,  Au  and  Cu  Wires 


Signal  Frequency,  f  (Hz) 


ac  Resistance,  Rac/I,  and  Internal  Inductance,  Lint/I,  for  1.0  mil  Diameter  Al,  Au  and  Cu  Wires 


O.OE+0  1.0E+0  2,0Ef8  3.0E+8  4.0E+8  5.0E+8  6.0E+8  7.0E+8  8.0E+8  9.0E+8  1.0E+9  1.1E+9  1.2E49  1.3E+9  1.4E+9  1.5E+9  1.6E+9  1.7E+9  l.OE+9  1.9E+9  2.0E+9  2.1E+9  2.2E+9  2.3E+9  2.4E+9  2.! 


1.5  mil  Aluminum  Wire  Round  Coil  Inductor  Q  vs.  f  for  various  numbers  of  turns, 


Frequency  Dependence  of  ac  Sheet  Resistance  of  Aluminum  Plates  with  One  or  Both  Sides  Conducting 
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Inductance  vs.  Number  of  Turns  for  1.5  mil  Al  Wirebond  Inductor 
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Number  of  Complete  (Wirebond  +  Metal  Stripe)  Turns  in  Solenoid  Inductor, 


MCM  1.5  mil  Al  Wirebond  Inductor  Q  vs.  f  for  various  numbers  of  turns, 
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Number  of  Turns,  N,  in  MCM  Wirebond  Solenoid  Inductor 


Q  at  f=950  MHz  of  1 .5  mil  Al  Wirebond  Solenoid  Inductors  vs.  Number  of  Turns, 
for  Various  Thicknesses,  tm,  of  the  100  |j.m  Wide  Substrate  Metallization  Path  . 
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Simplified  Inductor  Model  and  Method  for  Calculating  True 
Inductance,  Internal  Capacitance,  Self-Resonant  Frequency  and  Q 


